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r attempting to speak in honor of Dr. Henry 
Crew I thought of trying to find a subject 
about which I might perhaps know nearly as 
much as he. However, I soon gave this up as 
impossible, and decided on a short survey of the 
semiastronomical subject of the light of the 
night sky. Dr. Crew was astronomer at the Lick 
Observatory from 1891 to 1892 and for many 
years was an editor of the Astrophysical Journal. 
Therefore, he is familiar with the night sky and 
will be appreciative of the efforts of the next two 
generations to learn about it. 

In places remote from artificial illumination at 
night with no moon the luminosity of the sky 
is due to several sources of light, all of which 
are at a considerable distance. The sources are 
(a) stars and nebulae in interstellar space; 
(b) zodiacal light, comets, possibly scattered 
sunlight (if there is something there to scatter 
the sunlight) in interplanetary space; and (c) 
emissions from the gases of the upper atmosphere 
and polar aurorae in the upper atmosphere. 
Omitting polar aurora from consideration, the 
high atmospheric emissions are the strongest 
source of night illumination, being four or five 
times as intense as all of the other sources com- 
bined. 


In more detail, the intensity is divided as 


























































































* Dr. Henry Crew, Senior Past President of the American 
Physical Society, was guest of honor at a dinner of the 
Society on May 17, 1949 during its Semi-Centennial 
Meeting in Cambridge, Mass. Dr. E. O. Hulburt’s paper 
was presented on the morning following. 
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follows: for the pinongpeghte penny region 
3700-4500A, about % due to the stars, $ due to 
the high sniniiinie emissions, and possibly ¢ 
from interplanetary space (the amount attributed 
to sources in interplanetary space has decreased 
as the measurements have improved; it is now 
down to about ¢ of the whole and may eventually 
become a few hundredths or even disappear) ; for 
the visible —r seen with the dark-adapted 
eye, about } due to the stars.and ¢ due to the 
high stemiilions. 


Spectrum 


The night sky light is feeble and even with 
spectrographs of the highest light-gathering 
power relatively low dispersion and long ex- 
posures are necessary to obtain spectra. We need 
to refer only to recent work and mention that 
the best spectra at present available are those 
obtained by Barbier* in 1942-44 at the Observa- 
tory of Haute Provence, France, and by Elvey, 
Swings, and Linke‘ in 1939 at the McDonald 
Observatory, University of Texas. 

The dispersion of the spectrograph used by 
Barbier was 150 and 630A mm at 3200 and 
5000A, and the length of the spectrum from 3200 
to 5000A was about 7 mm; exposures from 100 
to 200 hours were used. The McDonald spectro- 

1D. Barbier, ig Rendus 224, 635 (1947). 

2E. O. Hulburt, J. Opt. Soc. Am. 39, 211-215 (1949). 


3 D. Barbier, Ann. d. Geophysik 1, 224-232 (1945). 


4 —— Swings, and Linke, Astrophysical J. 93, 337-348 
(1941) 
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graph had slightly higher dispersion and much 
greater light-gathering power; the dispersion 
was 115 and 500A mm at 3200 and 5000A, the 
length of the spectrum from 3200 to 5000A was 
about 10 mm, and exposures of 9 hours were used. 

The night sky spectrum is very complex and 
is composed of many lines and bands, for the 
most part incompletely resolved. In Fig. 1 is 
shown a spectrogram of the night sky taken with 
the McDonald spectrograph.‘ For wavelengths 
below 5000A the intensity is rather uniform with 
no emissions of outstanding strength. Above 
5000A the green line 5577, the yellow pair 5890, 
5896, and the red group at 6300, 6364, are con- 
spicuous features on the photographic spectrum. 
The dispersion of present spectra in the red and 
near infra-red is very low and the region is very 


Fic. 1. The spectrum of the night sky, McDonald 
Observatory (Ref. 4). 
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poorly resolved. A. B. Meinel® of the Lick Ob- 
servatory is beginning to obtain better data for 
wavelengths greater than 6000 by the use of a 
7500 line transmission grating and an f/1 camera, 
with a dispersion of 180A mm7—'; he has published 
only a few preliminary results. 

The main features of the spectrum of the 
night sky light, which are considered certain, are: 


Atomic lines. The atomic spectrum con- 
sists of the lines 5577, 6363.9, and 6300.2, 
forbidden transitions of oxygen, with 5896 
and 5890, sodium. 

Molecular lines and bands. Among the 
most intense are 10,440 and 6560, ist posi- 
tive system of Ne; the 10,440 emission is 
by far the strongest emission thus far ob- 
served, being about 30 times as intense as 
the green line 5577, the next most intense 
emission. It was discovered in 1940 by 
Stebbins, Whitford, and Swings,* by means 
of a photoelectric cell and filters, who de- 
termined the wavelength within +25A, and 
independently by Herman, Herman, and 
Gauzit.” 

The strongest blue-violet bands consist of 
the Vegard-Kaplan system of Nz; about 35 
bands are identified. 

The strongest ultraviolet bands are the 
Herzberg system of Oz, about 20 bands are 
identified. 

Some weaker bands occur in violet and 
ultraviolet regions of the spectrum due to 
unknown systems. 

A number of unidentified lines or bands, 
some partially resolved, occur from 5000 to 
10,000A. 

A sort of continuum and unresolved bands 
are found, which are fairly weak in the ultra- 
violet and which increase in intensity for 
wavelengths longer than 4000A. 

A weak continuous background with trace 
of Fraunhofer absorptions which was men- 
tioned in earlier work was not observed in 
the better more recent spectrograms. 


5 A.B. Meinel, Pub. Ast. Soc. Pacific 60, 373-377 (1948). 

6 Stebbins, Whitford, and Swings, Astrophysical J. 101, 
39-46 (1945). 

™ Herman, Herman, and Gauzit, J. de Physique, Ser. 8 
6, 182 (1945). 
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The Schumann-Runge system of O2 and the 
Lyman-Birge system of Nz have been looked for, 
but no certain identifications have been made; 
similarly for CO and COs, ozone, oxides of 
nitrogen, NH and N;, CH, CN, OH, H.0O, and 
compounds of Na, Si, and S. 

The spectral intensity distribution received at 
the surface of the earth is approximately as 
follows:® 





Wavelength Emitter Intensity 


10440A Ne 190 X 10-4 erg cm sec™ 
6300, 6363 O 5x10-4 
5890, 5896 Na 0.7 10-4 
5577 O 7X10 
Main band systems Na, O2 10X10-. 


The only reported spectral energy distribution 
of the night sky light is that of Babcock and 
Johnson® derived from measurements of a small 
photographic spectrum of dispersion 1100A mm7 
at H,. Their curve is given in Fig. 2; it rises with 
increasing wavelength throughout the extent of 
the spectrum from 3800 to 6500A, with bumps at 
the various lines and bands. A color temperature 
of 3450°K was ascribed to it. Thus the night 
sky light may be said to be yellowish in color. 


Variations 


Excluding polar aurorae the variations in the 
total intensity of night sky emission are relatively 
small, rarely as much as a factor of 2. There 
appear to be during the full nighttime hours 
many small irregular and complex intensity 
variations which vary with the place in the sky 
and which are different for the different wave- 
lengths. There are strong intensity changes 
during twilight in certain lines. 

Diurnal variations. Dufay and Mao-Lin! at 
Haute Provence Observatory obtained 588 spec- 
trograms during 189 nights from October 1940 
to January 1944. From these spectrograms sys- 
tematic measurements were made of the intensity 
variations of the oxygen lines 5577 and 6300A 
and of the sodium pair 5893. The 5577 line 
usually showed a weak maximum around mid- 
night (not more than 40 percent) often obscured 


8D. R. Bates, Gassiot Committee Report, Physical Soc. 
London, 21-33 (1948), 

®*H. W. Babcock and J. J. Johnson, Astrophysical J. 94, 
271-275 (1945). 


10 J. Dufay and T. Mao-Lin, Ann. d. Geaphysik 2, 189- 
230 (1946). 
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Fic. 2. Spectral distribution of night sky light (Ref. 9). 


by fluctuations; the 6300 and 5893 lines merely 
weakened slowly during the night, the enfeeble- 
ment being greater for 6300 than that for 5893 
which only amounted to 10 or 15 percent. 

Twilight effects. The green line 5577 was ob- 
served to exhibit no twilight enhancement. On 
the other hand the sodium yellow lines 5893 
dropped suddenly in intensity to 1 percent of 
their former value at about 10 minutes after 
sunset" with a corresponding recrudescence near 
dawn. Likewise the red oxygen lines 6300 showed 
a similar but slower change during twilight.” 
The strong N.* flash at twilight discovered by 
Slipher™ is intense in aurorae but absent from 
the night sky. 

Annual variations. In a temperate latitude 
5577 has slight maxima in October and February, 
and a 27-day variation in intensity, hence a 
faint correlation with solar phenoména; 6300 
and 5892 have maxima near midwinter and 
minima in summer of amplitude factors of about 
2.5 and 5, respectively; they seem indifferent to 
solar phenomena.!° 

Latitude variations. Measurements with cali- 
brated and standardized visual photometers indi- 
cated no changes, which could not be attributed 
to variations in haze, in the visual brightness 


1 R. Bernard, Zeits. f. Physik. 110, 291-302 (1938). 
2 C, T. Elvey and A. H. Farnsworth, Astrophysical J. 96, 
451-467 (1942). 


1%V. M. Slipher, Roy. Astr. Soc. M.N. 93, 657-668 
(1933). 
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of the night sky with latitude when no polar 
aurorae were visible.” “ Figure 3 is a photograph 
of one of the visual photometers. The observa- 
tions were made at Bocaiuva (Brazil), Bikini, 
Maryland, Whiteface Mountain (New York 
State), and Greenland, at latitudes —17°, +12°, 
+39°, +45°, and +63°, respectively. Average 
data for these stations are plotted in Fig. 4. 
Some older visual measurements" of the sky 
near Polaris indicated that the brightness in- 
creased by a factor of 4 from latitudes +44° to 
+80°, but this result may have been influenced 
by polar aurorae. 

Visual measurements of the night sky lumi- 
nosity refer mainly to the strong oxygen atomic 
line 5577. Because, by multiplying the night sky 
spectral energy curve by the sensitivity curve of 
the dark-adapted eye, one finds that 5577 has 
the main effect (38;) with only minor influence of 
the oxygen red lines 6300, 6363 the sodium yellow 


Fic. 3. Low brightness visual photometer, Naval 
Research Laboratory (Ref. 2). 


4 Data for latitude +63° were observed in February 
and March, 1949, by personnel at an Air Force Weather 
station in Greenland, using an NRL visual low brightness 
photometer. 

6B. Fessenkoff, C.R. (Doklady) Acad. Sci. URSS 32, 
320-322 (1941). 
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lines 5896, 5890 and the molecular bands in the 
blue and green. 

Variation with zenith angle. The variation of 
the visual brightness of the night sky with zenith 
angle was found to be about the same for lati- 
tudes from —17° to +63°, and increased from 
about 130 millimicrolambert (muL; 1 muL =10-° 
lambert = 2.96 X10~? candle per square foot) at 
the zenith to a maximum of about 250 muL at 
about 15° above the horizon for a clear atmos- 
phere; with slight haze the maximum near the 
horizon became less pronounced, and with more 
haze disappeared.? The variation was observed 
to be somewhat different for red and green 
wavelengths,'* but has not been thoroughly ex- 
plored for all wavelengths. 


Altitude 


The altitude of the high atmospheric emissions 
has been determined from the intensity variation 
with zenith angle, and erratic and uncertain 
values from 60 km to several hundred km were 
found. The uncertainties were due to nonuni- 
formity of the emissions and to incompletely 
worked-out corrections for atmospheric attenu- 
ations. But even with better correction formulas 
no improved or more trustworthy values of the 
altitudes were determined ;? and in no case was 
the atmospheric attenuation measured at the 
same time that the sky intensity was observed. 
Swings!’ stated the situation thus: ‘There are 
great irregularities in the distribution of the 
emission over the sky. These irregularities ap- 
pear consistently when simultaneous exposures 
in different azimuths at the same zenith distance 
are compared. No layer of uniform brightness 
exists, but rather a set of bright clouds which 
move around and change in brightness. The only 
remaining hope is that, by taking a sufficiently 
large number of observations, the erratic fluctua- 
tions will average out.” 

The altitude of the yellow emissions 5890, 5896 
of sodium was determined in another way based 
on the fact that these radiations decreased 
suddenly in intensity at twilight when the region 
in which they originated was shielded from the 


16 Abadie, Vassy, and Vassy, Ann. d. Geophysik 1, 189- 
224 (1944). 


17P, Swings, The atmospheres of the earth and planets 
(University of Chicago Press, 1948), Chapter 6. 
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direct ultraviolet rays of the sun. The calcula- 
tions were made by Penndorf!* from the obser- 
vations of Vegard and Tonsberg!® in February 
and March, 1949, at Tromso and Oslo. The 
altitude of the sodium emission at twilight was 
determined unmistakably to be in the region 
from 80 to 116 km. At present it seems reasonable 
to keep in mind the probability that the various 
night sky radiations may originate in different 
levels which may change in the course of the 
night and which may vary with the latitude 
and season. 


Theory 


The spectral identifications show that most of 
the night sky emissions come from the two most 
abundant gases of the atmosphere, oxygen and 
nitrogen; a small portion comes from sodium. 
No one has. questioned the general notion that 
the nocturnal emissions derive their energy from 
the sun, but no processes or quantitative details 
have been agreed upon because information is 
lacking about radiation transitions and atmos- 
pheric composition. In fact, without exception, 
all theoretical processes, qualitative or quantita- 
tive, thus far proposed, have been weighted 
down with a wealth of criticism.* 

The excitation energy of the line systems 
5577 OI, 6300 OI, 5893 Nal are 4.2, 2.0, and 
2.1 ev, respectively; and of the band systems 
Vegard-Kaplan Ne, Herzberg Oz, Schumann- 
Runge Oz, and Lyman N; are 7.0, 4.7, 6.2, and 
8.7 ev, respectively. It has been suggested that 
the energy is provided by (a) the association 
energies of atomic O and N, which are 5.09 and 
7.38 ev, respectively; (b) the energies of ioniza- 
tion, which are about 13.5 and 14.5 ev for the 
first ionization potentials of O and N, respec- 
tively; (c) particles of solar origin. It appeared 
that (a) was sufficient for the line systems but 
for none of the band systems except Herzberg; 


18R, Penndorf, J. of Meteor. 5, 152-160 (1948). 


19],, Vegard and E. Tonsberg, Geophysical Publ. 13, 
1-22 (1940). 
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Fic. 4. Night sky brightness against zenith angle 
at several stations. 


(b) seemed energetically satisfactory, but not 
completely in accord with the Frank-Condon 
probabilities ; (c) is speculative and has not been 
put in quantitative terms. 

The relatively small amount of sodium in the 
upper atmosphere necessary to account for the 
5893 emissions may come either from the surface 
of the earth, as the sea, or from interplanetary 
space. Both views have been suggested and the 
correct one is not known. 

In conclusion, we have attempted to exclude 
the aurora from the foregoing survey of the 
night sky light knowing that this cannot be done, 
and in the end should not be done. It is probable 
that a complete understanding of the night sky 
cannot be reached without at the same time 
having a complete understanding of the aurora. 
Both, it is believed, derive their primal energy 
from the sun. The auroral emissions are from 
the same familiar gases of our atmosphere, 
oxygen and nitrogen, as are those of the night 
sky, but with different spectral distributions, 
with enormously greater intensities and in quite 
different geometrical patterns. The auroral energy 
is received in concentrated form from explosions 
or outbursts on the sun, in contrast to the 
relatively gentle and steady solar energy stream 
which supports the normal night sky light. The 
kinship between the aurora and the night sky 
light may be analogous to that between a 
hurricane and a gentle rain. 


History is simply a piece of paper covered with print; the main thing is still to make history, not 


to write it—BISMARCK. 





On the Mechanism of Brownian Motion in Liquids 


FREDERIK J. BELINFANTE 
Purdue University, Lafayette, Indiana 


1. Einstein’s Law 


HE following derivation of the well-known 

formula for the mean square path traveled 
in Brownian motion by a spherical particle sus- 
pended in a liquid has the advantage of being 
more illuminating than the usual derivations.! In 
fact, it is based upon the qualitative discussion 
used in many introductory courses of modern 
physics, here made quantitative by a few simple 
formulas. 

Let a spherical particle of radius r, mass m, be 
suspended in a liquid of viscosity 7. Occasionally 
a molecule may hit the particle with exceptional 
vigor and set the particle into motion through the 
liquid. As a consequence, the probability of head- 
on collisions against other molecules is slightly 
increased, while the probability of a collision by a 
molecule overtaking the particle is slightly de- 
creased. Also, the particle will, in subsequent 
collisions, give back some of its gained mo- 
mentum to the newly colliding molecules.” Alto- 
gether, there are correlations between the effects 
of preceding collisions and of subsequent colli- 
sions. These correlations favor a retardation of 
the particle in its motion, an effect known as 
friction. From the above it is easily seen that in 
first approximation friction may be described as a 
force proportional to the velocity of the particle, 
so that 


mdv /dt = — Sv. (1) 


For big spheres one may use Stokes’ law, putting 
S=6rnr, (2) 
but it is well known that this law breaks down for 


1A derivation showing similarity to the one given here 
can be found in Miiller-Pouillet’s Lehrbuch der Physik, vol. 
III/2. 11. Auflage, pp. 327-328. (F. Vieweg & Sohn, 
Braunschweig, 1925.) There, however, the relaxation time 
t= m/S replaces the long time / used in our derivation, and 
as a consequence in Eq. (353), loc. cit., a factor remains 
unexplained, which does not appear in our derivation. 
Moreover, relaxation time is assumed there to be of the 
order of magnitude of the time between two collisions 
effectively disturbing the rectilinearity of the path of the 
rticle, quite contrary to our result expressed immediately 
ese our Eq. (24). 
I wish to thank Dr. H. M. James for valuable discus- 
sions about these and some other statistical questions. 


spheres of very small radius.* As we shall see 
below, it breaks down sooner in liquids of high 
viscosity than in liquids of low viscosity. 

Having taken into account the correlations be- 
tween subsequent collisions in this way, we may 
further consider them as independent of each 
other. At random times ¢; they will give the 
particle random impulses mV;. Due to these 
random forces as well as to the friction, the total 
velocity v of the particle at a time ¢ will be the 
vector sum of the velocities v; left over from V; 
after an interval (t—t;). Hence 


V(t) = Dia Vilé); (3) 


where v;,(¢) denotes the x-component of v;(¢) and 
where the #’s between brackets indicate that the 
velocities are functions of t., The summation is to 
be extended over all collisions (numbered by the 
label z) prior to the time ¢. 

When we calculate the average of the square 
[ v(t) ]* of the velocity at a time ¢, and we take 
v(t) from the first part of Eq. (3), then all terms 
(v;°v,;) for «#7 can be omitted, being equally 
often negative as positive. Taking also the time 
average of this mean square velocity over some 
long time interval from 0 to ¢o, we find, say for the 
mean square of the x-component 


Viz(t) = Vise 


—(S/m)( t—ti) Y 


Lw 


(0< ti< to) 


{Lve(t) }?} av =to f [v..(t) ]2dt 


+ Law 


(ti<0) 


f j [v:2(t) ]*dt , (4) 
0 


We substitute here v;2(¢) from the second part 
of Eq. (3), and perform the integrations. As- 
suming that there are f collisions per second 
giving appreciable impulses to the particle, we 
may further replace the summation signs by their 


3 Miiller-Pouillet, see reference 1, pp. 83-89. 

4 That is, we describe the forces exerted by molecules 
against a moving sphere as the superposition of the force of 
friction described by Eq. (1) and independent completely 
random forces giving the particle sudden impulses at 
irregular intervals. 
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equivalents 


to 
Li =f f dt;; 

(0< ti< to) 0 
Performing these last integrations too, while we 
replace V;,2 by its average, we obtain 
{Loz(t) }?} av = fm{ V iz”) av/2S, 


a result which we might also have obtained with 
less integrating but more thinking’ from 


La -ff dt; (5) 


(ti<0) 


(6) 


C(O ")av= (Ito) 4 f | [ott Pa. (7) 


5 If vi( —nto)x(t) is a symbolic notation for the x-compo- 
nent of the velocity at ¢ left over from the velocity acquired 
by a collision approximately at the time (¢;—mto) (that is, 
from a collision that took place about mto seconds before 
collision number 7), then in Eq. (4) we might replace the 
second sum (over ¢; <0) by 


oo to 
2 2 f [vi( —nto)x(t) dt, 
0 


(0 <ti <to) ™=! 


since obviously for each collision taking place before t=0 
there will be some other collision taking place between t=0 
and ¢=¢o just about a multiple of fo seconds later. Now, for 
purposes of calculating the mean square, we may write 
Viz(t-+-nto) instead of vi(—nto)x(t), because the mean square 
of v;2(t’) will depend on the time difference (¢’—#;) only. 
Also, we may write 


- to ie (x +1)to 
2 f [viz (t-+-nto) Pdt= 2 f viz(t’)2dt’ 
n=l. 0 n=1 


nto 
= f 0%2(t’) dt’. 
to 


Combining the result thus obtained with the first sum in 
Eq. (4) we find that 


{Lvz(t) By av=to? Zip f [vsz(t) }2dt. 


0 <ti <to 


As the integral over ¢ here can be written as an integral over 
(t—t;) and then becomes clearly independent of 7, we may 
now replace the summation over 7 simply by a multiplica- 
tion by the total number of collisions ft in the time interval 
from 0 to éo, and thus obtain the result of Eq. (7). Thus, the 
extension of the integration over ¢ in Eq. (7) to © (instead 
of to to as in Eq. (4)) here accurately takes care of the after- 
effects from all the collisions at times before =0, which in 
Eq. (4) were represented by the second sum, but which in 
Eq. (7) do not appear explicitly. For the rest, instead of 
taking these after-effects of all those collisions in the far 
past accurately into account, we might as well have 
neglected them completely, as well as we might have 
written integrals from ¢; to © for the original integrals from 
t; to to simply by neglecting the difference between these 
integrals, at least, if we choose to>>m/S. Thus, we might as 
well replace the words “more thinking’’ in the text by 
“more neglecting,’”’ and then have by Eq. (7) a simple 
means of obtaining the result shown in Eq. (6), by per- 
forming the one integration over ¢ in Eq. (7) after substitu- 
tion of v;2(t) from Eq. (3) into it. : 
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The law of equipartition of energy for the one 
degree of freedom vz now gives 


$m {[v.(t) ]*}av=3kT, (8) 
(and similarly for v, and v.), so that, by Eq. (6) 
fm? V iz*) ay = 2SRT. (9) 


The total displacement x(to) of the particle 
during the time interval from ¢=0 to t=? is again 
a vector sum 


X(to) = Deciy Ki(to) ; 


(ti< to) 


Xi (to) -{ vi(t)dt. (10) 


In the integral, we put v;(¢t) =0 for ¢<¢;. We shall 
calculate the mean square of the x-component of 
x(to). Again the cross terms (x; * x;) have a mean 
value zero, and we find the following result 
showing some analogy with Eq. (4): 


(lt) a= Doo | f ou(at} 


+E f cutbat} (11) 
0 


(ti<0) 


Into this expression, we substitute v;2(¢) from 
Eq. (3), integrate over ¢ and replace V;.? by its 
average. Next, we make use of Eq. (5). Per- 
forming the integrations over ¢; and putting 
m/S=r, we find 


{Lx(to) FY av 


= fm Vatad* {to— r(1 —e-tolr) } “ (1 2) 


Assuming that fo is chosen to be large in compari- 
son with the relaxation time 7=(m/S), we find 


{ [x(to) ]?} av = ftom Viz?) av/S*, (13) 


a result which we might again have obtained in 
an easier way by more thinking and neglecting® 
and less integrating, by means of 


(L(t) a (ft) | f | ou(dat} . (14) 


6 If r<to, the motion of the particle due to momentum 
acquired prior to ¢=0 is negligible during the time interval 
from t=0 to t=to, while the displacements x;(éo) due to the 
collisions during this interval are practically equal to x;(~ ) 
with the exception of the displacements due to only very 
few of these (fto) collisions near the very end of this long 
interval. Thus 2i)(xi(to)*) av leads to Eq. (14). 
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TABLE I. Critical constants of liquids. 


Liquid Chemical formula 
CH;COCH; 

CsH;NH: 

CceHe 

CCl, 

CsH1,0OH 
CH;COOC:H; 
CH;CH,0H 
C:H,OC2H; 
CH:OHCHOHCH:0H 


Hg 
CH:Br:. 


Acetone 

Aniline 

Benzene 

Carbon tetrachloride 

Cyclohexanol 

Ethyl acetate 

Ethyl alcohol 

Ethyl ether 

Glycerin 

Mercury 

Methylene bromide 

Phenol CsH,OH 

Propylacetate CH,COOC;H; 
ater H:.O 


Now substituting fm?(V;.")av from Eq. (9) into 
Eq. (13) we obtain Einstein’s law for the mean 
square displacement in the x direction: 


{[(to) }?/to} av =2kT/S (=kT/3xnr). (15) 


The expression between brackets is valid, if 
Eq. (2) is true. We may remark that both f and 
(Viz) ay drop out in reaching Eq. (15). This makes 
the result independent of special assumptions 
about the number of collisions and their direct 
effect upon the particle, as long as we can believe 
the equipartition law, Eq. (8). 


2. Discussion 


It is now interesting to estimate at least the 
orders of magnitude of the impulses mV; and of 
their frequency f. An upper limit for the latter is 
given by the number of thermal collisions per 


second of molecules against the spherical particle, 
which is 


F=4nr?-2(c)ay:N’ =2N'r(2ekT/yu)'. (16) 


Here, (c)ay is the average thermal velocity of 
molecules in the liquid; the factor } comes in 
through averaging over angles. N’ is a kind of 
“effective” density of the liquid molecules. If the 
liquid were a rarified gas, N’ would be simply the 
number of molecules per cm? (= JV). In our case, 
N’ is increased due to the crowdedness of the 
molecules in a liquid, by a factor C, which may be 
estimated’ as being of the order of magnitude 


7 The occurrence of the factor (U—bd) instead of VU in 
Van der Waals’ equation is roughly explained by stating 
that, because of the finite dimensions of the molecules, they 
do not have a volume V at their disposal for free motion, 
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U/(U—b), where VU is the molal volume, while } 
according to Van der Waals’ theory of fluids 
would be one-third of the critical (molal) volume. 
Introducing specific gravities p instead of molal 
volumes by p= M/1, we find by this first method 
of estimating C 


Sperit 
C..=———, 
SPecrit— p 
thence 


SpPerit 


N/=C,N= (17:1) 


#(3pcrit—p) 
Here, u is the mass of a molecule of the liquid. 
A second method of estimating C is given by 


Us 
G-——, 
3(C) av 


pus ( TH ) 
w \2RT] © 
This guess is based on the following ideas: Ac- 
cording to the well-known formula for the velocity 
of sound v,=(yp/p)! (with y=c,/c, for an ideal 


gas and y ~1 for a liquid), this velocity would be 
a constant independent of the frequency »v or the 


(17:2) 


but a volume that is smaller due to the fact that each 
molecule has to stay with its center at a distance of at least 
twice the molecular radius from the centers of other 
molecules. As a consequence, the molecules hit the walls 
more frequently and thus exert a larger pressure, which is 
now proportional to 1/(0—b) instead of to 1/VU, thus 
being increased by a factor C,= U/(U—d). Exactly the 
same reasoning may be considered applicable to the 
calculation of the number of collisions against a suspended 
particle. 
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wavelength A of the vibrations. From the deriva- 
tion of this formula it is clear that v, is meant 
here to be the phase velocity of the sound waves, 
so that v=v,/d. The group velocity of such sound 
waves then can be calculated from this by the 
well-known formula v,=dv/d(1/d). But since v, 
was constant in our equation for », this gives 
simply v,=v,. On the other hand, the group 
velocity is the velocity with which a short signal 
would be propagated through the liquid. As the 
wave front of such a signal moves forward 
through the liquid by the motion of molecules, 
one would then expect this velocity to be of the 
same order of magnitude as the average com- 
ponent (in the direction of propagation of the 
wave) of the thermal velocity c of the molecules. 
That is, one expects v,~3(c)ay. For an ideal gas 
this is actually the case (with v, = (ya/2)!-3(c)av). 
If the molecules are very crowded, however, the 
velocity of sound can become much larger than 
the thermal velocities; as each time, when a 
molecule at the front of a sound wave collides 
against another molecule in front of it, the,wave 
motion jumps forward over a distance up to 
twice the radius of a molecule. This causes an 
increase of the velocity of sound by the same 
factor, by which also the frequency of collisions 
of molecules against the walls is increased by the 
crowdedness of the molecules. (The reasoning in 
both cases is almost identical.) The latter ratio 
was just C,=U/(U—)), so that one might ex- 
pect now an approximate equality between 
C,=V/(U—)b) and Co=9,/3(C)av. 

It must be granted that in practice the agree- 
ment between these two expressions for various 
liquids is far from striking, as is shown in Table II, 
which has been calculated with use of some of the 
data listed in Table I. Therefore, in the following, 
both expressions (17:1) and (17:2) are used from 
time to time (compare for instance Eq. (26:11) 
with Eq. (26:12)). For many liquids, it is simpler 
to use Eq. (17:2) than to use Eq. (17:1) for often 
the velocity of sound is known, whereas the 
critical density (from which } is estimated) is not 
known. 

A third estimate of N’ can be made on the 
grounds that the frequency F of thermal colli- 
sions also determines the internal pressure in a 
liquid, so that 
Pint = N’RT. 
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TABLE II. Crowdedness factor for liquids. 


Carbon 
tetra- Ethyl Ethyl 
Liquid Acetone Benzene chloride acetate alcohol Water 
Ci 67 2.23 21 40 22.2 6 
C: 7.6 9.2 9.2 9.0 6.5 4.93 


By Van der Waals’ theory, with U~d for a 
liquid, the internal pressure Pine=p+(a/V*) 
=a/V? ~a/b?=27 Perit, SO that 


N3! =27derit/RT. (17:3) 


Effective molecular densities as calculated by 
these three completely different methods do not 
differ among each other more than by factors 
smaller than 100. For various liquids at room 
temperature, they all give for (F/r?), by Eq. (16), 
values between 107? and 107° sec~! cm-*. (See 
Table III for the results obtained by use of Eqs. 
(17:2) and (17:3) respectively.) 

If such thermal collisions against the particle 
occurred one at a time, one could calculate the 
impulse by the law of conservation of momentum. 
The impulse mV, which the particle would receive 
by such an unhindered collision, would be given 
by mV = 2uu, where u is the vector component of 
the thermal molecular velocity c normal to the 
surface of the particle. Averaging the squares of 
these impulses as the collisions take place,’ we 
have 


mXV?) ay = 4u2(U2) av. 


In the calculation of (u?),v it should be remem- 
bered here that the frequency of collision itself is 
proportional to u, so that 


(w= (f ensrariean) / * 
0 
( f rmsintudu ) = 2BT /p. 
0 


mV?) av = SukT ’ m V.) w= Suk T/3. (18) 


The momentum acquired by the particle by 
one such collision is rather small and much less 
than the total momentum mv of the particle re- 
quired for validity of the equipartition law, Eq. 
(8), since the mass » of one molecule is much 
smaller than the mass m of the suspended par- 


Thence 
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ticle. Yet, these collisions will be rather effective 
in causing a jerky motion of the particle, due to 
their high frequency. This frequency, however, is 
in fact not given by. F as calculated by Eqs. 
(16)-(17), but is much lower. That the frequency 
of such one-at-a-time collisions, at least for 
sufficiently big particles, cannot possibly be given 
by Eq. (16), is easily seen as follows: In Sec. 1, 
we derived Eq. (9) from the law of equipartition 
of energy, Eq. (8), together with the law of fric- 
tion, Eq. (1), without specifying at all what the 
value of S might be.’ Now, if collisions, imparting 
momenta mV;z to the particle as large as the 
quantity mV, given by Eq. (18), occurred with a 
frequency f as large as the quantity F from 
Eqs. (16)-(17), then Eq. (9) would give us an 
expression for S proportional to the square of r 
(through the factor f=F). For liquids of mo- 
lecular weight M+90 at room temperature 
(20°C) we should thus find 


S = (10°*!7?) (4u/3) ~ 2X 10°*7?, 


(19) 


This result must obviously be wrong for par- 
ticles that are not very small, because for such 
particles we know that Stokes’ law, Eq. (2), is 
valid, by which S is proportional to the first 
power of 7 only. That is, for large particles the 
friction is actually smaller than one would calcu- 
late on the assumption that the number of one-at- 
a-time collisions per second were as large as the 
number of thermal collisions F. 

The reason why f for not-too-small particles 
actually is much smaller than F, becomes clear, 
if we remember how extremely large is the num- 
ber of thermal collisions F against the particle in 
a liquid. Even against a particle as small as 
r=10-* cm we found 107 to 10'° collisions per 
second. (See below Eq. (17:3).) As each collision 
takes a finite time, this means that very often one 
collision has not yet been completed before other 
collisions start. Since the molecules of the liquid 
have a finite size, those new collisions cannot take 
place from exactly the same direction, so the 
average effect of simultaneity of collisions will be 
a partial neutralization of the effects of the single 
collisions. This neutralizing effect should not be 
confused with the friction discussed in the first 
section. The friction there considered was a con- 


8 Stokes’ law, Eq. (2), was not used before the second 
part of Eq. (15). 
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sequence of the velocity which the particle had 
already obtained, while here we discuss an effect 
that is not due to the velocity of the particle, but 
which, on the contrary, acts to prevent the mo- 
tion from starting at all. 

Thus we should not try here to use the sta- 
tistical considerations leading to Eqs. (16)—(18) 
for a calculation of the friction by kinetic con- 
siderations, as one could do in the case of a 
rarified gas,* for which f indeed would be equal to 
F for small particles. Such a derivation would 
simply establish statistically Eq. (19), which we 
obtained above from Eggs. (1) and (8). But Eq. 
(19) did not apply to a liquid. Indeed, the friction 
of a particle in a liquid is not caused by collisions 
of single molecules. As just explained, the effects 
of these single collisions neutralize each other 
almost completely. The fluctuations, in which this 
neutralization is incomplete, may be responsible 
for all resultant forces on the particle—thus 
causing the friction as well as the jerky Brownian 
motion itself—and it may be true? that the 
number of these fluctuations should be pro- 
portional to F* and thus to the first power of r. 
Yet, Stokes’ law does not follow from this by any 
simple kinetic considerations. In fact, as will be 
clear from the hydrodynamical derivation of 
Stokes’ law, the main effects causing the apparent 
friction are the following: (A) There will be a 
difference in internal pressure between front and 
back of the particle, due to a difference in N’ at 
these two spots. This is the main reason for the 
preference for head-on collisions mentioned in the 
argument leading to Eq. (1), and has little 
directly to do with the difference in average 
relative velocity of molecules and particle at 
these two spots. (B) Although the friction on the 
sides of the particle is greatly reduced by the fact 
that the liquid is taken along by the particle 
and thus is practically not moving with respect 
to its surface at all, yet there is a shearing force 
exerted by the liquid on the particle due to the 
internal friction of the liquid itself. Thus 7 enters 
Eq. (2). Therefore the shearing force depends on 
the velocity of the particle only indirectly through 
the rate of change of velocity of the liquid, with 
increasing distance from the particle. 

From the above it will be clear that the 
fluctuations in thermal collisions cannot simply 
be discussed statistically without taking into 
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TABLE III. Radii and frequencies computed from Stokes’ law.* 


105 -ryqy 105 -ry2 105 -rq3 105 -rq11+p0 1010 -y2 «9/7 10-5 -f’/r 
Liquid (cm) (cm) (cm) (cm~ g) (cm=! sec™ g) (cm=! sec“) 


Acetone 0.0009 0.0076 0.015 0.048 
Aniline ? 0.061 : 0.50 
Benzene 0.053 0.013 i 0.090 
Carbon tetrachloride 0.0064 0.0147 ; 0.187 
Cyclohexanol ? 1.0 : 7.65 
Ethyl acetate 0.0021 0.0094 ; 0.068 
Ethyl alcohol 0.0084 0.0285 i 0.18 
Ethyl ether 0.0047 0.0070 ; 0.040 
Glycerin ? 14 140 
Mercury ? 0.00178 .009t 0.193 
Methylene bromide ? 0.00097 0.0193 
Phenol ? 0.18t 0.33 1.54f 
Propyl acetate ? 0.011} 0.028 0.076t 
Water 0.013 0.015 0.017 H 


10-19 -F/r2 10710 -F3/r? 
(cm sec™!) 


X10!" 3 =3.3 10!" 
1017 
1017 
x 1017 
x 1017 
x 1017 
X< 1017 
X10" 
X10!" 
xX 108 
x 10!" 
X 10!7t 
X 10!" 
xX 108 


4.8X 10% 
4.0X 10° 
7.1X10% 
5.4X 105 
(5.8 X 10!) 


xX 
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ee 
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x 1017 

Xx 10!" 
X 1017} 
X 10847 
x 1017+ 
x10" 
x10" 
7X10 
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* The first four columns of figures give rt and rir in units 10-5 cm for particles with pp =1. For particles smaller than r1 or rj11, Stokes’ law is 
not valid. (See text.) The last four columns give trequencies in sec™ for particles \ with r =10-5 cm and po=1, calculated as if Stokes’ law were 
— valid. Here (1/r)*frequency of complete change of direction of particle; f’Sfrequency of effective collisions contributing to jerky motion; 

F= frequency of thermal collisions against particle. The subscripts 1, 2, 3 refer to calculation by means of Eq. (17:1), (17:2), or (17:3); for instance, 
F3 is obtained by substitution of Eq. (17:3) into Eq. (16). Where Stokes’ law gives wrongly f’>F for r =10-5 cm, we obviously have a case with 


r1 >1075 cm, that is, where = calculation of f’ or r we should not have applied Stokes’ law for ry =10-5 cm. (In the table, this is indicated by 
putting the values for f’ and for r between brackets.) 


+ Guessed on basis of the arbitrary assumption that perit~50 atmos. 
tt Guessed on basis of the arbitrary assumption perit~ 200 atmos. 
t Guessed on basis of the arbitrary assumption v.~ 1.4 105 cm/sec. 


account complicated correlations between the correct mean square velocity required by the law 
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velocities of molecules at different places in the 
liquid, if we want to explain the friction law, 
Egs. (1) and (2). As this is no simple problem, we 
shall here not enter into a further detailed dis- 
cussion of the mechanism of these effects. We 
shall, however, later discuss under what circum- 
stances the mechanisms (A) and (B), in the way 
in which they are taken into account in the usual 
hydrodynamical derivation of Stokes’ law, must 
fall short of giving an adequate description of the 
friction. 

Before we do so, let us first consider the case of 
sufficiently large particles, for which Stokes’ law, 
Eq. (2), with satisfactory accuracy gives the 
magnitude of the friction that is effective in 
damping the jerky Brownian motion. 

We could define the number of effective colli- 
sions per second f’ by 


f'm'& V2)ay=2SkT (20) 


with (V)ay given by Eq. (18); or, if Stokes’ law, 
Eq. (2), is valid, by 


f’m' V2) ay =12ekTor. (21) 


Comparing this with Eq. (9), we see that f’ is 
thus defined as the number of completely un- 
hindered collisions per second that is just suffi- 
cient to give the particle, under the influence of 
the friction law, Eq. (1) or Eqs. (1) and (2); the 


of equipartition of energy. We may add that /f’ is 
a rather conservative estimate of f, as actually 
the effective collisions (the ones that in the 
beginning of Sec. 1 were called ‘“‘collisions of 
exceptional vigor’) may be /ess-than-usually 
hindered rather than being completely unhin- 
dered. Then, of course, more effective collisions 
are needed than are given by f’; that is, f’<f. 

Substituting Eq. (18) into Eq. (21), we find 
that, for particles big enough for application of 
Stokes’ law, 

f' =9xrn/ 2p. (22) 
Thus, the larger the viscosity of the liquid, the 
greater number of thermal collisions should be 
effective in setting the particle into motion, in 
order to keep things going, if equipartition of 
energy is to hold for the particle. 

Let ¢ be the time after which the particle in its 
irregular motion has completely forgotten about 
its original velocity. This loss of memory is 
partially due to the friction in the liquid (which 
leads to t= 7). We can also say that the resultant 
of all impulses during the time ¢ should have a 
mean square comparable with the average square 
of the momentum of the particle.? Thus 


{ & mV ;} ad (ft) -mV >?) ay 
= (mv?) ay = (2m) (3kT/2) 
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or, by three times Eq. (9), 


t=m/2S = 37 =m/12ar. (23) 


Thus the motion of the particle can be con- 
sidered as rectilinear only during time intervals 
much shorter than the relaxation time 7. Com- 
paring the frequency (1/t), with which the 
particle is changing its direction, with the fre- 
quency of effective collisions conservatively esti- 
mated in Eq. (22), we see that the latter is larger 
than (1/t) by a factor 


f't=3m/8u, 


which is obviously very large. 

By Eq. (23), this also means that the relaxation 
time rt is many orders of magnitude longer than 
the time between two subsequent effective col- 
lisions,® contrary to the assumption made in 
Miiller-Pouillet’s Lehrbuch.! Even for a particle 
near the limit of visibility under a microscope 
using ordinary light (with r~10-® cm), in a 
liquid with molecular weight M=90, we find 
f'r=~2X10"%po, where po is the specific gravity 
of the particle. (Compare also Table III.) 

As a consequence of this fact, together with 
(V?) av + (vay =8u/3m<K1, (compare Eq. (18) 
with Eq. (8)), a particle in Brownian motion does 
not really describe a zigzag motion as is often 
pictured, but describes a path curved irregularly 
though not at one point excessively more than at 
another point. An apparent zigzag motion is ob- 
tained only if positions observed at definite time 
intervals are drawn in a diagram and then con- 
nected by straight lines.?° 

While f’ as given by Eqs. (21)-(22) is very 
large compared with (1/7), it must inevitably be 
smaller—often very much smaller—than F. This 
follows from the above qualitative arguments. 
For, indeed, the number f of collisions per second 
that cause the Brownian motion cannot be larger 
than the total number of collisions F given by 
Eq. (16), and due to the neutralizing effect 
probably will be much less than F. Our f’, defined 
by Eq. (20), will then obviously be still smaller, 
so that f’<F. Also, the square, m?V;,7, of the 
x-component of impulse picked up by the particle 
from a single relatively unhindered collision will, 


® Compare also M. von Smoluchowski, Ann. der Physik 
21, 756 (1906). 

10For physiological reasons the motion may perhaps 
make a similar appearance to the human eye. 


(24) 
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on the average, be not more than m*V,2),, as 
given by Eq. (18). Therefore, 


fm V i2) av < Fm V2) av. 


This relation, which we have already discussed 
qualitatively in the argument leading to Eq. (19), 
is now going to give us one condition for the 
validity of Stokes’ law for the friction damping 
the Brownian motion. For this purpose, we sub- 
stitute again Eqs. (16) and (18) into its right- 
hand member, and then use it for finding an 
upper limit for the right-hand member of Eq. (9). 
We thus find 


2SRkT <4nr®-3(c) ay: N’ + (8ukT/3), 
or, inevitably, 


S<4arN' ule) av/3. (25) 


Thence, Stokes’ law, Eq. (2), can be used in 
Eqs. (1)-(9) only if 6anr <42r?N’u(c)av/3, that is, 
a first condition for the validity of Stokes’ law for 
friction in Brownian motion is 


THM, 
where 


71 =9n/2N' UC) av. (26: 1) 


For particles smaller than this limit 73, Eq. (25) 
obviously gives S<4arr;N’y(c)av/3, or, by Eq. 
(26:1), S<6znr, so that for such small particles 
the friction in Brownian motion will be less than 
would follow from Stokes’ law. 

For the Brownian motion of smoke particles 
suspended in a gas, for which N’=N, and for 
which the viscosity 7 can be expressed in terms of 
the mean free path / of the molecules roughly by 
Enskog and Chapman’s equation" 7 ~0.5.Nyu(c) avl, 
our condition, Eq. (26:1), gives the well-known 
result” that Stokes’ law breaks down, if 7 is not 
much larger than /. (We actually find r>9//4.) If 
r is much smaller than /, Stokes’ law, Eq. (2), is 
indeed to be replaced by Epstein’s equation" 
S=4nr?Nu(c)av/3, which is in perfect agreement 
with our Eq. (25) with N’=N, and with < re- 

11For references, see, for instance, E. H. Kennard, 
Kinetic theory of gases (McGraw-Hill, New York, 1938), 
p. 147, Eq. (126b). 

12 See, for instance, E. H. Kennard, reference 11, p. 310. 

3 P. S. Epstein, Physical Rev. 23, 710 (1924), finds a 
slightly higher value for S, if he assumes that the colliding 
molecules are reflected diffusely by the particle. As we did 
not consider such possibility in our calculation of (Vz*)av 


in Eq. (18), the apparent disagreement between such find- 
ing and our inequality (25) is not surprising. 








placed by =, as is to be expected; for f= F and 
( Viz?) av = ( V2) av 

In the case of liquids, we shall take for N’ in 
Eq. (26:1) the value from one of the Eggs. 
(17:1, 2, 3). Thus, we find for 7; in cm 


3n(3perit —p) 








rn” 
2(C) avPPerit 
_— 1(3peric—p) M3 
= 0.60 x 10° x ———__——_;; (26:11) 
PPecrit 
9n 
rio= : (26:12) 
Apu, 
wkT 1600 X10 Xn 
Ti3>= a “) = ocean cmeeaemeat” (26:13) 
12Derit Mp’ 


Here, p’ is the critical pressure measured in 
atmospheres; 7 is the viscosity in cgs units 
(poises). With p=1 g/cm’, dUsouna~1.4X105 
cm/sec, p’~50 atmos and M#=90, this gives 
rie 1.6y X 10-5 cm or 113 ~3.4n X 10-5 cm. Here, 9 
will be very different for different liquids. For 
instance, at room temperature it is 0.001 for 
methylene bromide and about 15 for glycerin. 
For glycerin, therefore, Stokes’ law has already 
broken down for Brownian motion of particles of 
(more than) 10~* cm radius, while for methylene 
bromide Stokes’ law may be valid for the smallest 
particles visible through the microscope. In 
general, the applicability of Stokes’ law for the 
friction in Brownian motion will break down 
sooner (that is, for. larger size particles) in liquids 
of high viscosity than in liquids of low viscosity. 
Values for ri, 712, and 73, calculated from the 
data given in Table I, are listed in Table III. 

The reasons for our conclusion that Stokes’ law 
breaks down may be reformulated and summa- 
rized as follows: By Eq. (16), the number of 
collisions against a given particle is limited; by 
Eq. (18), the impulse per collision is limited. If 
Stokes’ law were always applicable to this 
Brownian motion, the friction would not decrease 
sufficiently rapidly with decreasing size of the 
particles (as compared to the decreasing amount 
of impulses available), so that the impulses would 
then, for very small particles, not suffice to give 
them the average speed required for validity of 
the equipartition law. Since, for Brownian’ mo- 
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tion, we believe more strongly in this equipartition 
law than in Stokes’ friction law, we must conclude 
that for such small particles Stokes’ law will not 
be valid. 

There are some more conditions limiting the 
validity of Stokes’ law for particles in Brownian 
motion. While the limitation to particles with 
r>ry given above was obtained here rather 
indirectly and, for liquids,“ without an attempt 
to understand the cause of the breakdown of 
Stokes’ law at r <1, the following two limitations 
will be derived directly by consideration of the 
hydrodynamical mechanism that explains Stokes’ 
law. 

First of all, we note that according to Stokes 
the pressure in the liquid at the front and at the 
back of the moving spherical particle would 


deviate from the average pressure in the liquid 
by an amount 


Ap = +3nv/2r. 


As the total internal pressure due to thermal 
motion will not be negative, the (—Ap) at the 
back of the particle cannot exceed the average 
internal pressure Pint ~27Perit in the liquid. We 
thus find that the mechanism leading to Stokes’ 
friction law cannot be effective, unless the condi- 
tion (3nv/2r) <27perit is satisfied, or, if we put 
v= (8kT/xm)} for the particle, we must have 





r>rny, 
with 
n 2kT\? 
| Hino (—) : (26:11) 
OPecrit \ mM 
As obviously!® 
riu/t13= (8/3) (u/m)*<1, * (27) 


we find that this new condition, r >7r1, is satisfied 
automatically as long as r>71, so that it does not 


really restrict the validity of Stokes’ law any 
further. 


The following condition is more important. It 


14 Note that for a liquid r; need not be as simply related 
to a mean free path / as for a gas, so that for a liquid there 
may be reasons for a limitation of the validity of Stokes’ 
law, different from those given by Epstein for a gas (see 
reference 13). 

16For further evaluation of r11, one should substi- 
tute m=(4r/3)ri1*p0 into Eq. (26:II), and then solve 
for ru. This gives r= (kT? /S4x? poPerit?) /§ = 0.237 
X10-5n2/5(p’)-2/59-1/5, or, for p’=50 and pol, rir 
=0. 059%/8X<10- 5 cm. 





476 FREDERIK J. 
requires that the hydrodynamical pattern of cur- 
rents in the liquid surrounding the particle, in 
case of validity of Stokes’ law, shall be able to 
establish itself within the time interval which 
elapses before the velocity of the particle changes 
its direction again. 

As seen above, this time ¢ is about half the re- 
laxation time r. Even in a time short compared 
with this ¢ (say ina time ~}3t~m/4S) the pattern 
should be able to propagate at least over a dis- 
tance =r. The speed of this propagation we 
identify with the velocity of sound in the liquid. 
This gives the condition }tv,>r, so that for 
Stokes’ law to be valid- for the friction in 
Brownian motion, we find ultimately the con- 
dition 

r>rnn, 
with 


ri = 18n/pov.. (26: 111) 
Now, since 


ri1/r12 = 8p/ po, 


(28) 


we find that (with p and po usually of the same 
order of magnitude) the limitation to the validity 


of Stokes’ law for friction in Brownian motion 
obtained here is slightly stricter than that given 
by Eq. (26:1). In Table III, values for rz in units 
10-5 cm are listed for particles with pp=1 in 
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various liquids. A comparison with the values 
obtained for 7; shows that in all cases 7 is indeed 
a little higher than 7, while 77/71 varies much 
less from one liquid to another than ry and 7; 
vary themselves. This is, of course, due to the 
common factor 7 in Eqs. (26:1) and (26:III). 

This approximate proportionality of ry and rm 
might seem to suggest that the mechanism of 
breakdown of Stokes’ law considered in Eq. 
(26: 111) could well be the very mechanism that, 
by Eq. (26:1), we tacitly assumed to exist in 
order to save the equipartition law for small 
particles. However, the fact that rm depends on 
po, While 7; does not, shows that this is not 
generally the case. In principle, we might find 
m>rmu for particles of extremely high specific 
gravity in extremely light liquids (with po>8p). 
For such highly exceptional cases, therefore, we 
have not yet found the mechanism responsible 
for the first breakdown of Stokes’ law for 
Brownian motion occurring at 7 ~7;.'* However, 
for the more usual case that p ~ pp we still can say 
that the validity of the equipartition law at least 
happens to be saved by the fact that hydro- 
dynamical current patterns causing Stokes’ fric- 
tion are prevented from being set up around small 
particles due to the lack of time for the pattern to 
establish itself within the relaxation time of such 
particles. 


Availability of Radioantimony-Beryllium Neutron Sources 


Oak Ridge National Laboratory has developed a rela- 
tively inexpensive radioantimony-beryllium neutron source 
which is now available for distribution to users in the 
United States. These sources will be allocated through the 
Isotopes Division, U. S. Atomic Energy Commission. 

The sources use a core of 32 grams of antimony metal 
surrounded by a ;y-in. thick beryllium metal cup, the 
assembly being enclosed in an aluminum jacket. The 
outside dimensions of this jacket are as follows: 


Cylinder containing source 1,%-in. diameter 
X14¢-in. long 

Aluminum plate protruding from 3 in. 
one end 

Hole in plate for handling 


Over-all length of assembly 


1 in. 
3 in. 
The entire unit is placed in the pile and the antimony 


activated to the desired level. Antimony (half-life 2.8 
days) and Antimony™ (half-life 60 days, gamma-ray 


. 


energy 1.72 Mev) are both formed in this irradiation, but 
the source can stand for two weeks before shipment to 
allow the Sb™ to decay. Saturation activity of the Sb'™ 
is approximately 3.5 curies; however, normally an activa- 
tion time of 120 days will be used which will give an activity 
of about 2.6 curies. A source activated for 69 days and 
cooled for 16 days gave an unshielded gamma-radiation 
reading of about 20 roentgens/hour at one foot. Other 
readings followed the inverse square very closely. Shorter 
activation periods may be arranged if a lower activity is 
satisfactory. After the activity of the Sb'* has decayed, 
it may be reactivated by reinsertion in the pile. 

It has been found that the saturated source will emit a 
total of about 8X10 neutrons per second or 9.810! 
neutrons/cm?/second at one inch, if it is assumed that the 
source behaves as a point source. The neutrons are approxi- 
mately monoenergetic. Arrangements must be made to 
provide a suitable storage space for the source; a lead pig 
with walls four inches thick would probably be sufficient. 





On Newton’s Third Law and the Conservation of Momentum 


E. Geryuoy 
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T is customary in both elementary and ad- 
vanced treatments of mechanics to derive the 
conservation laws of linear and angular momen- 
tum from Newton’s third law of motion.! On the 
other hand, the forces between moving charges 
do not in general obey Newton’s third law. One 
of the simpler instances of this fact, the problem 
of two charged particles moving in mutually 
perpendicular directions, has been discussed in a 
number of recent papers.?-* The usual answer to 
the difficulty seemingly raised by this failure of 
Newton’s third law has been given by Keller.’ 
The conservation laws remain valid when the 
momentum of the electromagnetic field is taken 
into account, from which it is concluded that the 
difficulty is only apparent and not real. 

It appears to this writer that the above answer 
is not wholly to the point. Mechanics is con- 
cerned with the motions of material bodies. If 
these bodies are assumed to consist of Newtonian 
particles, then Newton’s third law explains and 
elucidates the experimentally observed fact that 
material bodies conserve mechanical momentum. 
If the particles are charged, introduction of an 
electromagnetic field possessing momentum pre- 
serves the conservation laws. However, it is 
scarcely cogent to infer from this demonstration 
of what amounts to new conse:vation laws that 
material bodies, which are known to be a vast 
conglomeration of moving charges, will conserve 
mechanical momentum in the absence of external 
forces. As a matter of fact material bodies do not 
necessarily conserve momentum in the presence 


1 We follow Whittaker and explicitly include in Newton's 
third law the assumption that the forces between two 
particles are directed along the line joining them. This as- 
sumption of central forces was not included in Newton’s 
original formulation of the third law, but as is well known, 
is required for the usual proofs of the conservation law of 
angular momentum. Since throughout this paper we treat 
the conservation laws of linear and angular momentum on 
an equal footing, we are merely avoiding constant repeti- 
tion of the awkward phrase “‘and along the line joining 
them.” E. T. Whittaker, A treatise on the analytical 
dynamics of particles and rigid bodies (Dover Publications, 
1944), p. 27. 

2J. M. Keller, Am. J. Physics 10, 302 (1942). 

( “ne Robertson, Phil. Mag. 36, 32 (1945), and 39, 162 
1 fe 
*W. A. Tripp, Phil. Mag. 38, 149 (1947). 


of an electromagnetic field, as is evidenced, for 
example, by the observable effects of radiation 
pressure. Conservation of momentum is a fact of 
experience only because electromagnetic forces 
are very small, so small as to be usually negligible, 
and because in addition a body in thermal equi- 
librium with the radiation surrounding it radiates 
as much momentum as it receives, on the aver- 
age. If electromagnetic forces were very much 
larger, i.e., if the charge on the electron were very 
much larger, the linear momentum or angular 
momentum of isolated material bodies might 
very well show large fluctuations, even when 
averaged over relatively long time intervals. 

This brief discussion exemplifies the obvious 
point that using Newton’s third law to prove the 
conservation laws of linear and angular mo- 
mentum involves an assumption: that it makes 
good sense to regard material bodies as composed 
of Newtonian particles or, in other words, of mass 
points obeying Newton’s third law. Unfortu- 
nately, this assumption is rather difficult to 
defend. For this reason it is of interest to note 
that the conservation laws of linear and angular 
momentum can be justified by an argument 
which does not depend on special assumptions 
regarding the nature of the forces between the 
elementary particles composing material bodies. 
In particular, it is not necessary specifically to 
invoke Newton’s third law. Naturally this argu- 
ment also involves a physical assumption. It is 
sufficient to assume that the internal forces in 
material bodies do no work in any virtual dis- 
placement which is merely an infinitesimal trans- 
lation or rotation. This assumption is in many 
respects more in keeping with our present under- 
standing of the nature of matter than is the 
assumption that material bodies consist of 
Newtonian particles. The elucidation of the 
statements of this paragraph is presented in the 
following section. 

Ultimately the proof of the conservation laws 
of linear and angular momentum must be based 
on quantum mechanics. In order that a classical 
proof be meaningful it is necessary that the 
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fundamental assumptions on which classical me- 
chanics rests, especially the assumption of par- 
ticles following definite trajectories, make good 
sense from the standpoint of quantum mechanics. 
This paper attempts to demonstrate that while 
perhaps no defense of the assumptions of classical 
mechanics can be completely satisfactory, a 
better defense can be made if it is not specifically 
required that every one of these particles obey 
Newton’s third law. Although Newton’s third law 
may no longer be obeyed in detail, nonetheless 
systems of particles continue to behave as if 
Newton’s third law were obeyed, in the sense that 
the conservation laws of linear and angular mo- 
mentum remain valid. Moreover, justifying the 
conservation laws without detailed reference to 
Newton’s third law probably gives more insight 
into their physical significance than merely 
inferring them from the subsumed truth of 
Newton’s third law. 

It should be remarked that the result, that the 
invariance of the Lagrangian with respect to 
translations and rotations implies the conserva- 
tion laws of linear and angular momentum, has 
been proved long ago for systems describable by 
a Lagrangian.® (The corresponding invariance of 
the Hamiltonian guarantees the conservation 
laws of linear and angular momentum in quantum 
mechanics.) Basing the proofs of these conserva- 
tion laws on the existence of a Lagrangian seems, 
however, less significant physically than the 
proofs given below. It is by no means obvious 
that any system of particles has a Lagrangian, or, 
if there is a Lagrangian, that it depends only on 
the relative positions and velocities of the par- 
ticles, and does not depend on any other acces- 
sory quantities. To illustrate, the motion of any 
system of charged particles can be described by a 
Lagrangian, but the Lagrangian is an explicit 
function of the electromagnetic field variables, as 
well as of the particle positions and velocities. 
The Lagrangian is not invariant with respect to 
translations or rotations of the particles alone, 
and consequently mechanical momentum, the 
momentum of the particles only, is not conserved. 


Justification of the Conservation Laws 


The usual derivation of the laws governing the 
rate of change of linear and angular momentum 


Whittaker, see reference 1, p. 58. 
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of a system of particles proceeds as follows.* The 
equation of motion of the ith particle of a system 
is 

Fi+> fis=mii, 


jAi 


(1) 


where F; is the external force acting on the par- 
ticle and f,; is the internal force exerted on particle 
1 by particle 7. Since f;;= —f;;, summing Eq. (1) 
over all particles leads to 


a F;= ae m Xj. (2) 


Furthermore, since f;; is directed along the line 
between particles 7 and j, i.e., since fj; is parallel 
to r;—r; it also follows that 


p (r;XF;) = x m(r;X*¥;). (3) 


Equations (2) and (3) express the usual rules 
that the sum of the external forces equals the rate 
of change of momentum of the system, and the 
sum of the applied torques equals the rate of 
change of angular momentum. 

Now Newton’s third law can either be regarded 
as an independent postulate regarding the allowed 
types of forces between particles, or it can be 
deduced, following Mach,’ from the definition of 
mass. In either event, its use to deduce Eqs. (2) 
and (3) involves the assumptions (a) that any 
material body, such as a rigid body, is composed 
of a definite and immutable number of point 
particles, and (b) that the accelerations of these 
particles, when no external forces are acting, can 
always be understood in terms of pairs of forces 
between pairs of particles, i.e., in terms of two- 
body forces which obey Newton’s third law. 
When an attempt is made to identify these par- 
ticles, however, it is apparent that the above 
assumptions are inconsistent with present under- 
standing of the structure of matter. Thus, if the 
particles are supposed to be atoms or molecules, 
it can scarcely be maintained that the internal 
forces between such particles occur always in’ 
pairs which obey Newton’s third law. Rather, 
because atoms and molecules have structures 


6 Cf. e.g., G. Jom Theoretical physics (G. E. Stechert and 
Co., 1934), pp. 102 ff. 


Mach, The science of mechanics (Open Court 
Lindsay 
a physics (John Wiley and 

A. Simon, Phil. Mag. 38, 


Publishing Co., 1942), PP. 266 ff. Cf. also, R. B. 
ndation. 


and H. Margenau, Fou 
Sons, Inc., 1936), pp. 91 ff. and 
888 (194 7). 
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which can and do change, as when a molecule is 
polarized by the electric field of a neighboring 
atom, the forces between such particles depend 
not only on their relative distances but also on 
their relative orientations and on the positions 
and orientations of neighboring particles as well. 
Such forces are really many-body forces; the 
force acting on any one of three contiguous atoms 
depends on the simultaneous positions, orienta- 
tions, and polarizations of the other two, and can 
only with extreme artificiality be regarded as re- 
sulting from pairs of equal and opposite forces 
between the various atom pairs. In any event, 
the forces between such particles are not neces- 
sarily directed along the line connecting them as 
is evidenced by the simple examples of the elec- 
trostatic forces exerted by a dipole on a single 
charge or on another dipole. 

These difficulties are not removed by identi- 
fying the particles with electrons and the con- 
stituents of atomic nuclei. Electrons obey the 
laws of quantum mechanics, not classical me- 
chanics, and, as evidenced by interference and 
diffraction phenomena, are not at all well lo- 
calized in space in the sense demanded of a 
Newtonian particle. This makes relatively unim- 
portant the fact that, because of their spins, the 
forces between two electrons, like the forces be- 
tween two electric dipoles, are not necessarily 
directed along the line connecting them. Similar 
objections can be cited to identifying the particles 
of Eq. (1) with the protons and neutrons (and 
presumably also mesons) which compose atomic 
nuclei. 

The foregoing arguments have shown, not that 
it is impossible to carry through proofs of Eqs. (2) 
and (3) using Newton’s third law, but that such 
proofs are based on hypotheses which can only 
with difficulty and by unnatural artifices be 
reconciled with present knowledge of the struc- 
ture of matter. Equations (2) and (3) can, how- 
ever, be established in the following way. 


Let the equation of motion of the ith particle 
be 


F,+f;=m¥#;i, (4) 


where f; is the total internal force exerted on 
particle 7, and may be considered formally defined 
by Eq. (4) in the form 


f;=m?;—F ;=m,(¥;—8)). » (5) 
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In Eq. (5) 8; is the acceleration the ith particle 
would have if it were free to move under the 
influence of F; alone. The internal force f; is then 
just proportional to the difference between the 
acceleration ¥; which the particle 7 actually re- 
ceives and the acceleration it would have had if 
internal constraints, expressed through f;, were 
not effective. Equation (4) differs from Eq. (1) 
in that the internal force f; is not necessarily a 
sum of two-body forces f;;, obeying Newton’s 
third law. Because of this, many of the objections 
raised above to identifying the particles with 
atoms or molecules are no longer pertinent. Like 
electrons, atoms and molecules obey the laws of 
quantum mechanics and have wave properties, 
but there is good reason to believe that at ordi- 
nary temperatures and pressures the classical 
approximation of definite, well-localized tra- 
jectories is usually valid for atoms and molecules. 
Even where this approximation is not altogether 
valid, it is certainly far better for atoms or mole- 
cules than for electrons, at ordinary temperatures 
and pressures. Comparison of the theories of 
imperfect gases and of electronic conduction illus- 
trates the difference, in terms of approach to 
physical reality when identified as classical par- 
ticles, between atoms or molecules and electrons. 
It is this fact, that it seemingly makes good 
physical sense to call atoms and molecules par- 
ticles when the requirement of two-body forces 
obeying Newton’s third law is abandoned, that 
justifies the argumeni below based on Eqs. (4) 
and (5), whereas the use of Eq. (1) to obtain 
Eqs. (2) and (3) is not justifiable. These quali- 
tative considerations are formulated more pre- 
cisely at the end of this section. 

Granting the validity of Eqs. (4) and (5), the 
proof of the conservation laws proceeds as follows. 
Let dr; be any set of virtual variations of the 
positions of the particles compatible with the 
constraints, such that 


i f;-6r;=0. (6) 


Then it follows from Eq. (4) that for this set of 
virtual variations 


(7) 


a (F; —m#;) oT) 0. 


Equation (7) looks like D’Alembert’s principle, 
but differs from D’Alembert’s principle in that 
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Eq. (7) holds only for those variations for which 
Eq. (6) is satisfied. With arbitrary types of in- 
ternal forces Eq. (6), and therefore also Eq. (7), 
is not, in general, satisfied for all sets of virtual 
variations compatible with the constraints. For 
example, in the problem of two particles m, and 
ms moving in free space under the influence of 
their mutual gravitational attractions, the gravi- 
tational forces are the internal forces f; of Eq. (4), 
while the external forces F; are zero since the 
system is isolated. For these internal forces Eq. 
(6) is not satisfied for arbitrary sets of virtual 
variations of the positions of the two particles, 
say a set in which 6rz=0 but 6r: is along the 
direction of f. 

Consider now the two sets of virtual variations: 


(8a) 
(8b) 


ér;= du=a constant. 
ér;= (6¢ X13). 

Substitution of Eq. (8a) in Eq. (7) gives 
> (Fi—m#;)-du=0. (9) 


t 


Substituting Eq. (8b) in Eq. (7) gives 


= (F:—m#;) -(5¢ X1,) 
=> [r:X(F:—m#,) ]-5¢ 
=D (:XFi—maiX¥,)-66=0. (10) 


Equations (9) and (10) have been derived on the 
assumption that the variations (8a) and (8b) 
satisfy Eq. (6). Accepting this assumption for the 
moment, Eq. (9) implies that if the constraints 
permit an infinitesimal translation éu of the 
entire system in any given direction, then the 
sum of the components of the external forces 
along that direction equals the rate of change 
along that direction of the total momentum of the 
system. Equation (10) means that if the system 
can be rotated through an infinitesimal angle 5¢ 
about an axis through the origin, then the sum of 
the external torques about that axis equals the 
rate of change of total angular momentum along 
that axis. Since an isolated system can, in general, 
be translated or rotated in an arbitrary fashion, 
and since any constraints preventing such trans- 
lations or rotations in a nonisolated system can 
always be replaced by the imposition of suitable 
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external forces, it follows that Eqs. (9) and (10) 
are equivalent to Eqs. (2) and (3). 

It remains then to consider under what circum- 
stances the variations (8a) and (8b) satisfy Eq. 
(6). The variations (8a) and (8b) have the 
property that they leave unchanged the distances, 
angles and relative velocities between all parts of 
the system. In fact, even the relative accelera- 
tions and higher time derivatives are unchanged. 
The infinitesimal rotation and translation amount 
merely to setting the center of mass of the system 
at a different point in space, with a somewhat 
different orientation relative to fixed axes in 
space. The relative motions of the component 
particles of the system are not changed. There is 
no more effect on the motion of the system than if 
the origin had been translated and the coordinate 
axes rotated. Thus, in order that Eq. (6) be 
satisfied, it is sufficient that the internal forces do 
no work in any virtual displacement of the 
system that leaves all relative distances, angles, 
velocities, accelerations, etc., invariant. In fact, 
the equations >>; f;-du=)>);f;-6¢xXr;=0 are a 
necessary and sufficient condition for the con- 
servation laws, i.e., these equations are implied 
by and imply Eqs. (2) and (3). 

It is apparently a fact of experience that this 
condition is satisfied for ordinary material bodies, 
or, in other words, the internal forces of material 
bodies do no work in virtual displacements which 
merely translate and reorientate the body in 
space. Certainly there is no evidence either from 
mechanics or thermodynamics that a translation 
or rotation of an isolated macroscopic body can 
change its internal energy; presumably, if the 
internal energy is unchanged in such a translation 
or rotation the internal forces have done no work. 
If the internal energy of the body depends only 
on the relative distances, angles, velocities, etc., 
of its component parts, there is obviously no 
change in internal energy during translations or 
rotations. The mere requirement that the concept 
internal energy of a body have meaning is almost 
sufficient to guarantee that the internal energy 
depends only on relative distances, angles, ve- 
locities, etc. Nonetheless, the statement that the 
internal forces do no work, or the internal energy 
does not change during translations or rotations 
remains an assertion justifiable only on grounds 
of experience. Actually the particles of a body are 
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not at rest but are vibrating about positions of 
equilibrium, so that the internal forces may be 
doing work at any given instant, and there is a 
continual exchange of energy between the body 
and the vacuum. However, on the average, the 
internal energy does remain constant or nearly 
constant as long as the average relative distances, 
velocities, etc., of the particles are unchanged. 

The demonstration of the conservation laws 
therefore involves basically the assumption that 
space is isotropic; if space is isotropic, the in- 
ternal energy of an isolated system cannot depend 
on position or orientation relative to an external 
set of fixed axes, and the conservation laws then 
follow. The proofs of the conservation laws which 
rest on the invariance of the Lagrangian or 
Hamiltonian with respect to translations and 
rotations® also involve basically the assumption 
that space is isotropic. There is no need to justify 
the assumption of isotropy for truly empty space, 
which is isotropic almost by definition. Me- 
chanical systems are embedded, however, not in 
empty space but in a space containing electro- 
magnetic fields, neutrino fields, etc. Only if the 
‘‘system’’ includes all such fields can it be said to 
be embedded in an empty space whose isotropy 
implies generalized momentum conservation laws. 
If the fields are omitted from the system, i.e., if 
the system is thought to be defined completely by 
specifying the positions and velocities of its com- 
ponent mass-points, the mechanical momentum 
of an isolated system is conserved only under 
specific though widely occurring circumstances. 
The necessary and sufficient condition for such 
conservation, or alternatively, the necessary and 
sufficient condition that space appear isotropic to 
the mechanical system, is that the forces acting 
on the collection of mass-points which form the 
mechanical system do no net work in translations 
or rotations.*® 


8 Neither this statement nor any other part of this paper 
is really new. The essence of the derivation of the conserva- 
tion laws which has been given here is contained in 
Whittaker (reference 1). However, the details of the argu- 
ment, whereby the conservation laws are deduced from an 
assumption of isotropy without direct reference to a 
Lagrangian, are not given explicitly in any source known to 
the writer, (e.g., in such standard works as Whittaker, 
Appell, the Handbuch der physik, or in G. D. Birkhoff, 
Dynamical systems, American Mathematical Society Collo- 
quium Publications, v. 9, 1927). More importantly, the 
physical significance of the work done by the forces in trans- 
lations and rotations, and the relationship of this work to 
the assumption of isotropy, is nowhere sufficiently stressed. 
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The foregoing paragraphs lead to the following 
summarizing conclusion. It can be proved that 
the conservation laws of linear and angular mo- 
mentum are valid provided the internal forces of 
a system do no work during translations or rota- 
tions of the system. For an isolated body the 
assertion that the internal energy of the body 
does not change for translations or rotations 
implies that the internal forces do no work in 
translations or rotations. It is not inconceivable 
that the internal energy may change and the 
internal forces may therefore do work, but it 
appears from experience that the energy changes 
are negligible for ordinary macroscopic bodies 
(recalling that the internal forces always are 
defined by Eq. (5)). It follows then that such a 
body will behave as if the forces between its 
component parts obeyed Newton’s third law.° 
Similarly, in a collision between two macroscopic 
bodies, it can be assumed that the internal energy 
does not change during translations or rotations 
of the system as a whole, so that the bodies will 
interact as if the forces obeyed Newton's third 
law, i.e., linear and angular momentum will be 
conserved during the collision. 

To complete the argument, some remarks con- 
cerning the validity of Eq. (4) should perhaps be 
made. It is stressed that the problem of the 
applicability of Eqs. (4) or (1) is wholly divorced 
from the problem of the validity of Newton's 
third law. Whether atoms, molecules, or electrons 
are regarded as the fundamental classical par- 
ticles, Newton’s third law remains invalid. For 
Eq. (4) to be applicable toa system of quantum- 
mechanical particles, (a) the particles must be 
represented by wave packets which sptead suffi- 
ciently slowly that a trajectory can be inferred, 
and (b) the wave packets must not overlap, 
especially when the system is composed of 
identical particles. These requirements are gener- 
ally necessary. It can be shown that such packets 
can be constructed for the molecules of a gas 


®It should be noted that it is not possible to justify 
conservation of linear and angular momentum of a rigid 
body simply on the grounds that the isolated body at rest 
remains at rest, implying that the internal forces must form 
a system whose sum and resultant moment are zero. When- 
ever the body translates nonuniformly or rotates, strains 
are set up which cause the internal forces to differ from 
those obtaining when the body is at rest. 
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provided that 


Nh* 
———-— <1. 
(2xmkT)} 


(11) 


The inequality (11) is the well-known condition 
for the validity of classical statistics ;!° it is known 
to be satisfied for all gases at temperatures above 
1°K, but is not satisfied for assemblies of free 
electrons at temperatures below 2000°K. 

On the other hand, such packets cannot always 
be constructed for atoms or molecules. In dia- 


mond, for example, the failure of the law of . 


Dulong and Petit shows that the carbon atom is 
almost never in states of high quantum number. 
Moreover, even in gases, it is still necessary to 
demonstrate that the internal coordinates of the 
molecule can be neglected, so that the molecule 
can be regarded.as a quantum-mechanical par- 
ticle and characterized by a single coordinate, 
that of its center of mass. This means that the 
motion of the center of mass must be separable 
from the internal motions, which is just another 


10 R. H. Fowler, Statistical mechanics (Cambridge, 1936), 
pm. 35. 
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way of saying that the Hamiltonian describing 
the isolated molecule must be invariant to trans- 
lations. (It can, of course, be invariant on the 
average only, because of interaction with the 
electromagnetic field of the vacuum.) 

In fine, it is psychologically unsatisfying, from 
the standpoint of classical mechanics, simply to 
assume ab initio that the acceleration of the 
center of mass of an isolated macroscopic system 
is proportional to the external forces only. The 
classical-mechanical proof of the conservation 
laws divides the system into subaggregates, each 
of which presumably conserves momentum when 
isolated. The quantum-nfechanical justification 
of this assumption, if it can be given at all, in- 
volves an appeal to the invariance of the 
Hamiltonian or total energy with respect to 
translations. A classical proof of the conservation 
laws based on a similar invariance of the energy 
of the entire macroscopic system is more de- 
fensible and probably more significant than a 
proof based wholly on the assumption that the 
subaggregates obey Newton’s third law. 

I should like to thank Dr. Sidney Benson for 
his comments. 


Network Representation of Input and Output Admittances of Amplifiers 
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NPUT and output admittances are among the 

most important characteristics of amplifiers. 
They are defined as the ratios of the input 
(output) current to the input (output) voltage 
and are computed on the basis of the equivalent 
circuit. Nichols! and Miller? derived the ex- 
pression for the input admittance of a linear 
triode in 1919; Chaffee* in 1929 gave the first 
general derivation and a systematic analysis of 
both input and output admittances. 


1H. W. Nichols, “The audion as a circuit element,” 
Physical Rev. 13, 404-414 (1919). 

2J. M. Miller, ‘‘Dependence of input impedance of a 
three-electrode vacuum-tube upon the load in the plate 
circuit,” Bur. Standards Bulletin, N. 351 (1919). 

7E. L. Chaffee, “Equivalent circuits of an electron 
triode and the equivalent input and output admittances,”’ 
Proc. IRE 17, 1633-1648 (1929). 


This analysis, because of the complexity of the 
expressions, is usually obtained by considering 
separately the real and the imaginary compo- 
nents and introducing such simplifying or limit- 
ing assumptions as the case may suggest. It is 
the purpose of this paper to show that the input 
and output admittances may be represented as 
the admittances of linear passive networks, ob- 
tained from the equivalent circuit with certain 
additions of parallel or series branches. The 
consideration of these networks offers direct and 
physical interpretation of the impedance, show- 
ing the importance of the components and allow- 
ing the use of graphic phasor analysis. Through- 
out the following the assumptions of linearity 
and of negative grid bias will be made. 















Theory . 


The correspondence between the input vari- 
ables (Vi, J:) and the output variables (V2, 2) 
in a linear quadripole may be expressed by means 
of a bilinear transformation of one of the follow- 
ing types: 


Vy=Zulht+Ziels 
| : 


Th=Zal;+Zoel2 

Th= YurVit YieV2 
| I3= YoiVit YooV2 

Vi =AiuV2tArele 
| Iy;=AaV2+Aool>. 


(2) 


(3) 


If the polarities are defined as in Fig. 1, the 
coefficients of Eqs. (1), (2) and (3) can be com- 
puted or measured from the quadripole on the 
basis of relations of such type as 


Zu=(Vi/h)r=0, Yu=(hh/Vi)ve=0, 
Ai: =(Vi/ V2)12=0, 


For instance Z1(¥Yu) is the input impedance 
(admittance) at open (short) output circuit, 
—Z2(— Ye2) is the output impedance (admit- 
tance) at open (short) input circuit. 

The systems (1), (2), (3) are unequivocally 
determined by the matrices of their coefficients, 
which we shall indicate by ||Z||, || Y|!, ||Al] re- 
spectively. Among the elements of these matrices 
the following relations exist : 
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where |Z|, | ¥|, |A| are the determinants of 
Z|, || Vl], ||Al] respectively. The circuit analysis 
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Fic. 1. Definition of polarities in a quadripole. 
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of a linear quadripole may, therefore, be handled 
by means of the matrix theory, which presents 
some organizational advantages. Brown and Ben- 
nett‘ have shown recently the application of the 
theory to the solution of vacuum tube problems. 
A linear amplifier may be represented by its 
equivalent circuit which constitutes an active 
quadripole. To find the matrix of the latter, one 
can apply Kirchhoff’s equations or consider the 
quadripole as a combination of m simpler net- 
works. In the latter case, assuming for example 
m =2, one has ||A|| = ||A,||+||Aol| if the combina- 
tion is in cascade, ||Z|]=||Z:1||+||Z2|| if it is in 
series, || Y|| =|] ¥1||+|| Yel] if it is in parallel. For 
instance the equivalent circuit of a grounded- 
cathode amplifier may be considered as the 
parallel combination of the basic passive quadri- 
pole and of the active quadripole which repre- 
sents the internal constant current generator 
shown in Fig. 2. Its admittance matrix is written 


therefore 
0 0 
+ 
—pY3 0 


Vit. —Y2 
Y2 —Y.—Y3 
Yi+Y, —YF¥2 
coe —¥2— vt 
The input admittance of a quadripole closed 


on a load of admittance Yz, (impedance Zz) may 
be obtained from Eq. (3) as 
Vin . i _4 21 VitA aol _A 21 +A 22 Y, 
Vi AuVitAnl: AutAnwYz 
Using Eq. (6) there results 
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Fic. 2. Equivalent circuit of grounded-cathode am- 
plifier as parallel combination of an active and a passive 
quadripole. 


‘J. S. Brown and F. D. Bennett, “The application of 
ata) to vacuum tube circuits,’’ Proc. IRE 36, 844-852 
1948). 
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Fic. 3. Network representation of Eqs. (11) and (12). 


Similarly the output admittance of a quadripole 
with a source admittance Ys (impedance Zs) is 


—|Y|—YnYs ZstZu 1 
Yue FX (8) 
Yut Ys —ZZ5— |Z| Dis 


Division of numerator by denominator in Eqs. (7) 
and (8) gives expressions of type 


m+ pn 
Y=/+ 
p+uq 
ro , 
iis p +eq 
m’+ yn’ 


(9) 


(10) 


where /+m/p is the admittance (s+ p’/m’ being 
the impedance) of the corresponding passive 
quadripole obtained when p=0. 

The problem of finding a 2-terminal network 
with the given admittance (9) or impedance (10) 
may be solved in several ways. For instance we 
may write Eq. (9) in the form 


1 
+t, 
p/m+u(q/m) q/n+p/un 
or Eq. (10) in the form 
1 
PY aaadadaan atk 
m'/p'+-u(n'/p')  n'/q' +(m'/uq’) 


Y=l 


(12 


which are represented by the networks of 
Fig. 3(A, B) ; the circuit elements are represented 
as admittances in Fig. 3(A) and as impedances in 
Fig. 3(B). In this representation g together with 
the parallel branch 7 in (a) or with the series 
branch s in (b) constitutes the basic passive 
quadripole, to which the representation is re- 
duced when »=0. Simpler networks may be 
obtained in particular cases, when g or p=0 in 
Eq. (10), 2 or m=0 in Eq. (11), as will be shown 
in the following examples. 

In Table I are condensed the results obtained 
for the fundamental cases of the grounded- 


cathode, grounded-plate, grounded-grid ampli- 
fiers. The matrices of their respective equivalent 
circuits are: 

For the grounded-cathode amplifier: 


an | 
Y2—u¥; —¥2—Yill 
For the grounded-plate amplifier: 


YitY2 —Y¥2 
(Yil= 
YotuY3 — Yo—YV3—pnY3 


For the grounded-grid amplifier: 


Z:+Z, —Z2 
\Z\| , 
ZotuZ, —Z2—-Z3—puZe 


In the latter case impedances have been used 
instead of admittances for convenience of writing 
formulas. 

With reference to the notations of Table I and 
to the forma! representations (9) and (10), the 
following results are obtained : 

For the grounded-cathode amplifier : 
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For the grounded-plate amplifier : 
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For the grounded-grid amplifier : 
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Fic. 4. Network representation of the input admittance 
of a linear grounded-cathode amplifier (A) with resistive 
load, (B) with inductive load. 


Fic. 5. Network representation of the reactance tube 
and its output admittance. 


The reader may verify directly the derivation of 
the corresponding network representations as in- 
dicated in Table I. 

When the amplifier circuit is complicated with 
feedback branches the simplifications in Eqs. (9) 
and (10) may not take place and the network 
representation is of the general type of Fig. 3. 


Examples 


The above method of representation may be 
applied in many cases with advantage. For in- 
stance, in the classical case of the grounded-cath- 
ode amplifier with resistive load, the network rep- 
resentation assumes the form of Fig. 4(A), where 
Y,;=(1/R,)+jwCo, Yo=jwCop, Ya=(1/Ry), Vu 
=(1/Rz) and 


. -* 
— =— + 


1 R> 
(42 ; 
rr +“ 


j wh Cop Ry 


This shows how the input capacitance approaches 
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a maximum value of C,.+(u+1)C,, when the 
load impedance is very high.® 

When the load is inductive, the above network 
assumes the form of Fig. 4(B), where 


1 R, 1 1 

aa (1 “ )+ , 

. Be w*CopL jouCyp 
A negative resistance is then obtained, provided 
wCepl <1. 

Another interesting example is offered by the 
reactance tube, Fig. 5. Neglecting the capaci- 
tance C,,, we have ¥1=jwC, Yo=1/R, Y3=1/R,, 
Ys=0; consequently 


1 R, _ CRR, 
— == + jo——. 
. ae Mu 


The output admittance of the reactance tube is 
then represented by the network of Fig. 5(C), 
which gives complete information on the actual 
operation of the tube, often approximated by 
reducing the whole admittance to jwCRR,/u. 
Other examples could be worked out easily; 
for instance one can verify that the input ad- 
mittance of a cathode follower (grounded-plate 
amplifier) presents a negative conductance when 
the load is capacitive. In our equivalent net- 
work the admittance Y’ of the added branch is 


Yo Y, 1 jwC gx 
¥ -—(14+2 -( - Ja+R¥. 
Cd 3 MING M 
1 


; Cox Ry s 
=—+ jo—+ Yi] —+ jo 
BR, B 


a) 


BR, BL 


So the real component of Y’—(1/uR,) is positive 
or negative accordingly if the load is inductive 
(Yr=1/jwL) or capacitive (Yz=jwC). 


5F. E. Terman, Radio engineering (McGraw-Hill, New 
York, 1947), p. 365. 

®K. Schlesinger, ‘Cathode follower,’’ Proc. IRE 33, 
843-855 (1945). H. Reich, “Input admittance of cathode 
follower amplifiers,"’ Proc. TRE 35, 573-576 (1947). 


Anybody can make history; only a great man can write it.—OscaR WILDE. 





Albert A. Michelson at Case 


R. S. SHANKLAND 
Case Institute of Technology, Cleveland, Ohio 


LBERT ABRAHAM MICHELSON (Fig. 1) 

was the first professor of physics at Case 
Institute of Technology (then Case School of 
Applied Science), holding that position from 1881 
until 1889. During these years he performed some 
of his most important researches in optics. 
Michelson was born in Strelno, Germany, now 
Strezelno, Poland, on December 19, 1852. The 
family moved to Virginia City, Nevada when he 
was a small boy and later to San Francisco where 
Albert attended high school. It was his father’s 
wish that Albert become a naval officer, so the 
young man took the examination for the U. S. 
Naval Academy which he passed satisfactorily, 
only to find that all the Congressional appoint- 
ments had already been filled. Nothing daunted, 
young Michelson traveled from California to 
Annapolis to apply personally to the Superin- 
tendent of the Academy for admission. He must 
have made a remarkable impression because an 
interview was arranged permitting Michelson to 
present his case in person to President Grant. 
Here too he met disappointment as Grant had 
already made the ten presidential appointments 
at his disposal. However, Michelson so impressed 
the President with his serious determination and 
abilities that Grant made an exception and ap- 
pointed him as an extra man for the Naval 
Academy that year. 

Michelson was graduated from the Naval 
Academy in the class of 1873 and after spending 
the usual two years at sea asa passed midshipman 
returned to the Academy as instructor in physics 
and chemistry from 1875 until 1879. It was 
during these years that he was toid to “explain 
the velocity of light to the midshipmen.” To 
accomplish this “laboratory demonstration,” 
Michelson set up in November 1877,' simple 
equipment following Foucault’s rotating mirror 
method? of 1862 and obtained results accurate to 
an unexpected degree. This led him to refine the 
method using a longer light path along the north 
sea-wall at the Naval Academy, and in 1879 
he could report the value of 299,910+50 kilom- 
eters per second for the velocity of light in 

1A. A. Michelson, St. Louis Meeting of A.A.A.S., 1878; 


Am. J. Sci. 15, 394 (1878); 17, 324 (1879); 18, 390 (1879). 
2L. Foucault, Compt. Rend. 58, 501 (1862); 58, 792 
(1862). 


vacuum.’ At the same time Simon Newcomb, 
then superintendent of the Nautical Almanac 
Office, was engaged in an elaborate determination 
of the velocity of light over paths extending from 
Fort Meyer in Virginia to the Washington 
Monument and the old Naval Observatory. 
He arranged for Michelson’s transfer from the 
Naval Academy to the Nautical Almanac Office 
and for one year, 1879-1880, they collaborated in 
these experiments. Michelson’s value for the 
velocity obtained while working for Newcomb 
was 299,860+30 km/sec. This result was re- 
ported in 18914 in the complete summary of 
Newcomb’s researches on the velocity of light. 
From 1880 until 1882 Michelson studied and 
carried on researches in Europe, principally at the 
Universities of Heidelberg and Berlin and at Le 


Fic. 1. ALBERT A. MICHELSON in 1881. 


3A. A. Michelson, Astr. Papers, U. S. Nautical Almanac 
Office, I, 141 (1882). 

4Simon Newcomb, Astr. Papers, U. S. Nautical Almanac 
Office, TI, 201 (1891). 
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Collége de France in Paris. He was enabled to 
resign his commission in the Navy and go to 
Europe because he was assured by the trustees of 
the newly established Case School of Applied 
Science that he would be made the first professor 
of physics at that institution. In fact he was 
supplied certain funds for travel and for the pur- 
chase of the physical laboratory equipment 
needed at Case. During this period abroad 
Michelson invented his famous interferometer. 
The first model of this instrument was built in 
Berlin in 1881 by Schmidt and Haensch with 
funds provided by Alexander Graham Bell. This 
instrument was demonstrated by Michelson in 
Paris before Cornu and other leading French 
scientists. Cornu was very outspoken in his 
skepticism of Michelson’s explanation of the per- 
formance of the instrument, asserting that the 
interference fringes were due to some phenom- 
enon analogous to a Lloyd mirror effect in the 
half-silvered plate. However, Michelson demon- 
strated that the fringes disappeared when a piece 
of glass was placed in one of the optical paths 
beyond the half-silvered plate, and Cornu was 
convinced. 

While in Europe in 1881, Michelson received 
his official appointment as professor of physics at 
Case School of Applied Science and was granted 
leave of absence to continue his researches in 
Europe during the academic year 1881-1882. 
During this period he made the first unsuccessful 
attempt to detect the motion of the earth through 
space by means of his interferometer. Several 
factors contributed to this result. Michelson had 
been led to try the experiment upon reading a 
suggestion by Clerk Maxwell’ that the motion of 
the earth might be detected by measures of the 
velocity of light if sufficiently refined techniques 
were available. Michelson’s new interferometer 
provided just the experimental method needed. 


5 James Clerk Maxwell; article on Ether in Vol. VIII of 
9th Ed. of the Encyclopedia Britannica: “‘ . . . If it were 
possible to determine the velocity of light by observing the 
time it takes to travel between one station and another on 
the earth’s surface, we might, by comparing the observed 
velocities in opposite directions determine the velocity of 
the aether with respect to these terrestrial stations. All 
methods, however, by which it is practicable to determine 
the velocity of light from terrestrial experiments depend on 
the measurement of the time required for .the double 
journey from one station to the other and back again, and 
the increase of this time on account of the relative velocity 
of the aether equal to that of the earth in its orbit would be 
only about one hundred millionth part of the whole time of 
transmission, and would therefore be quite invisible.” 
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However, Michelson seems to have expected a 
first-order effect due to the earth’s speed relative 
to that of light, and so used only a single transit 
and return of the light beams along the inter- 
ferometer arms which were 120 cm in length. Two 
trials of the experiment were made: the first at 
the Physikalisches Institut of the University of 
Berlin, where vibration due to street traffic 
interfered ; and the second in a hollow space in 
the brick pier below the large telescope of the 
Astrophysicalisches Observatorium at Potsdam, 
where the experiment gave a negative result.® 

In September of 1882 Michelson assumed his 
duties at Case in Cleveland and here began a 
notable series of researches which firmly estab- 
lished his position in physics. During this first 
college year of 1882-1883 he made a precise de- 
termination of the velocity of light, the accuracy 
of which was not improved upon until his own 
experiments were completed in California in 
1927. This determination was made at the sug- 
gestion of Simon Newcomb’ who obtained finan- 
cial support for the work from the Bache Fund of 
the National Academy of Sciences. An optical 
path nearly } mile long was surveyed along the 
north side of the tracks of the New York, Chicago 
and St. Louis Railroad (Nickel Plate) at the rear 
of the Case campus. This survey was made by the 
Case Professor of Civil Engineering, John Eisen- 
mann, along the north rail of the track and offset 
to two concrete piers upon which the optical 
equipment was mounted. The base line was meas- 
ured as 2049.5253 feet. The arc light, dynamo 
and gas engine used to obtain a powerful light 
source were supplied by the inventor Charles F. 
Brush, of Cleveland. The rotational speed of the 
mirror was determined with the aid of an elec- 
trically driven Koenig tuning fork of nominal 
frequency 128 cycles/sec. Michelson made 23 
complete determinations of the velocity of light 
in the period from October 12 until November 
14, 1882 and gave as the mean value 299,853+60 
km/sec for the velocity of light reduced to 
vacuum.’ This result was accepted as the best 
available measure for forty-five years. 

6A. A. Michelson, Am. J. Sci. 22, 120 (1881); Phil. Mag. 
13, 236 (1882); Am. J. Sci. 23, 395 (1882); J. de Physique 
1, 183 (1882). 

7 Simon Newcomb, Astr. Papers, U. S. Nautical Almanac 
Office, II, 235 (1891). 

8 Report of John Eisenmann to A.A.M., Oct. 14, 1882. 


9 A. A. Michelson, Astr. Papers, U. S. Nautical Almanac 
Office, I, 243 (1891); Am. J. Sci. 31, 62 (1886). 
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During Michelson’s second year at Case, 1883- 
1884, he measured the velocity of light in distilled 
water. The optical path between mirrors in this 
experiment was reduced to 17.63 feet, of which 
10.03 feet were in the tube containing the distilled 
water. The result expressed as the ratio of his 
determination of the velocity in vacuum to the 
velocity in distilled water was found to be 1.33." 
At that time this result was received as strong 
evidence for the wave theory of light. Foucault 
and also Fizeau and Breguet, at the suggestion 
of Arago, had previously shown qualitatively 
that light traveled more slowly in water than in 
air,” but no quantitative measurements were 
available until those of Michelson. 

Michelson next filled the tube with carbon 
disulphide and determined the velocity of white 
light as well as the velocity of red and blue light in 
this medium. Here he obtained results which he 
had in no way anticipated. The ratio of the 
velocity in vacuum to that in CS, he found to be 
1.758+0.02, whereas the index of refraction was 
known (Verdet) to be 1.64. He also found that 
the ratio of the velocities of red and blue light 
was 1.014, the first case he had found of dis- 
persion. In his report to Simon Newcomb” he 
wrote “. . . it is with great diffidence that I 
state that the velocity of light in carbon di- 
sulphide is to that in air as 1.00 to 1.76 .. .” 

These results, however, received an immediate 
explanation in terms of Lord Rayleigh’s® theory 
which predicted the difference between the group 
velocity, which Michelson’s method measured, 
and the phase velocity given by measurements of 
the index of refraction. These results attracted 
notice throughout the world of physics and 


10 Letter of A. A. Michelson to Simon Newcomb, Aug. 15, 
1883, printed in Astr. Papers, U. S. Nautical Almanac 
Office, II, 248 (1891). 

uF, Arago, Compt. Rend. 7, 954 (1838); H. Fizeau and 
L. Breguet, Compt. Rend. 30, 773 (1850); L. Foucault, 
Compt. Rend. 30, 551 (1850). 

#2 A, A. Michelson, Astr. Papers, U. S. Nautical Almanac 
Office II, 256-258 (1891). 

1% Lord Rayleigh; Nature 24, 382 (1881); 25, 52 (1881); 
J. Willard Gibbs, Am. J. Sci. 31, 62 (1886). Gibbs also 
showed that the group velocity measured by Michelson in 
distilled water should lead to a velocity ratio of 1.35 
instead of 1.33 found by Michelson. However, the experi- 
mental uncertainty in Michelson’s determination appears 
to be +0.03 so that his result is in agreement with theory. 
It was not until 1944 that any notable improvement was 
made on these early measurements of Michelson’s. See 
R. A. Houstoun, Proc. Roy. Soc. of Edin. 62, 58 (1944). 
See also a very complete discussion of this subject by 
Raymond T. Birge, Lond. Phys. Soc. Reports on Progress 
in Physics 8, 93-94 (1941). 
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Fic. 2. Michelson-Morley interferometer used in Cleve- 
land in 1887. The illustration is taken from The American 
Journal of Science, Series 3, 34, 337 (1887). This journal is 
frequently known as “‘Silliman’s Journal,” after its founder, 
Benjamin Silliman, Professor of Chemistry at Yale 
University. 


enhanced Michelson’s reputation as a leader in 
experimental optics." 

In the year 1885 began the famous collabora- 
tion between Michelson and Professor Edward 
W. Morley of Western Reserve University which 
continued until Michelson left Case four years 
later. Their first experiment was to measure the 
velocity of light in moving water to test the 
theory of Fresnel which predicted that the ve- 
locity would be altered by a fraction (m?—1)/mn? 
of the velocity of translation of the medium, 
where 1 is the index of refraction. The expected 
fraction for water was 0.438. For this experiment 
an adaptation of the Michelson interferometer 
was employed and the measurement recorded the 
difference between the velocity in the column 
of water when in motion and at rest; the experi- 
mental technique being far superior to that em- 
ployed by Fizeau® in 1851. Their experimental 
result gave a factor of 0.434+0.02 in remarkable 
agreement with theory.'® 

Following this work Michelson and Morley 
performed the interferometer experiment which 
bears their name and upon which their fame 
chiefly rests. H. A. Lorentz'? had shown that 
the effect to be expected in this experiment should 
be a second-order one, which partly explained 
why Michelson’s Berlin experiments designed to 
measure the larger first-order effect were failures. 
Michelson and Morley accordingly designed and 
built a large multiple path interferometer, with 
the best obtainable optical parts mounted on a 

\ British Association Report, 654 (1884). 

18H. Fizeau, Annales de Chimie 57, 385 (1859). 

6 A. A. Michelson and E. W. Morley, Am. J. Sci. 31, 377 


(1886). 
17H. A. Lorentz, Arch. Néerl. 21, 103 (1886). 
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rigid stone slab 150 centimeters square and 30 
centimeters thick which was floated on mercury. 
A perspective drawing of the Michelson-Morley 
interferometer is shown in Fig. 2. Light from an 
argand burner (not shown) was divided into two 
pencils which were reflected from four mirrors 
mounted at each corner of the stone. The re- 
united pencils were observed in a_ telescope 
mounted on the nearby left-hand corner of the 
stone as shown in the illustration. A wooden 
cover protected the interferometer from air cur- 
rents and too-rapid temperature changes. With 
this apparatus having effective optical paths in 
each arm of 1100 centimeters, they expected to 
find an easily measurable shift in the position of 
the white light interference fringes when the in- 
strument was rotated through an angle of 90° in 
the horizontal plane. The final observations with 
this equipment were made in July of 1887 and 
contrary to all expectations no shift in the inter- 
ference fringes was observed.' At first the signifi- 
cance of this now famous experiment was not 
appreciated and many early attempts at its 
explanation are now known to be irrelevant. 
Even Michelson himself in giving an address on 
optical research as Vice President of section B of 
the A.A.A.S. in Cleveland in 1888 failed to men- 
tion this result.!® In fact as late as 1900 Lord 
Kelvin in his address to the International Con- 
gress of Physics in Paris could only refer to the 
Michelson-Morley experiment as one of the “‘two 
small clouds” (the other being the lack of an 
adequate theoretical explanation for the black- 
body spectrum) which dimmed the otherwise 
brilliant sky of the new scientific century. It was 
only after the pioneer work of Lorentz and 
Fitzgerald was followed by Einstein’s develop- 
ment of the special theory of relativity in 1905 
that the Michelson-Morley experiment received 
its proper praise and clear explanation. 

The following year 1888-1889 was Michelson’s 
last at Case and during this period he began his 
study of the profiles of spectrum lines with the 
interferometer, having already made the very 
important discovery that the red hydrogen line 
is a doublet.” He also laid plans to calibrate the 
standard meter in terms of light waves. In fact 


18 A. A. Michelson and E. W. Morley, Am. J. Sci. 34, 333 
(1887); Phil. Mag. 24, 449 (1887); J. de Physique 1, 444 
(1888); Sidereal Messenger 6, 306 (1887). 

19 Article in The Cleveland Leader, Aug. 16, 1888. 

2 A. A. Michelson and E. W. Morley, Am. J. Sci. 34, 
427 (1887). : 
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the principal part of his A.A.A.S. address referred 
to above” concerned ‘‘A Description of Apparatus 
for Making the Light Wave the Standard of 
Length.” This work was also carried on in col- 
laboration with Professor Morley and plans were 
made to measure the standard meter in Paris in 
terms of light waves by the interferometer tech- 
nique. The selection of the red cadmium line 
followed the work in Cleveland, but it was there 
that the “‘visibility of the interference fringes”’ 
was established as the criterion for the sharpness 
of spectrum lines.”! The final work on the stand- 
ard meter Michelson accomplished at Paris at the 
International Bureau of Weights and Measures 
in 1892-1893. 

Michelson resigned from the Case Faculty in 
1889 and went first to Clark University and then 
to The University of Chicago. During his long 
career at the latter institution he added to his 
fame by extensive researches in spectra and the 
fine structure of spectrum lines to which work he 
contributed not only refinements in the inter- 
ferometer techniques begun in Cleveland, but 
also the echelon spectrograph (1898), and im- 
proved methods for the production of ruled 
gratings. Also at the University of Chicago in 
collaboration with Dr. Henry Gordon Gale he 
performed the earth-tide experiment (1916) on 
the grounds of the Yerkes observatory at Lake 
Geneva, and in 1925 measured the rotation of the 
earth by interferometer observations on light 
moving in a large rectangular system of evacu- 
ated pipes near the present site of the Chicago 
airport. Finally, in his later years, which he spent 
largely in California, he improved his velocity of 
light measurements and at the Mount Wilson 
Observatory used an interferometer technique 
first attempted at Yerkes to determine the di- 
ameters of red giant stars. 

He died on May 9, 1931, the first American to 
have received the Nobel Prize in physics, the 
second president of the American Physical So- 
ciety, and the acknowledged leader in the field of 
precise optical measurements. 

I am greatly indebted to Dr. Karl K. Darrow 
who suggested that this paper be written and who 
gave valuablecriticisms, and alsoto Professor Ray- 
mond T. Birge for helpful discussions especially 
relating to the work on the velocity of light. 


21A,. A. Michelson and E. W. Morley, Am. J. Sci. 38, 
181 (1889); Phil. Mag.24, 463 (1887); Sctence 10, 86 (1887). 
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5. Temperature Measurements on 
Nonisothermal Combustion Flames 


A. Historical Introduction.—All the optical 
techniques for measuring flame temperatures per- 
mit only limited resolution of the flame zone and 
therefore yield measurements corresponding to 
effective temperatures over nonisothermal re- 
gions of varying extent. In view of this limitation 
it is apparent that auxiliary considerations are 
required to permit the correlation of temperature 
inhomogeneities with experimentally determined 
results. Preliminary experimental and theoretical 
work on this problem was reported by Griffiths 
and Awberry in 1929.%5 However, a thorough 
study of temperature determinations on complex 
flames does not appear to have been carried out 
until very recently.** Considerations of tempera- 
ture measurements on complex flames appear to 
have been largely restricted to application of the 
line reversal method. 

B. Temperature Measurements through Two 
Isothermal Flame Regions Which Are at Different 
Temperatures.—Griffiths and Awberry* reported 
that measured flame temperatures were ap- 
parently independent of the thickness of the 
flame. They attributed this observation, which 
should be approximately valid for an isothermal 
flame region, to the simultaneous variation of 
absorption and emission resulting in an apparent 
uniform temperature. For nonuniform flame zones 
the temperature is not independent of thickness 
as can be demonstrated by 
argument. *5 57, 65 

In order to illustrate the problems which arise 
when the line reversal method is used for the de- 


an elementary 


* Sections, references, and equations are a continuation 
of Part I which appeared in Am. J. Physics 17, 422 (1949). 

** This paper presents the results of one phase of re- 
search carried out at the Jet Propulsion Laboratory, 
California Institute of Technology, under U. S. Army 
Ordnance Department Contract No. W-04-200-ORD-1482. 

66 F, P. Bundy and H. M. Strong, Project Squid, Tech- 
nical Memo. No. NYU-2, April, 1948. General Electric Co., 
Project Hermes, Report No. 55275, Sept., 1948. ? 


termination of effective flame temperatures 
through nonisothermal regions, we may consider 
the simplest possible case in which the light from 
a blackbody source at temperature Ts passes 
first through a flame of thickness x;, mass absorp- 
tion coefficient k1, density p1, emissivity ¢:, at the 
temperature J, and then through a flame of 
thickness x2, mass absorption coefficient ke, den- 
sity pe, emissivity €2, and temperature 7». At the 
point of reversal for this composite flame it is 
apparent that the following relation must be 
valid : 


J(A, Ts) =J(A, Tg)e7* iP it1¢—kapaze 
+J(x, Ti)(1 — e~*ipiz1) e—kapare 
+J(A, T2)(1—e7*29272), (15) 
where use has been made of Kirchhoff’s law. The 
effective blackbody temperature 7T,,2 of the 
composite flame is defined through the relation 
J(A, Ti, 2) =J(A, T1) (1 — em *191*1/) eta vaze 
+J(A, Ts)(1—e7*292*2). 


From Eggs. (15) and (16) it follows that 


Jr, Ts)/J(a, Ti, 2) 
=[1 —exp(—ipixi— 


(16) 


Ropoxe) }. (17) 


Since 
J(A, Ts)/JO, T1,2) 
=exp{(c2/A)[(1/T1,2) —(1/Ts) }}, 


it follows from Eq. (17) that the apparent tem- 
perature T;,2 of a flame consisting of two iso- 
thermal regions depends on the thickness of the 
isothermal zones and on their optical densities 
but is apparently independent of the order in 
which the flame zones are arranged. The inter- 
pretation of this conclusion requires some caution 
since it can be seen from Eq. (16) that the defini- 
tion of T;,2 itself depends on the order in which 
the flames are arranged and that therefore the 
effective temperature measured experimentally 
does depend on the relative order in which the 
isothermal regions are arranged. Thus, Griffiths 
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and Awberry** found that the observed reversal 
temperature approximated the temperature of 
the hotter region if the light from the comparison 
radiator passed first through a cool flame and 
then through a relatively hot flame. With the re- 
verse arrangement the observed temperature 
tended to be intermediate between the tempera- 
tures of the two flames. 

Recently Bundy and Strong® have devised a 
method for the reconstruction of the temperature 
pattern of a flame composed of two isothermal 
regions. One of their proposed methods involves 
application of the reversal equations to one of the 
component isothermal zones and to the com- 
posite flame examined as a whole. Thus the re- 
versal temperature of the blackbody comparison 
source is Ts; at the point of reversal of the iso- 
thermal flame region whose temperature is 71. 
Hence 


J(A, Ts1) =a J(r, Ti) +J(A, Ts1)(1—an), (18) 


where a1=€,=1—exp(—:pix1) represents the 
absorptivity (which is equal to the emissivity by 
Kirchhoff’s law) of the flame whose temperature 
is T;. Similarly, at the point of reversal for two 
isothermal flame regions of temperatures 7; and 
T2, respectively, the comparison source reads a 
temperature 7's;,2. Therefore, 


J(r, Ts, 2) =aeJ(X, T2)+a,J(A, T,)(1 — ae) 
+J(A, Ts1,2)(1—a1)(1—a@e2), (19) 


where a2=1—exp(—kep2x2). Using the relation 
defining 71,2 as given in Eq. (16) it can be seen 
that Eq. (19) reduces to 


J(A, Ts1,2)(a1+a2—a1a2) = J(d, 71,2). 
From Eqs. (18) and (20) it follows that 


J(A, Ts1)/J(A, Ts, 2) 
=[LaiJ(A, T1)/J(A, T1,2) (1 —aetar/ar). (21) 


By the use of Wien’s equation, Eq. (21) can be 
rewritten as 


exp { (¢2/A)((1/T si, 2) — (1/Ts1) + (1/T)) 
—(1/T1, 2) ]} =a1(1—a2)+ a2. 


The quantity a: is given from Eq. (16) as 
aiJ (A, T1) —J(A, T1, 2) 
ayJ (x, T,)-—J(, T2) 


(20) 


(22) 


ag= 


(23) 
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From the condition of reversal for a single iso- 
thermal flame zone it has been shown (see Sec. 3) 
that 


Ts1 = Ti, 
whence Eq. (22) becomes 


exp { (c2/A)[(1/T si, 2) — (1/T1, 2) J} 
=a;(1—a2)+a2. (24) 


By using the reversal condition for the flame 
whose temperature is T; it can similarly be shown 
that 


exp { (c2/A)[(1/Ts1,2) —(1/T1, 2) J} 
=ar(1 —ai)+ai. (25) 


It is apparent from Eq. (23) and Eqs. (24) and 
(25) that the temperature JT: can be determined 
from the experimentally determined values of 
T s1,2 provided T; and a: or a2 are known. In this 
manner it is, therefore, possible to ascertain the 
temperature of one of two isothermal regions 
from measurements on the other isothermal 
region and on the composite flame.®* In practice 
it is possible that JT, may be so low that a 
thermocouple measurement is feasible, thereby 
supplementing the temperature determination of 
the composite flame by the reversal method. 
Strong and Bundy® have made a thorough 
study of the application of intensity measure- 
ments on the sodium D-line contour to a determi- 
nation of absorptivities and reconstruction of the 
spectral brightness of two isothermal flame zones. 
This experimental procedure permits calculation 
of the unknown temperature 7: from the simple 


expression 
i+—- =) 
1 


=exp{(c2/A)[(1/T 1,2) — (1/71) J}. 


The procedure of Strong and Bundy requires a 
rather elaborate experimental arrangement since 
it involves a quantitative study of line intensity 
measurements on the contours produced by self- 
absorption of the sodium D-line during passage 
through the cooler flame region. 


a,J (Xr, Ti) 
J(A, T1,2) ( 


(26) 


66 In principle it is, of course, possible to extend this 
technique of the study of complex flames with the result 
that the temperature of the mth zone can be inferred from 
measurements on the composite flame if the temperatures 
of n—1 zones have been determined independently. 
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It is apparent from the foregoing discussion 
that the difficulties of determining even the tem- 
perature of one of two isothermal regions without 
making a direct determination on this region are 
such as to limit seriously the practical value of 
measurements of this type since flames in which 
two well-defined isothermal regions exist may be 
difficult to approximate. 

For a description of the experimental arrange- 
ment used for the determination of temperatures 
on composite flames and for a more complete 
discussion of the theory the reader is referred to 
the published literature. Bundy®® has made a 
theoretical study of line widths of the sodium line 
in order to facilitate interpretation of experi- 
mental results. The interpretation of temperature 
measurements made by use of the two-path 
method (see Sec. 6) if the flame consists of two 
concentric spherical isothermal regions has also 
been considered. 

It may be of interest to note in this connection 
that detailed studies of the effect of line width on 
experimental results were made some 40 years 
ago.®7 A recent and thorough treatment of the 
self-absorption of spectral lines may be found in 
an article by Cowan and Dieke.® 


6. The Two-Path Method 

The two-path method for the measurement of 
flame temperatures is an optical technique which 
permits temperature measurements on systems in 
which the intensity of radiation emitted varies 
rapidly (and aperiodically) with time. Flame 
temperature measurements by use of the tech- 
niques described previously are generally unsuit- 
able for applications of this type. For example, 
the reversal methods could be employed only if 
the comparison radiator also emits radiation the 
intensity of which varies with time in the same 
proportion as in the flame under study. The two- 
path method has been described by Hottel and 
Broughton*’ and is finding application at the 
present time to temperature determinations on 
pulse-jet flames and other systems in which rapid 
temperature fluctuations occur.® 7° 


| hs 3 Ladenburg and F. Reiche, Ann. der Physik 42, 181 
1913). 


68 R. D. Cowan and G. H. Dieke, Rev. Mod. Physics 20, 
418 (1948). 

69 J. H. Hett, Project Squid, Technical Memo. NYU-2, 
April, 1948. 

70J. A. Sanderson, J. A. Curcio, and D. V. Estes, 
Project Squid, Technical Memo. NYU-2, April, 1948. 


493 


The two-path method depends upon a com- 
parison of (spectral) brightness temperatures 
when two different path lengths through a flame 
are viewed separately. In practice it may prove to 
be desirable to make measurements at wave- 
lengths corresponding to emission from an added 
substance which is not normally present in the 
flame. The sodium D-line has been employed for 
this purpose.®® Other investigators recommend 
use of the 4.4u-emission band of carbon dioxide”® 
or of a narrow wavelength band in the visible 
region’? which can be isolated through the use of 
interference filters’ 7? and detected by means of 
photocells. It is evident that for the two-path 
method, as for other optical methods for measur- 
ing flame temperature, it is essential that the 
substance whose radiant emission is studied is in 
statistical equilibrium with the random transla- 
tional velocity distribution of the gas molecules. 

In the two-path method the path length can be 
conveniently doubled by placing a mirror at one 
end of the field of view and making appropriate 
corrections for the reflectivity of the mirror. It is 
apparent that a rather uncertain correction for 
(multiple) reflection from the chamber walls 


- should also be introduced. If the motor wall acts 


as a cool blackbody, then the radiant energy 
emitted by a path length L of the combustion 
region, neglecting transmission losses and the like 
through windows, in the wavelength interval be- 
tween \ and A+d\, is 


J(r, Tori) dd = J(A, Tr){1 —exp( —kypL) Jd, (27) 


where J(A, Tori) is the observed radiant energy, 
Tor1 the apparent brightness temperature corre- 
sponding to a path length L, J(A, Tr) the radiant 
energy emitted by a blackbody at the true flame 
temperature 7, k, the spectral mass emissivity, 
and p the density of the gaseous emitters. When 
the path length is doubled by placing a mirror at 
one end of the field of view, the observed radiant 
intensity corresponding to the apparent bright- 
ness temperature 7J>,2 is 


J(d, Tor2)dd = J(A, Tr) [1 —exp(—RypL) ] 


[1+ R, exp(—kyapL) ]dd, (28) 


where R, is the spectral reflection coefficient of 


mJ. A. Anderson and R. L. Ueberall, NRL Report 
H-2731, January, 1946. 

7 L. N. Hadley and D. M. Dennison, J. Opt. Soc. Am. 37, 
451 (1947). 
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the front-surface reflecting mirror. From Eggs. 
(27) and (28) it follows that 


J(A, Tor2)/J(A, Tort) =1+Raexp(—RapL). (29) 


From Eq. (29) it can be seen that the intensity 
ratio J(A, Tor2)/J(A, Tori) approaches 1+R, as 
kypL approaches zero whereas the intensity ratio 
tends toward unity as k,pL becomes large. Refer- 
ence to Eq. (27) shows that 


exp( —kypL) =1-— J(x, Tort) /J(A, Tr), 
whence 


1 —J(x, Tor2)/J(A, Tort) 


= —R,{1 —J(A, Tors) /J(A, Tr) ], (30) 


or 


{1—exp[(—¢2/A)(1/Tor2—1/Tor1) J} 
oad —R,{ 1 —exp[(—¢2/A)(1/Torn—1/T r) ]}, (30a) 


where ¢2 is the second radiation constant and the 
radiant intensities J(A, T) have been replaced by 
the use of Wien’s equation. It is apparent from 
Eqs. (30) and (30a) that the experimental de- 
termination of J(A, Tor) and J(A, Tir2) permits 
calculation of J(A, Tr) and, therefore, determines 
the gas temperature Tr uniquely according to 


Planck’s (or, to a close approximation, Wien’s) 
equation. 

The preceding summary of the principles in- 
volved in flame temperature measurements by 
use of the two-path method has been restricted to 
an isothermal region as source of radiation. As is 
also the case with other optical methods for flame 
temperature measurements, considerable com- 
plication is introduced whenever measurements 
are made on nonisothermal regions (see Sec. 5). 
Bundy® has considered the utility of the two- 
path method if the flame on which measurements 
are made consists of a hot flame zone surrounded 
by a cooler region. His idealized calculations led 
him to conclude that the two-path method would 
yield reasonable results (corresponding closely to 
the actual temperature of the hot central core of 
the flame) if the product k,)pZ was small enough 
so that the intensity ratio J(A, Tor2)/J(A, Tori) as 
given by Eq. (29) would be of the order of v2 if 
R, is equal to unity. 


7. The Compensated Hot-Wire Method 


A. Historical Introduction.—A compensated 
thermocouple technique for flame temperature 
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measurements was employed by Berkenbusch 
in 1899 at the suggestion of W. Nernst.7* 
Berkenbusch’s method effectively eliminates the 
troublesome radiation correction required for 
flame temperature measurements made with 
thermocouples. Schmidt** replaced the thermo- 
couple used by Berkenbusch by a hot-wire device. 
A method practically identical with that used by 
Schmidt was employed some years later at the 
National Physical Laboratory’! and also by 
Griffiths and Awberry,*® who do not appear to 
have been aware of the earlier work of Schmidt. 
Schmidt,** Griffiths and Awberry,*® and Kohn*’ 
found excellent agreement between results ob- 
tained by the compensated hot-wire technique 
and by the line-reversal method. More recent 
discussions of the compensated hot-wire method 
have been given, for example, by Lewis and von 
Elbe* in a review article on the measurement of 
flame temperature. 

B. The Hot-Wire Method.—The hot-wire 
method for the measurement of flame tempera- 
tures is included in the present survey because it 
is probably the simplest and‘ most useful non- 
optical technique for checking or supplementing 
results obtained through application of one of the 
optical methods. The hot-wire method is based 
on an energy balance for a radiating wire placed 
in a vacuum or exposed to a flame. Schmidt*® 
heated an extended wire electrically ina vacuum, 
and determined its temperature with an optical 
pyrometer. The rate of loss of energy by radiation 
per unit length of wire E was measured simul- 
taneously in order to determine the functional 
relation between £ and the wire temperature 7. 
Next the rate of energy input £’ into the wire 
required to raise it to the temperature T” while it 
was exposed to the flame zone was determined. 
At the point of intersection of a plot of E vs. T 
with a plot of EZ’ vs. T’, it is evident that T=T’ 
and E=E’. But the condition T=T’ and E=E’ 
assured the existence of thermal equilibrium’ 
between the wire and the flame since the wire 
was losing energy by radiation at the same rate in 
the flame as in a vacuum while it was at the same 
temperature in both. 

78 F, Berkenbusch, Ann. der Physik u. Chem. 67, 649 (1899). 

™ Report for 1926, p. 63. 

75 Equilibrium exists only if the opacity of the flame for 


the energy radiated by the wire is negligibly small. This 
condition seems to be fulfilled in practice in many cases. 
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The method of Griffiths and Awberry*® is 
identical with that used by Schmidt except that 
the current input into the wire is measured in- 
stead of the total energy lost by radiation. If I 
represents the current required to heat the wire 
to the temperature T in a vacuum, and J’ is the 
current required to heat the wire to the tempera- 
ture T’ in the combustion zone, then the condi- 
tion J=J' at T=T’ again assures the existence of 
thermal equilibrium between the wire and the 
flame. Since it is easier to measure the current 
input into a wire than to determine the energy 
lost by radiation from a wire, it seems preferable 
to employ the experimental arrangement of 
Griffiths and Awberry rather than that of 
Schmidt. 

C. Limitations of the Compensated Hot-Wire 
Method.—It is apparent that the compensated 
hot-wire method will provide a measure of the 
true gas temperature only if the wire does not by 
its presence disturb the combustion processes. 
Furthermore, the wire must retain the same 
emissivity per unit length in the flame as. it 
possessed in a vacuum. That the wire does, in 
general, cause some disturbance of the combus- 
tion process has been recognized for some time.*” 
Pitting of the wire surface also appears to be a 
source of error in some cases.*® The hot-wire 
technique is thus seen to involve some of the 
limitations of thermocouple measurements and 
probe techniques. As is the case with thermo- 
couple measurements, a point-by-point explora- 
tion of the flame region can be made. Compared 
with thermocouple measurements, the compen- 
sated hot-wire technique has the advantage of 
eliminating the radiation corrections which must 
be made whenever thermocouples are used at 
elevated temperatures. Furthermore, since only a 
single wire is used, no error can be introduced by 
preferential catalysis of the combustion process 
by one or the other of the component wires 
of a thermocouple as used ordinarily or of 
the compensated thermocouple introduced by 
Berkenbusch.?* 


8. Estimation of Flame Temperature from the 
Change in Spectral Line Intensity with 
Quantum Number 


A. Historical Introduction.—Flame tempera- 
ture measurements utilizing the change in in- 
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tensity of spectral lines with quantum number 
were introduced by Ornstein and his collaborators 
in 1928. Ornstein and van Wyk?*&78 used the 
O—O band of the negative nitrogen bands ob- 
served during electric discharges in order to de- 
termine the rotational temperature. The determi- 
nation of the rotational temperature yields, of 
course, a significant result only if the rotational 
energy levels are in statistical equilibrium with 
the random translational molecular energies or at 
least behave as though this statistical equilibrium 
were maintained. Lindh” showed that the distri- 
bution of rotational levels in nitrogen corre- 
sponded to statistical equilibrium during elec- 
tronic collision, in agreement with earlier findings 
by Harries*® who noted that electron impact does 
not appreciably alter the rotational distribution. 
Elliott® found a Boltzmann distribution of rota- 
tional levels corresponding to the temperature of 
active nitrogen during activation with active 
nitrogen. Temperature determinations based on 
the change in observed intensity for the vibra- 
tional levels of the *2 2y-bands of cyanogen 
during an electric discharge were found by 
Ornstein and Brinkman® to be in agreement with 
rotational temperatures. Temperature determi- 
nations on electric arcs using the rotational in- 
tensity variation have also been made on the 
spectra of gold** and hydrogen* with good re- 
sults. The early work on temperature measure- 
ments in arcs has been reviewed by Ornstein*®® 
who, incidentally, calls attention to the fact that 
different rotational temperatures were obtained 
for the bands of cyanogen and aluminum oxide, a 
fact which he attributes to emission from differ- 
ent positions in the arc (which are at different 
temperatures) for the two emitting, molecular 
species. 


76 L. S. Ornstein and W. R. van Wyk, Zeits. f. Physik 49, 
315 (1928). 
( 77 W. R. van Wyk, Proc. Amsterdam Acad. Sci. 32, 1243 
1929). 
78G. Herzberg, Molecular spectra and molecular struc- 
ture I, diatomic molecules (Prentice-Hall, Inc., New York, 
1939). 
79 A. E. Lindh, Zeits. f. Physik 67, 67 (1931). 
80 W. Harries, Zeits. f. Physik 42, 26 (1927). 
81 A, Elliott, Zeits. f. Physik 67, 75 (1931). 
QL. S. Ornstein and H. Brinkman, Proc. Amsterdam 
—. - 34, 33 (1931). 
Ornstein and S. Sambursky, Proc. Amsterdam 
Acad. “Soi 34, 339 (1931). 
Ss. Ornstein, J. G. aye and J. Wouda, Proc. 
haute Acad. Sci. 34, 505 (19. 
85 L. S. Ornstein, Physikalische Peis 32, 517 (1931). 
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Shortly after the pioneering work of Ornstein 
and his collaborators had been reported, it be- 
came apparent that the rotational line intensity 
method may lead to erroneous values for the 
temperature because of abnormal rotational 
energy distributions. Gaviola and Wood** found 
an abnormal population of the higher rotational 
energy levels of HgH as the result of excitation 
by Hg, a phenomenon which has since been 
studied extensively by Rieke*? who observed a 
more nearly normal rotational distribution of 
HgH in the presence of OH than in the presence 
of Nz. Abnormal rotation appears in the negative 
bands of CO and the negative bands of N2*+ ex- 
cited in a low voltage arc,** whereas a more nearly 
normal distribution of energy is observed for the 
same bands in the electrodeless discharge and in 
a neon transformer glow discharge. Abnormal 
excitation of rotational levels has been observed 
for C. and CH during electric discharges through 
acetylene,® for OH during an electric discharge 
through water,*® and for CuH in an electric arc.*! 
The abnormal rotation of OH has been the sub- 
ject of extensive study by Oldenberg* who found 
normal rotation of OH formed in an electric 
discharge through H.O during absorption and 
noted that the abnormal rotation observed 
during emission could be greatly reduced by the 
addition of large amounts of helium. More re- 
cently Lyman*®* has made temperature determi- 
nations using the 2 *TI-transition of OH in 
oxygen-hydrogen flames and in copper arcs. She 
found a thermal distribution of intensities from 
which temperatures of 3275°K and 6350°K, re- 
spectively, were calculated. The OH spectrum in 
a vacuum tube discharge showed abnormal rota- 


86 FE. Gaviola and R. W. Wood, Phil. Mag. 6, 1191 

toaey E. Gaviola, Phil. Mag. 6, 1167 (1928) ; ibid. 6, 1154 
1 . 

81 F. F. Rieke, Physical Rev. 42, 587 (1932); J. Chem. 
Physics 4, 513 (1936) and 5, 831 (1937 » 

80. S. Duffendack, V. Revans, and A. S. Roy, 
Physical Rev. 45, 807 (1934). 

89 W. Loch te-Holtgreven, Zeits. f. Physik 64, 443 (1930) ; 
ibid. 67, 590 (1931). 

K. F. Bonhoefer and T. G. Pearson, Zeits. f. Physika- 
pr, Chem. (B) 14, 1 (1931). 

* H. Schiiler and H. Gollnow, Zeits. f. Physik 108, 714 
(1938) ; ibid. 109, 432 (1938). The anomalous excitation of 
CuH has recently been ascribed to a catalytic effect pro- 
duced by the surface of the copper cathode: L. Reinebeck, 
Zeits. f. Naturforschung 2a, 251 (1947). 

% OQ. Oldenberg, Physical Rev. 46, 210 (1934); J. Chem. 
Physics 3, 272 (1935). 
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tion which was found to be dependent on the 
pressure of water vapor in the tube. A comparison 
of spectroscopically determined rotational tem- 
peratures in arcs and of thermometric tempera- 
tures has been given by Horie. 

As is evident from the preceding literature 
survey of some of the work on temperature de- 
terminations by means of the rotational or 
vibrational line intensity variation, most of the 
applications of this technique have been con- 
cerned with studies on arcs. Oldenberg” has given 
a discussion of some of the factors responsible for 
the production of abnormal rotation in molecules 
which will be considered later on with reference 
to their applicability to combustion flames. The 
occurrence of abnormal vibrational or rotational 
excitation in combustion flames cannot be ruled 
out a priori in view of some direct experimental 
findings of vibrational (and possibly of rotational) 
time lag of O2 and Ne during explosions in a 
spherical vessel.°* More recently Gaydon noted a 
significant change in the rotational temperature 
of OH produced as the result of the addition of 
acetylene.*® 

B. Principles Involved in Flame Temperature 
Measurements Based on the Variation in the In- 
tensity of Spectral Lines with Quantum Number.— 
The principles involved in flame temperature 
determinations by use of the variation in the in- 
tensity of spectral lines are given, for example, in 
one of the well-known textbooks on molecular 
spectra by Herzberg.”* For the present purposes a 
summary of the relevant formulas for determining 
the intensity of spectral lines, similar to that 
given recently by Chwalow,*’ seems in order. 

Since quadrupole transition probabilities for 
radiation of visible light are only about 10-* 
times the dipole transition probabilities, the 
following discussion will be restricted to dipole 
radiation because of the obvious difficulty of 
observing the weak emission or absorption lines 
corresponding to quadrupole radiation. For rea- 
sons which have been discussed by Oldenberg,” 
reliable temperature measurements are more 


% T. Horie, Proc. Physico-Math. Soc. of Japan (3), 22, 
787 (1940). 


%°B. Lewis and G. von Elbe, Combustion flames and 
= of gases (Cambridge University Press, 1938). 
G. Gaydon and H. G. Wolfhard, Proc. Roy. Soc. 
A194 165 (1948). 


E. Chwalow, Picatinny Arsenal Report No. 
1666, October, 1947. 





DETERMINATION OF FLAME TEMPERATURES 


difficult to make by the use of intensity variations 
of vibrational lines than of rotational lines (see, 
however, Ref. 82 for application of the vibration 
spectrum of cyanogen). Furthermore, since the 
theoretical considerations are similar, the present 
discussion will be restricted to the change in in- 
tensity of rotational lines. 

The intensities of rotational lines during ab- 
sorption J,», and emission J.» are given, to a good 
approximation, by the following two relations, 
respectively :78 


Taps =const. X (J’+ J" +1) 

Ke-Bs'I"+helkTyR2 (31) 
I om =const. X (J’+ J" +1) 

Ke BI I'+VhelkT ys R2. (32) 
Here J’ is the rotational quantum number of the 
upper state, J’ the rotational quantum number 
of the lower state, # is Planck’s constant, c is the 
velocity of light, represents the Boltzmann con- 
stant, T is the absolute temperature, v is the 
observed frequency of the spectral lines (which is 
very nearly constant for a given vibration-rota- 
tion band), R is the matrix element corresponding 
to the given dipole transition, and B’ and B” are 
defined through the relations B’=h/(8x’cI’) and 
B" =h/(8x*cI"’) where J’ and J" are the moments 
of inertia corresponding to the upper and lower 
states, respectively. The selection rules for vari- 
ous spectral transitions are dependent upon the 
vanishing or nonvanishing of the matrix element 
R”" between two states 7 and m. The vector R™ 
is defined through the relation 


R = f vavnX(Z er;)dy, 


where y,, and y,, are the respective eigenfunctions 
of the two combining states. The star denotes the 
complex conjugate, >>;e,r; is the total dipole 
moment (i.e., the dipole moment summed over all 
of the charged particles) with components 
be=DiexX, Py=Liewsi, Po=L ies, and the 
integral is extended over the 3N-dimensional con- 
figuration space of the N particles composing the 
molecule whose rotational intensities are being 
determined.’* The squares of the matrix elements 
correspond to the classical transition probabilities 
which can be evaluated either by the use of 
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Einstein transition probabilities®** or from Van 
Vleck’s correspondence principle for absorption. 
For a more satisfactory treatment of absorption 
and emission of radiation, reference should be 
made to the work of Dirac.!” 

For a given vibration-rotation band the factors 
vR? and y*R? are very nearly independent of the 
rotational energy levels whence Eqs. (31) and 
(32) reduce to Eqs. (33) and (34), respectively : 


ie B'J'"(J"+1)he 
Repcmicinigillaomomenen, (5S 
THI" 4+1 kT 


Ta B'J'(J'+1)he 
eg lag — ——______— 
THI HA 


where C=const. XvR? and C’ = r‘C are, approxi- 
mately, independent of the factor J(J+1). There- 
fore, a plot of logZa./(J’+J"+1) vs. J’(J" +1) 
or of logZem/(J’+ J" +1) vs. J’(J’+1) should be 
linear with the slope yielding the value of the 
rotational temperature.”® 

Equations (33) and (34) refer to '2:>——"2- 
transitions involving states with zero electronic 
angular momentum. For a discussion of other 
transitions reference should be made to Herzberg”® 
or to the quantitative treatment of Hénl and 
London.!™ A discussion of refinements of Eqs. 
(33) and (34) taking into account such factors as 
vibration-rotation interactions can be found in 
the literature.”® 

C. Application of the Variation of Rotational 
Line Intensity with Quantum Number to Flame 
Temperature Measurements.—For temperature 
determinations in combustion gases the investi- 
gator has usually a choice of molecular species 
which can be used for temperature measurements. 
For example, in a gasoline-oxygen flame, temper- 
ature determinations can be made, in principle, 
from the spectra of OH, Cz, CH, CO, Oo, and Noe. 
The choice of molecular species becomes thus a 
matter of convenience provided the rotational 
energy distribution of the selected molecular 
species fulfills the requirement of being in sta- 
tistical equilibrium with the random translational 
energy of the gas molecules. 


’ 


98 A, Einstein, Physikalische Zeits. 18, 121 (1917); Verh. 
der Deuts. Physikalische Ges. 18, 318 (1916). 

99 J. H. Van Vleck, Physical Rev. 24, 330, 347 (1924). 

100 P, A. M. Dirac, Proc. Roy. Soc. A112, 661 (1926) ; zbid. 
A114, 243 (1927). 

101 H, Hénl and F. London, Zeits. f. Physik 33, 803 (1925). 
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Since the theoretical treatments of temperature 
measurements from absorption and from emis- 
sion spectra are similar, there seems to be little 
preference for one or the other method from this 
point of view. However, attention should be 
called to the fact that emission spectra are much 
more intense than absorption spectra and are 
therefore more easily observed. The reason for 
this is the proportionality of I,m to v4 rather than 
to vas is the case for I», (see Eqs. (31) and (32)). 

Oldenberg* has discussed abnormal rotation in 
molecules and has given a review of some of the 
processes which may be responsible for abnormal 
rotational excitation. The processes recognized 
by Oldenberg are the following : 


(1) Simultaneous dissociation and excitation 
of polyatomic molecules. The production 
of abnormally excited OH (designated as 
OH*) during electric discharges through 
H.0 is a well known example of this type 
and may be described schematically as 


energy + HLXO-H+0OH*. 


Oldenberg’s work, proving that the addi- 
tion of helium, which is known not to 
alter electronic excitation,!°? produced a 
more nearly normal emission spectrum 
of OH, would appear to indicate no per- 
sistence of abnormal rotation during 
collision.” 

Excitation of diatomic molecules may 
result in abnormal vibration in accord 
with the Franck-Condon principle. Some 
of the excess energy may then be trans- 
ferred to rotational energy levels as the 
result of impacts during the life of the 
excited vibrational state. Rieke®’ has 
shown that this process is responsible for 
abnormal excitation of HgH. Oldenberg 
suggests that the intensity change in 
certain H, bands produced by the addi- 
tion of He may have a similar origin.*” 
Excitation by impact. Oldenberg cites 
two examples, one of which involves 
preferential excitation of Nz by a canal 
ray.!04 


12 TL, von Hamos, Zeits. f. Physik 74, 379 (1932). 

103 H. Beutler and E. Rabinowitsch, Zeits. f. Physikalische 
Chem. (B) 8, 403 (1930). 

10H. D. Smyth and E. G. F. Arnott, Physical Rev. 36, 
1023 (1930). 


(4) Abnormal excitation as the result of 
perturbations. Oldenberg indicates that 
the fact that different rotational temper- 
atures were found! for the strong and 
weak lines of F. may be attributed to 
perturbations as may also the similar 
effect!°* observed for CN. 


Of the four processes listed by Oldenberg as 
possible sources of abnormal energy distribution, 
only the first two would be expected to be of 
importance for studies on combustion flames. 
Perhaps the terminology should be changed by 
replacing Oldenberg’s excitation as the result of 
dissociation by the more inclusive designation of 
excitation associated with chemical reaction. 

The vibrational energy distribution is easily 
changed as the result of excitation by collision 
because of the operation of the Franck-Condon 
principle.*? This fact forms the basis for the 
statement made by Pearse and Gaydon" that 
band systems arising from states in which the 
configuration of the nuclei differs markedly from 
that of the normal state may be relatively strong 
in flames but not in arcs since excitation in flames 
occurs mainly as the result of collisions. Ab- 
normal excitation of vibrational levels may then 
lead to abnormal rotational energy distribution 
in accord with Oldenberg’s mechanism (2) listed 
above. However, it should be noted that energy 
exchange by collision will ultimately be respon- 
sible for a thermal distribution of energy. It 
follows, therefore, that at elevated pressures the 
rate of exchange of energy by collision may be so 
rapid that a thermal energy distribution is main- 
tained at all times during combustion, for all 
practical purposes. 

The rate of attainment of thermodynamic 
equilibrium after the occurrence of abnormal ex- 
citation can be analyzed, approximately, in terms 
of classical physical considerations. The exchange 
of energy by collision will evidently be limited 
because of the requirements of the laws of con- 
servation of energy and momentum." Therefore, 
progress toward equilibrium will be slow for col- 


105 H, G. Gale and G. S. Monk, Astrophysical J. 69, 77 
(1929). 

106 M. Fassbender, Zeits. f. Physik 30, 90 (1924). 

107 R. W. B. Pearse and A. G. Gaydon, The identification 
of molecular spectra (John Wiley and Sons, New York, 1941). 

108Q. Oldenberg and A. A. Frost, Chem. Rev. 20 99 
(1937) ; Physical Rev. 37, 194 (1931). 
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lisions between particles of greatly different 
masses. In a combustion chamber operating at 
elevated pressures it is not to be expected that 
appreciable excitation anomalies will persist for 
any length of time because of the frequent occur- 
rence of collisions between many different mo- 
lecular species, free radicals, and ions. However, 
the possibility of abnormal excitation in regions 
of active combustion cannot be ruled out a priori 
because the balance between rates of abnormal 
excitation by reaction and collision and rates of 
equalization of energy by collision is not known 
in sufficient detail to permit useful predictions. 

At elevated pressures even the vibrational 
energy levels of some of the molecules may have 
population densities corresponding to statistical 
equilibrium and, therefore, the intensity change 
of vibrational lines may be suitable for tempera- 
ture determinations. ** 1%. 110 Nevertheless, caution 
seems to be in order regarding the application of 
the line intensity technique for temperature 
measurements, particularly in view of some ex- 
perimental findings which indicate abnormal dis- 
tribution of energy.**° In addition to the 
methods described in other sections of this report, 
valid temperature measurements can frequently 
be obtained by using line intensity variations in 
absorption spectra,®?-"3 emission spectra ex- 
cited by heat,®?44"5 or emission spectra of 
molecules such as Ne, which are not readily ex- 
cited abnormally. 

Since the rate of establishment of thermal equi- 
librium at low pressures is relatively slow, it is not 
surprising that examples of abnormal rotational 
energy distribution in low pressure combustion 
flames are encountered. A very striking example 
of processes of this type has been reported by 
Gaydon** who noted a change in the rotational 
temperature of OH of several thousand degrees 
when some acetylene was introduced into the 
flame. Discussions of excitation lags in gases 
under various other conditions have been re- 
ported in the literature. Kneser and Knudsen 


109 G, Herzberg, Zeits. f. Physik 49, 761 (1928). 

10 N. R. Tawde, Proc. Physical Soc. 46, 324 (1934). 

11H. Kuhn, Zeits. f. Phystk 39, 77 (1926). 

12 E. Hutchisson, Physical Rev. 36, 416 (1930). 

13K, Wurm, Zeits. f. Astrophysik 5, 260 (1932). 

14L. S. Ornstein, H. Brinkman, and A. Beunes, Zeits. f. 
Physik 77, 72 (1932). 
1933) Th. J. Ter Horst and G. Krygsman, Physica 1, 114 
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found excitation lags in experiments with sound 
waves and have given numerous references to 
earlier work on the same subject."® Large rota- 
tional time lags have been reported during rapid 
compressions of diatomic gases"? and for the 
discharge of steam through an orifice."* Con- 
siderable work on the general subject of excitation 
lags in gas dynamics has been reported"*!*? and 
may ultimately prove to be of use for predicting 
the molecular species which can be used for 
reliable temperature determinations. 


9. Conclusions 


The survey of optical methods for the measure- 
ment of flame temperature presented in Secs. 1 to 
8 has been concerned primarily with a review of 
important principles involved in temperature de- 
terminations. The discussion was based on the 
assumption that a temperature measurement will 
be useful, from the operational point of view, only 
if it isa direct measure of the random translational 
energy of the molecules. As explained in the 
Introduction, this point of view was adopted in 
order to make the temperature a measure of 
combustion efficiency as defined ordinarily in 
practical applications. It should perhaps be em- 
phasized that for some purposes, such as the 
investigation of rotational excitation or the for- 
mation of molecules of abnormal energy content, 
a useful ‘‘temperature measurement”’ would be a 
temperature measurement under nonequilibrium 
conditions. Thus if we were interested in de- 
termining the distribution of hydroxy] radicals in 
the various rotational energy levels, then we 
would measure the intensity of the rotational 
emission lines and would not be particularly con- 
cerned with the equality of the rotational tem- 
perature of OH and the mean translational 
temperature of all of the gas molecules. This 
problem arises only if our objective is the de- 
termination of the mean translational tempera- 


6H, O. Kneser and V. O. Knudsen, Ann. der Physik 21, 
682 (1935). 
(1939). Lewis and G. von Elbe, J. Chem. Physics 7, 197 
(1939. V. Korvin-Kroukovsky, J. Franklin Inst. 227, 99 

u9., Landau and E. Teller, Physikalische Zeits. der 
Sowjetunion 1, 34 (1936). 

120 A, Kantrowitz, J. Chem. Physics 14, 150 (1946). 

121 P. W. Huber and A. Kantrowitz, J. Chem. Physics 15, 
275 (1947). 

12C, Zener, Eleventh Report of the Committee on 
Contact Catalysis, National Research Council, p. 103, 1935. 
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ture. It is evident that any optical measurement*“dence to the contrary, we are therefore justified 


of flame temperature will yield definite and 
useful information concerning some aspects of the 
flame and the gases. 

A middle-of-the-road course has been followed 
throughout the present discussion as regards the 
experimental measurement of the translational 
temperature. Thus attention has been called to 
the fact that erroneous temperature readings are 
obtained in some cases. At the same time an 
effort has been made to emphasize that the defi- 
nition of the translational temperature retains its 
usefulness under all conditions and that this 
temperature can be measured if not with one 
method then with another. Unfortunately, it is 
not possible to give a concise prescription for 
evaluating the translational temperature under 
all conditions. An attempt to present a conclusion 
of this type would be inconsistent not only with 
the thesis of the present survey but would also 
conflict with some experimentally observed 
results. 

We do not have a priori information concern- 
ing the approach to statistical equilibrium in a 
given combustion flame. In the absence of evi- 


in choosing any optical technique for which the 
equipment is available or which appears particu- 
larly interesting for other reasons. In many cases 
(e.g., the flame of a Meker burner) our first tem- 
perature measurement will probably determine 
the translational temperature. In general, how- 
ever, it is necessary to verify the equilibrium 
assumption by demonstrating that concordant 
results can be obtained by the use of two or more 
independent techniques. Detailed knowledge of 
the experimental findings of other investigators 
will permit the intelligent exclusion of some 
methods under certain conditions. For example, 
for a temperature determination in a region of 
active combustion we might prefer one of the 
reversal methods to a measurement based on the 
change of rotational line intensity with quantum 
number while, at the same time, we would expect 
absorption-emission pyrometry to be less likely 
to give large errors than the line reversal method. 
Nevertheless, it is entirely possible that the 
translational temperature of the flame under 
study can be determined equally well by all three 
methods. 
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HE determination of the normal coefficient 

of absorption and acoustical impedance of 

materials is an important and interesting type of 

experiment which may be performed in the usual 

acoustical laboratory by the method of standing 
waves in a tube. 

A loudspeaker may be used to set up standing 
waves at one end of a tight, heavy tube with the 
acoustical material at the opposite end as shown 
in Fig. 1. The normal absorption coefficient is 
found by measurement of the difference between 
the maximum H;2 and minimum [, intensity in 
decibels (i.e., H2—Z1) using a traveling pressure 
microphone. The resistive and reactive compo- 


nents of impedance are found by measurement of 
the ratio D,/Dz in addition to the above two 
measurements for the normal absorption coeffi- 
cient. The distance from the acoustical material 
to the first excess pressure minimum is D,, while 
the distance from the first excess pressure mini- 
mum to the second excess pressure minimum is Ds. 


Theory 


The following formula has been developed? 
for the normal coefficient of absorption a, when 


1S. L. Brown and M. Y. Colby, Laboratory manual, 
electricity and magnetism, light, and sound (University Co- 
operative Society, Austin, Texas, 1944), p. 45. 

2H. J. Sabine, J. Acous. Soc. Am. 14, 143 (1942). 





.) oh 


, c= CC ae VW SS 


ACOUSTICAL IMPEDANCE AND 


the tube method is used 


logio 1L/20—1\? 
an= i- (———) (1) 
logio'Z/20+1 


where L(=H2—L)) is the difference between the 
maximum and minimum intensities in decibels. 
This equation may also be written as 


dn=1—k? (2) 


where & is the pressure reflection coefficient. 

Expressions for the resistive R and reactive X 
components of the acoustical impedance have 
been derived by Davis,’ Olsen,t and Morse.§ 
They are as follows 


R 1—,? 


cc a 3 
pe 1+k?+2k cos(2rD;/D2) i 


X 2k sin(27D,/Dz) (4) 
pc 1+k?+2k cos(2eD,/D») 


where p is the density of the air, c is the velocity 
of sound in air, & is the pressure reflection coeffi- 
cient as defined in Eq. (2) and D;/Dz is the ratio 
previously defined and shown in Fig. 1. 

For air at sea level the specific acoustic resist- 
ance pc is approximately 42.7 g/(sec cm?). Thus 
Rand X are found readily from Eqs. (3) and (4), 
and the acoustic impedance is given by 


Za=RatjXa (5) 
with magnitude 


|Za| =(R?+X?)}. (6) 


Procedure 


There are two impedance tubes in the acoustic 
laboratory at William Jewell College. Both are 
nine feet long with square cross sections and 
inside areas 1.0 and 0.25 ft? respectively. The 
former is used for absorbing materials such as 
diaphragmatic tiles that may not be examined in 
smaller surface units. The latter may be used at 
much higher frequencies and is, therefore, more 
useful when materials of area 0.25 ft? may be 


3 A. H. Davis, Modern acoustics (Macmillan, New York, 
1934), p. 262. 

4H. F. Olsen, Elements of acoustical engineering (D. van 
Nostrand, New York, 1940), p. 384. 

5P. M. Morse, Vibration and sound (McGraw-Hill, New 
York, 1936), p. 192. 
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Fic. 1. Arrangement of apparatus used to determine acous- 
tical impedance and absorption coefficients. 


tested. Both tubes are constructed of well- 
seasoned two-inch lumber which was sanded, 
covered with wood filler and then given two coats 
of shellac. 

Other items of equipment include a variable 
oscillator, a high quality amplifier, speaker, small 
traveling microphone of good quality, and some 
type of output meter such as a power meter, 
voltmeter, or sound level meter. 

The apparatus is arranged as shown in Fig. 1. 
The sound source consists of a variable oscillator 
connected to the speaker through an amplifier. 
The acoustical material is attached to a backing 
material in the same way as used on the wall and 
the whole is secured tightly to the tube by the 
end plate which is bolted to the tube. Vacuum 
wax or the equivalent may be used to secure 
tightness if in doubt. The small traveling pressure 
microphone serves as a detector of the nodes and 
loops of pressure as it moves up and down the 
tube on two tiny rollers. The cross-sectional area 
of about 16.2 percent intercepted by the micro- 
phone and roller mechanism does not change the 
results appreciably. The maximum and minimum 
intensity of the standing wave pattern as meas- 
ured by the excess pressure amplitude is read on 
the output meter and their difference is given by 


L=H.—I}. 


The harmonic components of the tube may be 
computed provisionally by use of the formula, 


fm=mc/2l, 


where c and / represent the velocity of sound and 
the length of the tube respectively, and m is the 
order of the harmonic. The second harmonic 
where m=2 represents the lowest tone for which 
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Fic. 2. Resistive and reactive components of acoustical 
impedance of a single diaphragmatic acoustical tile with 
cross section of 1.0 ft?. 


resonant measurements are practical. In our tube 
the second harmonic has a frequency about 130 
cycles per second. Overtones corresponding to 
frequencies of about 260, 390 and 520 have also 
been used. 

Since there is attenuation within the tube, 
there is some increase in the intensity at suc- 
cessive minima as one measures these minima 
away from the sample. In general, Le will be 
greater than Zi, but if the tube is well con- 
structed, the error will not be serious unless high 
accuracy is desired. In this case, a method of 
correction is given in Sabine’s? article. 

The real and reactive components of the im- 
pedance may be plotted as shown in Fig. 2. The 
reactive component may decrease in magnitude 
and pass through the zero value. The real part 
may be expected to be more nearly straight ex- 
cept for factors not accounted for in the above 
theory. The absolute value of the impedance may 
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be drawn by some distinguishing curve on the 
same graph. A graph of the companion normal 
absorption coefficient coupled with the absolute 
value of the impedance gives the most useful type 
of information on acoustical absorption material. 

The experiment is important in that it links the 
ideas of the absorbing characteristics with the 
impedance characteristics of materials. Thus, the 
factor & used in the determination of the normal 
absorption coefficient is used in the calculation of 
the resistive and reactive components of impe- 
dance of the material. 

The experiment does not indicate a method for 
separating the two components of reactance, 
namely, the spring and the mass factors. Infor- 
mation on the magnitudes of these two com- 
ponents may be estimated by making certain 
simplifying assumptions especially if the ab- 
sorbing material has a marked absorption maxi- 
mum. In diaphramatic materials,*® this absorption 
maximum may be related to a minimum re- 
actance. By equating the mass reactance to the 
spring factor reactance, the frequency for maxi- 
mum absorption may be approximated. 

Absorption may be expected to be high? when 
the acoustical impedance is equal approximately 
to the impedance of the air. Material may there- 
fore be made up so as to have reactances (X) of 
zero magnitude, or near zero magnitude at the 
frequency desired for maximum absorption. The 
dissipative (R) factor may also be varied in a 
manner to approach the value of that for air. 
Usually, the problem is to design acoustical 
absorption material so as to obtain a reasonably 
constant impedance over a wide frequency range 
in order to produce a reasonably uniform absorp- 
tion coefficient. 

6° L. B. Ham and H. T. Darracott, ‘ 
istics of vibrafram, 

7H. Sabine, 


acoustical materials,” 
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It is when the pirates count their booty that they become mere thieves.—WILLIAM BOLITHO. 








T is likely that teachers of the first-year general 
course in physics would hesitate to give to 
their students the problem of solving an alter- 
nating current circuit similar to that shown in 
Fig. 1. This type of circuit would ordinarily be 
one for advanced students in electrical engi- 
neering. The equivalent direct current circuit, 
however, would not be considered too difficult 
for students of elementary physics. The purpose 
of this article is to describe a method for solution 
of some kinds of steady state a.c. circuit problems 
which is, in principle, no more difficult than the 
solution of equivalent direct current circuit 
problems. 

For the d.c. circuit equivalent to Fig. 1, the 
straightforward application of Kirchhoff’s laws 
results in eight simultaneous equations not all 
independent to be solved for a number of un- 
known quantities. The solution of these equations 
represents much more of a mathematical effort 
than is usually required of first-year physics 
students and the method of solving such circuits 
by direct application of Kirchhoff’s laws is 
usually avoided. The method of simplification of 
circuit is the one most easily learned. This method 
can be applied to a.c. circuit problems in nearly 
the same manner as for d.c. circuit problems. The 
only difference between the two solutions is that 
for a.c. circuit problems additions are vector 
rather than scalar. 

To determine the effective potential drop, 
effective current, and power delivered between 
any two points of the circuit in Fig. 1, the circuit 
is first represented as made up of impedances. 
Impedance Z=V/JI, where V is the effective 
potential difference between two points and J is 
the effective current between those two points. 
Impedances are inductive, resistive, capacitive, 
and combinations of these. The relation between 
inductive impedance and inductance is Z=2zfL, 
where L is the inductance and f is the frequency ; 
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that between resistive impedance and resistance 
is Z=R; and that between capacitive impedance 
and capacitance is Z=1/2zfC, where C is the 
capacitance. 

The laws which describe alternating current 
circuits may be considered as laws discovered by 
experiment. Suppose the potential drops in the 
series circuit of Fig. 2 are measured with a 
voltmeter which reads effective potential differ- 
ences. It is discovered that with switch K closed 
Vap¥ Vas+Vact+ Ven where Van, Vaz, Vac, 
Vep represent voltmeter readings. When, how- 
ever, the voltages are added as vectors, as shown 
in Fig. 3, it is found that Vap= Vas+ Vact Veo. 
In taking the vector sum of potential drops, if the 
potential drop across inductance is considered as 
a vector vertically upward, the potential drop 
across resistance is a vector horizontally to the 
right, and the potential drop across capacitance 
is a vector vertically downward. Further meas- 
urements of electromotive force between D and 
G, G and A and D and A, Fig. 2, made with 
switch K open so that no current flows through 
the generators, shows that Exap#Eact+Eep. If, 
however, the electromotive forces are considered 
as vectors such that the angle between Eg and 
Egp is the phase difference (¢2.—¢1) between the 
two generators it is found that Eap=Eac+Eep 
or E= E,+ £2, as shown in Fig. 4. Any convenient 
direction is chosen for E;, and £2 is drawn in a 
direction making an angle (¢2:—¢1) with EF. 

These measurements and others like them are 
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Fic. 1. Alternating current circuit. 
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described by Kirchhoff’s first law for alternating 
current circuits: (1) In any closed path the vector 
sum of the potential drops equals the vector sum of 
the electromotive forces. 

With an ammeter which reads effective cur- 
rent, measurements of the total current through 
a group of impedances in parallel and the separate 
currents through each of the individual impe- 
dances, as for example in the circuit of Fig. 1, 
show that J#J,+J2.+J3. If, however, current is 
considered as a vector, then J=J,+J2.+TJ; as 
shown in Fig. 5. In making the vector sum, if 
current through inductance is considered as a 
vector upwards, current through resistance is a 
vector horizontally to the right, and current 
through capacitance is a vector vertically down. 
Kirchhoff’s second law for alternating current 
circuits may therefore be stated : (2) At any point 
in a network the vector sum of the currents ap- 
proaching that point equals the vector sum of the 
currents leaving that point. 

For impedances in series, from Kirchhoff’s first 
law for alternating current circuits one may write, 
for any two points in the circuit, V=Vi+V2 
+V3;+--- where V is the potential difference 
between the two points and Vi, V2, V3--- are the 
separate potential differences for impedances or 
groups of impedances in series between the two 
points. Replacing each potential difference by ZI 
and accepting as an experimental observation the 
fact that the current is the same for all impedances 
in series, one obtains Z7=Z,I+Z.J+2Z3I. Can- 
celling J leaves Z=Z1:+Z2+2Z;3;+---. This is the 
same as the equation for resistances in series. For 
impedances, however, the equation represents a 
vector addition. If inductive impedance is con- 
sidered as a vector in a direction vertically up, it 
is horizontally to the right if resistive, and 
vertically down if capacitive. The resultant impe- 
dance may be obtained by graphical addition. 

For impedances in parallel, the second of 
Kirchhoff’s laws for alternating current circuits 
may be written J=J,+J2+Js, where J;, I2, and 


Vec 


Veo Fic. 3. Vector addition 
of voltages for series a.c. 
circuit. 


and 


I; are vectors to be added graphically. Experi- 
ment shows that the potential difference across 
impedances in parallel is the same for all impe- 
dances. Since J=V/Z, then (V/Z)=(V/2Z;) 
+(V/Z2)+(V/Zs). Cancelling V leaves (1/Z) 
=(1/Z1)+(1/Z2)+(1/Z3) where again one finds 
that the equation for impedances in parallel is the 
same as that for resistances in parallel except that 
for impedances in parallel the addition must be 
made by the vector method. In making this 
addition if 1/Z is considered as a vector vertically 
up for an inductive impedance, 1/Z is horizontally 
to the right for a resistive impedance and verti- 
cally down for a capacitive impedance. 


end 
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Fic. 4. Vector addition of electromotive forces 
for series a.c. circuit. 

The circuit of Fig. 1 can now be solved by the 
method of simplification of circuit following the 
same procedure as that followed for direct current 
circuits except for the substitution of vector 
addition for scalar addition. Each group of impe- 
dances in parallel is first reduced to a single 
impedance using the equation (1/Z»)=(1/Z1) 
+(1/Z:2)+(1/Z3). The vector addition of these 
reciprocals is shown in Fig. 6 and the resultant 
vector, 1/Z» shows the direction to be ascribed to 
Zo, the impedance of the group. Repeating the 
procedure for the second group in parallel, the 
vector addition is shown in Fig. 7. The resultant 
impedance Zo and its direction are thus obtained. 
The circuit has now been reduced to one with a 
number of impedances all in series, as shown in 
Fig. 8. With application of the equation for im- 
pedances in series Z=Z,1+Z2+---, one obtains 
the resultant impedance from the vector addition 
shown in Fig. 9. The internal impedances of the 
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generators Z; and Zy are considered inductive in 
this problem. The resultant impedance Z and its 
direction are determined from measurements on 
the diagram. The two alternating current genera- 
tors are replaced by a single alternating current 
generator by the vector addition of the electro- 
motive forces, as shown in Fig. 4. The circuit is 
thus simplified to a single a.c. generator and 
single impedance as in Fig. 10. For this simplified 
circuit the current J = E/Z. This effective current 
for the simplified circuit is also the current for the 
series circuit, Fig. 8. From the current and indi- 
vidual impedances in the series circuit, all of the 
potential drops can be determined with the equa- 
tion V=ZI. These potential drops are recognized 
to be the same as the potential drops between the 
corresponding points of the original circuit, Fig. 1. 
Since now all potential drops are known for the 
original circuit, the individual currents between 
any two points can be determined by the equa- 
tion T= V/Z. 

As a check on the accuracy of the results, the 
individual currents for several impedances in 


end 


Fic. 5. Vector addition start 
of currents for a group of 
impedances in parallel. l 


parallel may be added vectorially as in Fig. 5. 
The vector sum should equal in magnitude the 
current calculated for the simplified circuit. 

The power delivered by the generators is 
P=EIk, where k is the power factor. The power 
factor for this circuit is Zre/Z where Ze is the 
horizontal component of the resultant impedance 
Z, as shown in Fig. 9. 

Power expended is given by P=EIZ,/Z. 
From Ohm’s law, however, E/Z=I where Z is 
the resultant impedance for the series circuit. 
Thus P=ZprlI*. Here Zp is the horizontal com- 
ponent of the entire impedance in the series 
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Fic. 7. Vector addition 
of reciprocal impedances 
for parallel group consist- 
ing of Z, and Zs. 


circuit, J is the current and P the power delivered 
by the generators. 

The power delivered between any two points 
such as A and B in the series circuit, Fig. 8, is also 
P=VIXpower factor, where V is the voltage 
between the two points. The power factor A to B 
is (Zo)r/Zo, in which (Zo)p is the horizontal com- 
ponent of Zo, as shown in Fig. 9. 

Thus 


| (Zor 
Pas= abl . 
Z 


0 
But by Ohms law Vaz/Zo=TI, so Pas=(Zo)r/’. 

In general, the power delivered between any 
two points in a series circuit is P=Z,rl*, where 
Zr is the horizontal component of the impedance 
between the two points. 

The power delivered to any one impedance Z 
in a group of impedances connected in parallel is 
P=VIk, where in this case the power factor is 
given by (1/Z)r/(1/Z), where (1/Z)r is the 
horizontal component of 1/Z, V is the potential 
difference across Z, and J is the current 
through Z. By Ohm’s law J=V/Z. Thus 
P=(V?/Z)[(1/Z)r/1/Z]. The quantity 1/Z can- 
cels and P= V*(1/Z)r, where P is the power de- 
livered to any one impedance Z of a group in 
parallel and V is the potential drop across the 
group. 

For either series or parallel connections the 
power is zero for pure inductive or capacitive 
impedances. ‘ 

When all currents, potential drops, and powers 
are known, the a.c. circuit may be considered to 
be as completely solved as the equivalent d.c. 
circuit. For a.c. circuits, however, it is important 
to determine also the phase relations. 


1 
Fic. 8. Reduced circuit 


equivalent to circuit of 
Fig. 1. 
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Fic. 9. Vector addition of impedances for series a.c. circuit 
of Fig. 8. 


The phase relation between current and voltage 
for any part of a series circuit such as that in 
Fig. 8 is shown in Fig. 9. The current for any part 
of the series circuit is represented by a vector in 
the direction of Zz; and the voltage between any 
two points, as a vector in the direction of the 
impedance between the two points. For the cir- 
cuit as a whole the phase difference between 
current and voltage equals the angle between Zr 
and Z. For the impedance Zo, the phase differ- 
ence between voltage and current is the angle 
between Zr and Zo. 

To show lead and lag, vectors in Fig. 9 for 
parts of a series circuit may be considered as all 
rotating counterclockwise. For the given circuit 
as a whole the voltage is ahead of the current. 
For Zo, the voltage is behind the current. 

The phase relation between voltage and current 
for any one impedance Z in a group of impedances 
in parallel is shown by Fig. 6. The voltage for all 
impedances in the group is represented by a 
horizontal line parallel to (1/Zo)z and the current 
through any impedance is considered to have the 
same direction as 1/Z. Here Zo is the combined 
impedance of the group and (1/Zo)z is the hori- 
zontal component of 1/Z». The angle between the 
vector 1/Z and (1/Zo)r is the phase difference 
between the current and voltage for the impe- 
dance Z. 

The lead and lag are shown if, for impedances 
in parallel, the voltage, current, and 1/Z are 
considered to be vectors all rotating clockwise. 
Thus for the group of combined impedance Zo, 
the voltage leads the current, but in the capaci- 
tive impedance Zi, the current leads the voltage. 


In the use of vector additions to determine resultant 
impedances for several impedances in series, or in parallel, 
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it is recognized that impedance is not an ordinary vector. 
Voltage and current are ordinary vectors but sometimes 
called phasors, and impedance is a multiplier such that 
V=ZI. If I is not parallel to V, the impedance is such a 
multiplier that ZI is a vector parallel to V. When several 
impedances are connected in series, the resultant impedance 
is obtained by the vector addition of the separate impe- 
dances. What are really being added are the potential 
differences, not the impedances. The equation Z=Zi+Z:2 
+--+ was obtained from ZJ=Z,I+Z2I+---+ by cancelling 
I. The current has the same magnitude and direction for 
all, but 2:7, Z2I, and Z;3I are vectors in directions different 
from J. The magnitude of the current is cancelled. There 
remains a unit current vector multiplier for each impe- 
dance, so that each impedance times its unit current vector 
is a vector in the direction of the voltage across the impe- 
dance. It is convenient to take this into account by treating 


Fic. 10. Simplified series a.c. 
circuit equivalent to circuit of 
Fig. 1. 


impedance as though it were a vector and thus obtain the 
resultant impedance for the group in series. 

The same situation occurs when impedances in parallel 
are combined. The vector addition of (1/Z)=(1/Z:) 
+(1/Z:2) is really an addition of currents (V/Z)=(V/Z1) 
+(V/Z:2)--- with Z as a divider for the potential differ- 
ences. The potential differences all have the same magni- 
tude and direction. The impedance is of such a character 
that V/Z:1, V/Ze, «++ are vectors in the directions of the 
currents. The magnitudes of the potential differences are 
cancelled. The directions, however, still persist so that 
(1/Z) = (1/21) +(1/Z:2) represents unit potential difference 
vectors divided by Z, Z:, or Z2. This is taken into account 
by considering that 1/Z for an impedance is a vector in the 
same direction as the current through that impedance. 


In summary, the solution of an a.c. circuit 
problem like that of Fig. 1 consists of representing 
impedance, and the reciprocal of impedance, as 
vector quantities having directions vertically up 
if pure inductive, horizontally to the right if pure 
resistive, and vertically down if pure capacitive. 
With this convention, the procedure of simplifica- 
tion of circuit to a single electromotive force and 
impedance, followed by expansion back to the 
original circuit can be followed for both a.c. 
and d.c. circuits. This procedure enables a stu- 
dent of elementary physics to learn how to solve 
a.c. circuit problems with a minimum of diff- 
culty. Those who go on to courses in electrical 
engineering will likely have little difficulty in 
learning a new convention for directions. 

The authors wish to express thanks to Pro- 
fessor R. R. Benedict and Dr. Stanley C. Snow- 
don for helpful discussions. 





The Laplace Transform as a Form of Curve Fitting 


Horace M. TRENT 
Naval Research Laboratory, Washington 20, D. C. 


HE importance of the Laplace transforma- 

tion technique for the analysis of linear 
systems needs no justification in. the minds of 
physicists and engineers. The ever increasing 
volume of literature appearing in scientific jour- 
nals and employing these methods is sufficient 
motivation for teaching the fundamentals of the 
subject to physics students at an intermediate 
level. 

The conventional presentation of the subject 
begins with the formal definition of the Laplace 
transform of a function of time; that is, with the 
formula, 


(1) 


Lf(t) =f f(t)e-*dt = F(s). 
0 


It is pointed out that this transformation, in 
general, converts an integro-differential equation 
in real variables to an algebraic equation in 
complex variables. The resulting simplification 
in mathematical manipulation is obvious to all. 
Now the more sophisticated students are able to 
accept this presentation without reservation. The 
majority, however, find great difficulty in de- 
veloping an early confidence in the method since 
they see no reason why this definition was 
adopted other than, ‘it works.’’ In short, they 
have no physical insight into the process and 
have no mental picture to accompany their 
manipulations. The situation becomes even more 
confused after the introduction of the inverse 
transformation formula, 


1 ct+jo 
f(t) -—j F(s)e**ds. (2) 
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This relation is presented supposedly to complete 
the theory of the process but is seldom used in 
physical problems. The unusual limits are con- 
fusing, to say the least. 

The author has found that the basic Laplace 
formulas can be developed in a few easy stages 
from the well-known process of curve fitting and 
that a knowledge of this development leads to 
clearer interpretations of the Laplace techniques 


when applied to physical problems. It seems 
appropriate, therefore, to present such a treat- 
ment of the Laplace transformations to the 
readers of the Journal. The presentation will 
begin with a brief review of the curve fitting 
process. 


Elements of Curve Fitting 


Suppose there is given a function of time such 
as curve A of Fig. 1. It is well known that an 
assumed arbitrary function can be made to pass 
through specified points of the given curve 
provided there are constants in the assumed 
function that can be adjusted. For example, the 
function f(t)(=)ado+ait, in which the symbol 
(=) is used to indicate equality at a finite 
number of points, can be made to pass through 
the points f(t:) and f(t,) as indicated by the 
curve B since there are two adjustable constants, 
@ and a;. On the other hand, the function 
f(2)(=)ao+ait+aet?+a;3t? can be made to pass 
through the four points, f(t:), f(te), f(ts), and 
f(ts) as shown in the curve C. The improvement 
in the accuracy of the fitting process comes from 
an increase in the number of points of inter- 
section. If this idea is extended, an exact fit can 
be obtained provided the assumed function 
agrees with the given curve at all points. This 
implies that there must be an infinite number of 
constants which can be adjusted. 

It is important to note that the form of the 
assumed function is not unique. For example, 
the curve of Fig. 1 could be fitted at the same 
four points by the function f(t)(=)ao+ait+a2l* 


Fic. 1. Curves approximating a given function. 
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+a; sin2¢. Hence, the choice of the fitting func- 
tion must be dictated by other considerations, 
chief of which is ease in evaluating the constants 
in the assumed expression. Such a choice will 
now be made which is a series expression and 
which includes the Fourier series as a special case. 


A General Expression for Fitting a Curve 


A simple, yet versatile form for fitting a curve 
exactly in the interval 0 to T is the series ex- 
pression, 


2 [ao 2r 2an 
f(t) =e ase COS—+#: - -,e°° Cos——t: : 
TL2 ie 


2a 2rn 
+by,e* sin—t- - -b,e°* sin, (3) 
Fr r 


where ¢ is an arbitrary constant. This series can 
be written more compactly as 


fo) =— [f+z ane 
t)=—e*"| — a, Cos——t 
T \2 T 


n=1 


w 2rn 
+> bd, sin] (4) 


n=l 


It will now be shown that the constants, a, 
and b,, can be computed from simple formulas 
and that the basic Laplace transformations (Eqs. 
(1) and (2)) follow from this expression for curve 
fitting. 

If Eq. (4) is multiplied by e~°‘dt and integrated 
from 0 to T, it follows that, 


= 
ao= f f(t)e-etdt. (5) 


Likewise, by multiplying by e~* cos(2rn/T)tdt 
and e~* sin(2an/T)idt in turn and integrating, 
the two relations, 


‘ 7 
a,= f f(t)e-* cos(2xn/T)tdt, (6) 
0 


T 
n= f f(t}e-* sin(2xn/T)tdt, (7) 


are obtained. 
These last three equations can be made more 
compact. Let 21rm/T =w and let j represent the 
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imaginary number (—1)#. An easy calculation 
will show that 


T 
an—jba= f f(he-*te-'dt = F(c, nwo), (8) 
0 


of which ¢ is an assumed constant. This is the 
basic formula for computing the constants of 
Eq. (4) and is true for any value of including 
zero. A similar calculation shows that 


7 
ant jba= f f(the-teimotdt = Fc, nwo). (9) 
0 


It is to be noted that F(c, nwo) is the complex 
conjugate of F(c, nwo) and that 
P(c, nwo) = F(c, —nwo), 


Qn =4([F(c, nwo) + F(c, nao) ], 


(10) 
(11) 


1 
bn =—[F(c, nwo) — F(c, nwo) ]. (12) 
2j 


If the basic curve fitting relation (4) is recon- 
structed in the light of Eqs. (11) and (12), it is 
found that , 


f(t) =—e'L F(c, 0)+ > F(c, nwo)eirot 
2r n=1 


+ > P(c, nwo)e™°"], (13) 


n=1 


By Eq. (9), this becomes 


f(t) = Sel, 0)+ > F(e, Nw) eirwot 


n=l 


+ re F(c, —mwo)e— i" ], 


n=l 


Let 1 be replaced by —x in the last summation, 
yielding 


fle) =e L Fle, 0)+ Fle, maxon 


n=l 


+ x F(c, nmwo)e?™°* J. 


n=—1 


It is obvious that the right-hand side of the last 
equation can be represented by the single sum- 
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mation 


w oo 
JQ) =e DX Fle, meo)eineot, 


wv n=—D 


(14) 


Let us recapitulate the steps that have been 
taken. We have adopted the expression 


2 [a -@ 2mm = x _ aan 
—e*| —+ ¥ a, cos—t+ > b, sin——t 
T 2 n=1 n=1 rT 
as particularly suitable for fitting exactly a func- 
tion of time between the limits =0 and t=T. 
The expression involves an arbitrary constant c, 
the need for which will appear later, and a set 
of constants, a, and b,. It was shown that any 


pair of constants a, and b, can be calculated from 
the relation 


T 
a,—jda= f f(t)e—*te— 0'dt = F(c, nwo). (15) 
0 


It was also shown that if the constants are 
known, then the expression can be constructed 
by means of the relation 


wo 90 ; 
f(t)=—e* DS Fic, nwo)e*, 


T  n=—% 


(16) 


These relations (15) and (16) can be viewed as 
direct and inverse transformations, the first 
generating a set of coefficients when a function 
is known, the second generating a function when 
the coefficients are given. Doubtless, the mathe- 
matically inclined will be bothered at this stage 
with existence considerations. These have been 
passed over roughshod in an effort to establish a 
simple sequence of ideas. It should be obvious, 
however, that theorems for the existence of 
Fourier series apply now to the function e~“f(é) 
rather than f(t). This leads to greater freedom as 
will be mentioned later. 


Extension of the Range of Representation 


So far, the given function f(#) has been fitted 
only over the finite time interval 0 to T. Suppose 
the interval is enlarged until T approaches 
infinity? Let us assume this can be done. It is to 
be observed that wo=(22/T) approaches zero. 
Now wo is the common difference in angular 
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frequency of the harmonic factors between con- 
secutive values of m. Hence, in the limit, F(c, nwo) 
or a, and 6b, exist at all frequencies. Let us 
simplify our expression by writing w for nwo. 
The expression (15) becomes in the limit 


F(e, =f f(the-“e- “dt. (17) 
0 

Two remarks are in order. In the first place, even 
though a, and b,, as defined by Eq. (4), may be 
finite, the coefficient of each term in the basic 
series is an infinitesimal since each term is 
multiplied by 2/T=wo/z. In the second place, 
Eq. (17) is just like a Fourier transformation 
except that f(t)e-*t replaces the usual f(t) as 
found in the Fourier integral. Hence our processes 
are valid for the class of functions f(#)e~*t for 
which a Fourier integral representation is pos- 
sible; hence any existence theorem on such 
integrals can be applied directly to our situation. 
Let it be recalled that ¢ is an arbitrary constant. 
It is just this freedom in the choice of ¢ which 
allows us to make some functions well-behaved 
which would not be otherwise. For example, the 
common excitation function, 


fi) =0, t<0 


=Asinpt, t>0, 


has no Fourier representation but can be repre- 
sented by our relation (4) provided c>0. 

The role played by the arbitrary constant c 
can be sensed in a general way by examining our 
curve fitting process from another point of view. 
Equation (4) shows that f(#) is fitted exactly by 
an infinite set of sinusoids which may decay, 
expand, or be of constant amplitude depending 
on the value of c. Consider fitting the function, 
f(j)=Ae-*'. Equation (17) can exist only if 
c>—pP, that is, if the individual terms in Eq. (4) 
decay more slowly that the original function. 
Again, if f(t) = Be, then c>k; i.e., the terms of 
Eq. (4) expand at a faster rate than e*'. As a 
third example, consider the common step-func- 


tion 
f() =0, 
f®=1, 
This function has no Fourier representation but 
can be fitted by Eq. (4) provided c>0. Since the 


t<0 
t>0. 
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original function neither decays nor expands, it 
can be fitted provided our terms expand at any 
finite rate. Similar considerations will show that 
our basic relation cannot be made to fit the 
function f(t) = D(e”), since no finite value of c can 
be found which will make the terms in our rela- 
tion expand at a faster rate than the original 
function. Curve fitting with harmonic terms of 
constant amplitude is well known from the 
Fourier series but fitting by use of expanding 
terms is not so well known even though they 
are more versatile. 
If T approaches infinity, Eq. (16) becomes 


1 «@ 
aie d eF(c, w)e#*two. 


T w= 


But wo becomes infinitesimally small in the limit; 
hence the sum can be represented by the integral, 


1 2 
fy=— f e*' F(c, we dw. 


2 ¥ _« 


(18) 


Equations (17) and (18) can be called direct and 
inverse transformation formulas in the same 
fashion as for Eqs. (15) and (16); useful when 
the range of representation extends from t=0 
tot=o, 


The Laplace Transformation Formulas 


An examination of the relations (14)—(18) dis- 
closes that c and w always occur in the combina- 
tion c+jw. Therefore, let a new variable s be 
defined by the relation 


s=c+ jw. (19) 
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It follows that ds = jdw, s=c—jo whenw=—o, 


and s=c+jo when w=+. Using these rela- 
tions, Eqs. (17) and (18) become 


(20) 


F(s) -f f(t)e~*“dt, 


1 ct+ jx 
f(t) -— f F(s)e*ds, (21) 


TJ] “c—jn 


which are the direct and inverse Laplace trans- 
formation formulas. 


Summary 


It is well known that there are numerous sets 
of functions which can be employed to approxi- 
mate a more or less arbitrary function to any 
desired accuracy over a finite real interval. If 
then a set is chosen, each term of which is an 
exponentially decaying or expanding sinusoid, 
the sum of such a set can be made to fit a broad 
class of physical functions. Formulas for com- 
puting the coefficients of the terms (direct trans- 
formation) and for generating the function from 
a set of the coefficients (inverse transformation) 
can be established. These formulas are simple 
and compact. When the range of representation 
is extended from t=0 to t= ~, the formulas still 
remain simple and compact and are valid for 
most physical functions by virtue of an arbitrary 
constant which determines the rate of decay or 
growth of the individual terms in the series 
representation. The Laplace transformation for- 
mulas follow naturally. 


There is another point which is important with regard to this country’s position in particular. 
If we adopt policies that are too restrictive about scientific knowledge, other countries will also do the 
same. As CHARLES KETTERING of the General Motors Corporation has put it, ‘‘When you lock 
the doors of the laboratory, you lock out more than you lock in.” This result could be especially 
weakening to America for we have not thus far been very productive in truly fundamental science. 
On the other hand we are vevy well equipped to make rapid progress in exploiting technological uses 
of such fundamental knowledge. For these reasons we have to guard against policies that are so 
narrowly restrictive that we stifle our own development.—E. U. Conpon (1948). 





The Action Principles 


H. M. Dapourran 
Trinity College, Hartford, Connecticut 


TUDENTS of the foundations of mechanics, 
particularly Heinrich Hertz,! Ernst Mach,? 
and Karl Pearson,’ have criticized Newton’s laws 
of motion and have pointed out their limitations. 
That these laws do not form a sufficiently satis- 
factory foundation of mechanics is brought out 
by the fact that, for about sixty years after 
the publication of Newton’s Principia, eminent 
physicists and mathematicians had to have re- 
course to other laws in order to solve many 
problems of mechanics. 


The postulate of Huyghens, “that if any 
weights are put in motion by the force of 
gravity they cannot move so that the center 
of gravity of them all shall rise higher than 
the place from which it descended,” was 
generally one of the principles of the solu- 
tion: but other principles were always 
needed in addition to this, and it required 
the exercise of ingenuity and skill to detect 
the most suitable in each case. Such prob- 
lems were for some time a sort of trial of 
strength among mathematicians. The Traité 
de Dynamique published by D’Alembert in 
1743 put an end to this kind of challenge by 
supplying a direct and general method of 
resolving, or at least throwing into equa- 
tions, any imaginable problem. The me- 
chanical difficulties were in this way reduced 
to difficulties of pure mathematics. 


The limitations of Newton’s laws are due to 
the facts that the first law is a special case of 
the second; the second law applies directly to a 
body which may be considered as a particle; and 
the third law applies to two interacting particles. 
It must be stated that laws applicable to more 
complex systems can be derived from Newton’s 
laws ; but the presentation of mechanics becomes 
simpler when the subject is developed from a 


1 Heinrich Hertz, Principles of mechanics (Macmillan). 
a Mach, Science of mechanics (Open Court, ed. 5, 
1 ; 

3 Karl Pearson, Grammar of science (A. & C. Black, 
London, ed. 2, 1900). 

*Edward J. Routh, Elementary rigid dynamics (Mac- 
millan, ed. 7, 1905). j 


broad general principle from which laws appli- 
cable to special cases are derived, than when this 
process is reversed. For this reason, the presenta- 
tion of advanced mechanics is far more coherent 
than that of elementary mechanics. 

Since the publication of D’Alembert’s prin- 
ciple, a number of other general principles have 
been formulated and made the basis of advanced 
texts on mechanics. Unfortunately, none of these 
principles is suited to elementary mechanics. 
Even the simplest of them, D’Alembert’s prin- 
ciple, is too formidable to be placed in the hands 
of the beginner in its original form. This principle 
can, however, be stated and is stated, in my 
Analytical Mechanics,® in a form which is suited 
to elementary texts. 

The new form consists of twin principles which 
I have called the linear action principle and the 
angular action principle. The linear action prin- 
ciple states : 


The vector sum of all the external linear 
actions to which a system of particles or 
any part of it is subject at any instant 
equals zero. }A=0. 


The angular action principle, which is derived 
from the linear action principle, states: 


The vector sum of all the external angular 
actions to which a system of particles or 
any part of it is subject at any instant 
equals zero. }-Aa=0. 


Linear actions are divided into two classes, 
namely, forces and kinetic reactions. A force is a 
vector which represents the pull or push of one 
particle on another. A kinetic reaction is a 
vector which represents the action of the ma- 
terial universe (the inertial field) on a particle. 
It acts only while the particle is accelerated. It is 
proportional to and opposed to the acceleration. 
The constant of proportionality is the mass of 
the particle. With these definitions the linear 


5H. M. Dadourian, Analytical mechanics, (Van Nostrand, 
ed. 3, 1931). 
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Fic. 1. Motion of 
a water-filled spherical 
shell on an_ inclined 
plane. Rotation is con- 
sidered about the point 


Hn. 


Fic. 2. Motion of 
a water-filled spherical 
shell on an_ inclined 
plane. Rotation is con- 
sidered about the cen- 


of contact between 


ter of the shell. 
sphere and plane. 


action principle may be written in the form, 


> (F—mv) =0, (1) 


where F is the force acting on a particle of mass m 
that possesses a velocity v. 

Angular actions are divided into two classes, 
namely, torques and angular kinetic reactions. 
A torque is the moment of a force; and an 
angular kinetic reaction is the moment of a 
(linear) kinetic reaction. Hence the angular ac- 


tion principle may be written in the form 
X(G—mvp) =0, (2) 


where G denotes torques and p the lever arm 
of mv. (The vector product sign is omitted in 
Eq. (2) for pedagogical reasons and also because 
rotational motions considered in elementary me- 
chanics are uniplanar motions, in which case the 
vector product notation need not be introduced.) 

When resolved into component equations, 
Eq. (1) yields the first three, and Eq. (2) the 
last three of the six differential equations of 
D’Alembert’s principle. 

If the system consists of a single particle, or 
many particles, each of which has the same 
acceleration, Eq. (1) reduces to 


> F-—mv=0, ; (1’) 


which is Newton’s second law of motion. 
If the system consists of a rigid body in uni- 
planar motion, Eq. (2) reduces to 


rG—Ie=0, (2') 


when J remains constant during the motion; 


DADOURIAN 


and to 
> G —d/dt(Iw) =0, (2’’) 


when J varies. 

The conditions for the equilibrium of a par- 
ticle, Newton’s third law, the law on the sum 
of internal forces, the law of the motion of the 
center of mass, the principle of the conservation 
of linear momentum, etc., can be derived from 
Eq. (1). Analogous conditions, laws and prin- 
ciples which relate to rotation can be derived 
from Eq. (2). 


Advantages of the Action Principles 


The action principles, in spite of their broad 
generality, can be understood by the student, to 
the extent of his needs, at any stage of his 
study of mechanics. The principles lend them- 
selves to a gradual unfolding of their implications 
and, consequently, to the presentation of me- 
chanics in a logically coherent and gently graded 
way. The almost identical wording of the two 
principles brings out clearly the fact that there is 
a perfect analogy between the magnitudes in- 
volved in motion of translation and motion of 
rotation. 

The concept of kinetic reaction affords the 
operational definition of kinetic energy as the 
result of the work done by a force against the 
kinetic reaction, as potential energy is defined 
as the result of the work done against a con- 
servative force. 

The wording of each principle brings out 
clearly certain points of great pedagogical im- 
portance. The action principle makes clear, for 
instance, the fact that internal forces should not 
be considered, and the fact that forces which act 
on (not those exerted by) the body whose equi- 
librium or motion is discussed, should be con- 
sidered. 

The clear-cut distinction between force and 
kinetic reaction eliminates the contradictions and 
inconsistencies which result from the use of the 
terms “‘fictitious force,” ‘imaginary force,” “‘cen- 
trifugal force,’ and ‘‘reversed effective force.”’ 
The contradiction between the condition of 
equilibrium, according to which there is no 
acceleration when the sum of the forces equals 
zero, and the condition of motion, in which there 
is acceleration although the sum of the “forces’’ 
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equals zero is removed by the application of the 
action principles. The inconsistency of calling 
the normal component of the kinetic reaction 
centrifugal force without giving the tangential 
component a special name, say, afugal force 
(away-fleeing force) and the inconsistency of 
designating by the word force both the cause 
and the effect of acceleration cannot be avoided 
by the device of attaching an adjective to the 
word. 


Can we without confusing our ideas, 
suddenly begin to speak of forces which arise 
through motion, which are a consequence of 
motion? Can we behave as if we had already 
asserted anything about forces of this new 
kind in our laws, as if by calling them 
forces we could invest them with the proper- 
ties of forces? These questions must be 
clearly answered in the negative. The only 
possible explanation is that, properly speak- 
ing, centrifugal force is not a force at all. 


Problems, such as the following, that would 
appear difficult to a student who has learned 
only Newton’s laws, should be relatively simple 
to one who is trained in the application of the 
action principles. 

Discuss the motion of a thin spherical shell, 
which rolls down an inclined plane; the shell is 
full of water and its inner surface is perfectly 
smooth. 

Solution: (a) Taking the axis of rotation 
through the point of contact, and applying the 
angular action principle to Fig. 1, we get 


(m+m’')gr sina—mir—Ipo=0, 


PRINCIPLES 


or 
(m+m’') gr sina — mir —5m'ri/3 =0, 


where m is the mass of the water and m’ the 
mass of the shell, 7 is the radius of the shell, 
a is the angle of the plane to the horizontal, 
® is the linear acceleration of the center of 
mass, » is the angular acceleration of the 
spherical shell, and Jo is the moment of inertia 
of the shell relative to the axis through point O, 
perpendicular to the plane of motion. Hence 


_ _(m+m'’) sina 
= $——_-_9 


v= 
3m+5m’ 


Since the acceleration is constant, we have 
v=Votit, and s=vot+ ir. 

(b) Taking the axis of rotation through the 
center of the shell, and applying both action 
principles to Fig. 2, we have 


(m+m’)g sina— F—(m+m’)i=0, 
N—(m+m’)g cosa=0, 


Fr—I,@=0, 


where J, is the moment of inertia relative to 
the axis through the center of the shell. Solving 
these equations for +, we obtain the same result 
as that in Part (a). 

Discussion: If m=0, that is, if the shell is 
empty, ’=2g sina, as it should. If m’=0, that 
is, if the mass of the shell is negligible, i=g sina; 
hence the motion is the same as that of a body 
which slides down a smooth inclined plane. 


He is still at play, save only that his play is such as manhood stops to watch, and his playthings 
are those which the gods gave their children. The universe is his box of toys. He dabbles his fingers 
in the day-fall. He is gold-dusty with trembling amidst the stars. He makes bright mischief with 
the moon. The meteors nuzzle their noses in his hand. He teases into growling the kennelled thunder, 
and laughs at the shaking of its fiery chain. He dances in and out of the gates of heaven; its floor is 
littered with his broken fancies. He runs wild over the fields of ether. He chases the rolling world. 
He gets between the feet of the horses of the sun. He stands in the lap of patient Nature, and twines 
her loosened tresses after a hundred wilful fashions, to see how she will look nicest in his song.— 
FRANCIS THOMPSON, describing SHELLEY. Quoted by A. S. Eve as applicable to LorD RUTHER- 


FORD. 
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A Simple Fresnel Diffraction Experiment 


L. A. SANDERMAN AND R. S. BRADFORD 
University of Washington, Seattle, Washington 


HE following laboratory experiment involving Fresnel 
diffraction has been devised to enable our students 

to study the phenomena of diffraction in the elementary 
laboratory and to measure the wavelength of sodium light 
with a fair degree of accuracy. The necessary equipment is 
available in any elementary optics laboratory and the 
experiment may be performed in a partially darkened room. 
The essential parts of the apparatus are shown in Fig. 1. 
The point source of light L consists of a sodium lamp with 
a pinhole screen placed in front of it. The circular aperture 
D, greatly magnified in Fig. 1, is a 2-mm hole drilled in a 
sheet of thin metal and mounted on an optical bench. The 
observation of the diffraction pattern is accomplished by a 
simple lens or instrument eyepiece P mounted on the 
optical bench behind the aperture in such a way that the 
eyepiece may be moved along the bench. A set of cross 
hairs or an indicating pointer is mounted at the focal plane 
of the eyepiece lens to provide a reference point for the 
center of the diffraction pattern. As the eyepiece is moved 
away from the circular opening, fewer Fresnel zones of the 
circular aperture D contribute to the illumination of the 
center of the diffraction pattern. Therefore, the center of the 


Fic. 1. Diagram showing the relationship between the wavelength of the 
light used and the constants of the apparatus. 


pattern becomes alternately dark and bright as an even or 
odd number of zones contribute to the illumination. 

In performing the experiment the student places the 
eyepiece near to the circular aperture D and slowly moves 
it away, noting the distances to the focal plane of the 
eyepiece from the circular aperture at the consecutive 
positions of maximum brightness and darkness of the center 
of the observed pattern. For a source distance of 2 meters 
and an optical bench of the same length, approximately 15 
positions of light and dark centers can readily be identified. 

The following theory will show how the data obtained 
may be used to determine the wavelength of the light used. 
As shown in Fig. 1, let a be the distance from observed 
pattern to circular opening, b be the distance from point 
source to circular aperture, 7 be the radius of circular 
aperture, 2 be the number of Fresnel zones contributing to 
a given diffraction pattern, x be the sagitta of the wave 
front arc ACB, and ) be the wavelength of light. 

In the right triangle AEP 


a+r=(a—x+n/2)? (1) 


r? = —2ax+anrd\+x?—xnd+n7r?2/4. (2) 


Since a>x and a>, we may neglect the last three terms 
of Eq. (2) and obtain 


r?= —2ax+any. (3) 


The approximate sagitta equation for the condition when 
b>x is 


x=r?/2b. (4) 
Eliminating x between Eqs. (3) and (4) we obtain 
1/b+1/a=ny/r*. (5) 


It is interesting to note that Eq. (5) has the same form as 
the usual lens or mirror equation. Equation (5) may be 
written in the following form 


n[ab/(a+b) J=r?/n. (6) 


Since the right hand member of Eq. (6) is a constant, the 
quantity n[ab/(a+b)] is also a constant for the corres- 
ponding values of ” and a. 

There are several conditions imposed upon the series of 
whole numbers 2 which will help the student to correlate 
the different values of a with the appropriate 7. For a light- 
centered diffraction pattern m is an even number and an odd 
number for a dark-centered pattern. Also, the values of 
are consecutive whole numbers for consecutive light and 
dark centers, the smallest value of » corresponding to the 
greatest distance CP. 


TABLE I. Typical data.* 


ab ab 
Corresponding 1 
n a+b 


150.5 
150.8 
152.6 
154.3 
149.9 
152.0 
150.3 
151.2 
152.7 
150.7 
148.2 
149.6 
154.2 


Mean 151.3 


* b =232.0 cm; r =0.0947 cm. Calculated wavelength 5920A. Wave- 
length of light used 5893A. 


Thus by correlating each value of a with the appropriate 
n to obtain the best agreement with Eq. (6), a mean value 
of r?/X may be obtained. The diameter of the aperture 
having been determined by means of a measuring micro- 
scope, a value of the wavelength of the light may be ob- 
tained by the use of Eq. (6) and compared with the ac- 
cepted value. 

The data in Table I were chosen at random from student 
reports. 
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A Self-Sustaining Electronichord 


W. Davip BEMMELS 
Ottawa University, Ottawa, Kansas 


HE name “electronichord” was used by Urquhart! in 

describing a demonstration device for making audible 
the different modes of vibration of a plucked string. His 
apparatus consisted of a stretched wire whose ends were 
connected, through a matching transformer, to the input 
of an audio amplifier. Horizontal magnetic fields were 
applied to the wire by setting horseshoe magnets astraddle 
the wire. Plucking the wire in a vertical plane would then 
induce in it an alternating emf which was amplified and 
made audible by the amplifier and a loudspeaker. 

With one magnet placed at the center of the wire, the 
audible tone would consist principally of the fundamental 
frequency of the vibrating string. Two magnets placed so 
that each was located at an antinode for the second 
harmonic, and oriented in opposite sense across the wire, 
would produce a tone consisting principally of the second 
harmonic. In similar fashion, any one of the modes of 
vibration of the plucked string might be made audible by 
placing the equivalent number of magnets across the wire 
in alternate orientation, each magnet being located at an 
antinode for the desired harinonic. 

The writer has developed an elaboration of the original 
design and has called it the ‘‘self-sustaining electronichord.” 
The modification consists of the addition of a transducer 
which is coupled mechanically to the wire and connected 
electrically to the output of the amplifier (see Fig. 1). 


ae | 
AMPLIFIER ourpur 


Fic. 1. Self-sustaining electronichord showing positions of magnets for 
vibration of the string in two loops. 


Thus, when the wire is caused to vibrate in a vertical plane, 
the alternating emf induced in it is amplified and produces 
a similar vibration in the piston arm of the transducer. 
Since this arm is coupled to the wire, energy from the 
transducer is fed back into the wire. If this feedback vibra- 
tion is in phase with the vibration of the wire, a sustained 
vibration of the string results. 

The amplitude of the vibration is controlled by the 
setting of the gain control of the amplifier. The number of 
loops in which the string vibrates is determined by the 
number and the location of the magnets. For example, a 
single magnet placed at the center of the wire will cause it 
to vibrate in one loop, while two magnets placed at the 
antinodes for the second harmonic (as in Fig. 1) will cause 
vibration in two loops. 

The transducer used by the author is the driving unit 
taken from a magnetic speaker. It is mounted above the 
stretched steel wire and mechanically coupled to it by a 
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piece of No. 18 copper wire. The latter is soldered to the 
piston arm of the driver, carried down to the steel wire, 
wrapped around the steel wire and soldered. The location 
of the point of contact with the steel wire is fairly critical, 
as it determines the mechanical impedance matching be- 
tween the vibrating string and the driving unit. In the 
author’s apparatus it is located one-half inch from one end 
of the one-meter vibrating wire. 

The amplifier normally used with this apparatus has two 
stages of voltage gain followed by an output stage which 
delivers up to 8 watts of power. With all of this amplifica- 
tion, it is not necessary to match the wire to the electrical 
input impedance of the amplifier. By matching both input 
and output electrical impedances properly, an amplifier 
consisting of one stage of voltage gain followed by a single 
6F6 power output stage has operated the apparatus very 
nicely. The strength of the horseshoe magnets is not critical, 
and successful operation has been obtained with a wide 
variety of magnets found in the laboratory. 

When the apparatus is connected for the first time, the 
gain control of the amplifier is turned full on and the string 
is gently plucked. If the transducer is connected in proper 
phase for positive feedback, the vibrations of the string will 
immediately build up. If this does not occur, it is usually 
found that the transducer has been connected in inverse 
phase. Its connections to the output of the amplifier are 
simply reversed, and the wire will then vibrate. Whenever 
the number and the location of the magnets are changed to 
cause the string to vibrate in some other harmonic mode, 
the magnet nearest the transducer must always have the 
same orientation; that is, its N pole must always be on the 
same side of the wire. 

The self-sustaining feature has two advantages over the 
original electronichord. First, the selective action of positive 
feedback singles out the desired harmonic more effectively. 
When three magnets, for example, are placed at the 
antinodes for the third harmonic and the string is merely 
plucked, some of the other harmonics are heard along with 
the third. With the self-sustaining action, however, only the 
desired harmonic is produced. Second, the continuously 
sustained vibrations may be made with an amplitude so 
large that a class may see the string vibrate in addition to 
hearing the tone which it produces. 

The author wishes to thank Mr. James W. Grooks for his 
help in building and testing this piece of equipment. 


1N. Urquhart, Am. J. Physics (Am. Physics Teacher) 2, 29 (1934). 


Experiments in Static Electricity 


PauL Roop 
Western Michigan College, Kalamazoo, Michigan 


N the discussion of electrostatic induction the introduc- 
tion of the concept of potential adds materially to the 
effectiveness of the presentation. The simple idea of the 
separation of charges in a conductor under the influence of 
an electric field is inadequate. In the Journal for October 
1940, it was pointed out by the writer that a metal-pointed 





516 


rod, grounded, held near three charged balloons would have 
a potential of about —17,000 volts and that the resulting 
discharge from the point would neutralize the charge on the 
balloons causing them to fall together. Several additional 
experiments, which are best explained using the concept of 
potential, are described below. 


1. If a pointed conducting rod is held by a handle 
made of a good insulator close to three balloons which 
have been charged by rubbing them with fur, no point 
discharge will take place. The potential difference be- 
tween the rod and the balloons is much less than it is 
when the rod is grounded and the electric field at the 
point is not large enough to produce point discharge. 

2. If the pointed end of the rod is now capped by a 
wooden cylinder, grounded, and again held near the 
balloons still no point discharge will take place. The 
potential difference between balloons and rod is the 
same as with an uncapped rod but the wooden cylinder 
prevents any repelled positive ions in the air from 
reaching the balloons and discharging the negative 
charge on them. 

3. If the wooden cap is removed and the grounded 
rod held near the balloons, the latter are quickly dis- 
charged. If, however, a rod having many points at one 
end is grounded and held near, the point discharge will 
be very weak and the neutralization process very slow, 
since the field intensity at each point is much smaller 
than before. 

4. If one rod is held by an insulating handle near the 
balloons and a second rod also held by an insulating 
handle is touched to the first and then separated, an 
induced charge appears on the first rod. The amount 
of the induced charge depends on the potential differ- 
ence between the balloons and the rod. This may be 
measured by transferring some charge from the rod to 
an electroscope with a proof plane. If the experiment is 
repeated with the two rods grounded by touching one 
with the fingers, the induced charge will be found to be 
much greater than before. Using only the simple idea 
of the separation of charges by the inducing action of a 
charged object, the charge appearing on the rod in each 
of the above cases would be the same. Since, however, 
the potential difference is greater in the second case 
than in the first the induced charge appearing on the 
rod is correspondingly greater. 


! Paul Rood, “An experiment in static electricity,"" Am. J. Physics 
8, 320 (1940). 


The Vibrating String Experiment 


HowarpD N. MAXWELL* AND WAYNE GREEN** 
Kalamazoo College, Kalamazoo, Michigan 


RDINARILY, students in our elementary physics 
courses do only three laboratory experiments in 
sound: viz., the Siren Disk, Melde’s Experiment, and the 
Kundt’s Tube. This year, when seeking for an additional 
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TABLE I. Effect of length on the frequency of vibration of a 
stretched wire. 


Length 
ratio 


Reciprocal of 


Percent 
length ratio 


Frequency E 
difference 


ratio 


1.0 1.00 ‘ 0 
0.509 1.97 o 1 
0. 


5 
0.2505 3.99 5 


experiment, it was decided to try the Vibrating Strings 
Experiment using a sonometer. However, since much diffi- 
culty was encountered in attempting to tune the sonometer 
by the usual method of absence of beats, the following ex- 
perimental arrangement and procedure were tried: A 
sonometer having three different types of strings was em- 
ployed. The sounds produced by the strings were picked 
up by a crystal microphone, and converted to a variable 
electrical potential difference which was then amplified 
and applied to the horizontal deflecting plates of an oscil- 
loscope. To the vertical deflecting plates of the oscilloscope 
there was applied a sinusoidal potential difference sup- 
plied by a Hewlett-Packard variable oscillator. The result- 
ing Lissajous figures were then used to determine when the 
strings were properly tuned. 

I. To investigate the dependence of frequency on length, 
the oscillator was first set to 128 vib/sec by holding a 
tuning fork of that frequency in front of the microphone 
and adjusting the oscillator until a stable pattern, an 
ellipse, resulted. Then one string, a meter long, was tuned 
to this same frequency by adjusting its tension. The oscil- 
lator was then set at a frequency of 256 vib/sec with the 
aid of a tuning fork. Using the same wire and tension, a 
movable bridge was adjusted until the wire’s frequency 
became 256 vib/sec. The length of wire was measured with 
a meter stick. The same procedure was repeated for 512 
vib/sec. Some typical data are given in Table I. 

II. The dependence of frequency on tension was found 
by setting the tension on this same string at four pounds 
and adjusting the length with the movable bridge until 
the string was tuned to 128 vib/sec. At this constant length 
the tension on the string was increased until the string was 
tuned to 256 vib/sec. (In both cases the tuning was per- 
formed as previously described.) Table II gives some 
typical data. 


TABLE II. Effect of tension on the frequency of vibration of a 
stretched wire. 


Square root of 
tension 


Square root of 


I E Frequency 
tension ratio 


ratio 


Percent 
difference 


2.00 1.0 1.0 0 
3.93 1.97 2.0 1.5 


TABLE III. Effect of linear density on the frequency of a 
stretched wire. 


Square root of 


Reciprocal of 
density ratio 


Percent 
frequency ratio 


Frequency difference 
256 1.0 1.0 0 
127 2.03 2.02 0. 
121 2.16 2.12 1 


5 
38 
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III. To find the dependence of frequency on the linear 
density, when the length and tension were constant, each 
of the three strings had 16 pounds tension applied. The 
length of one was varied until the string was tuned to a 
frequency of 256 vib/sec. Then using this same length for 
each of the other two strings, the frequency of the oscillator 
was varied until a stable elliptical pattern resulted on the 
oscilloscope screen. The corresponding frequencies were 
read from the calibrated dial of the oscillator. Table III 
shows some typical data. As the data show, the accuracy 
is in the neighborhood of 1.5 percent, which is much better 
than could be obtained by the average student by the beat 
method of tuning. 

With the apparatus used here, the chief error was in 
setting the spring balances at the end of the wires to the 
correct tension. A sonometer using weight pans and pulleys 
to supply the tension might possibly increase the accuracy 
of the results. However, taken as a whole, the experiment 
was very successful in teaching the laws of vibrating strings 
to our class in Liberal Arts Physics. 


* Now at Ohio Wesleyan University. 
** Now at Ohio State University. 


Use of Motion Pictures in Laboratory 
Dynamical Studies 


Carv L. BaILey, JAAN JURISSON, AND M. EUGENE Rupp 
Concordia College, Moorhead, Minnesota 


N student laboratory experiments in dynamics, as, for 

example, Newton’s second law of motion or the con- 
servation of momentum, it is, of course, highly desirable 
to allow motion with the least possible constraint and to 
observe the displacements of moving parts without using 
perturbing markers. This circumstance has led us to use 
motion pictures in the student laboratory for quantitative 
observation of rapid dynamic processes, a technique which 
is used in a well-known textbook.! The photographic 
method combines maximum freedom of movement with 
accurate observation of position as a function of time, if 
the camera frame speed is sufficiently high. 

A double pendulum apparatus, with which we have 
studied the conservation of momentum in the student 
laboratory, is shown in Fig. 1. It consists of two wooden 
blocks, each suspended by a steel rod from a common free- 
swing pivot and counterweighted in such a way that the 
center of mass of each pendulum is at the center of its 
wooden block. A spring is compressed between the two 
blocks, and they are tied together with a thread. When this 
thread is burned, the pendulums recoil. The mass of each 
block can be varied at will. 

The motion is photographed with a motion picture 
camera at 64 frames per second. The camera used in this 
work was a war-surplus gun-point aiming camera, adapted 
only to the extent of adding a viewfinder. Since there is a 
coordinate screen behind the pendulums, the student can 
determine the heights to which the blocks rise by examin- 
ing the film with a magnifying viewer. 


Fic. 1. A frame from movie film of the recoil pendulums flying apart 
after the burning of the binding thread. In this frame the pendulums are 
not near the end of their swing and are not completely stopped by the 
camera. The numbers are retouched for clarity in reproduction; they 
show clearly in the original. 


Since each pendulum receives the same impulse, the 
initial angular momenta are equal, giving 


(11+ Mir)o,= (2+ Mor*)wo. (1) 


Here w; and we are the initial angular velocities of the two 
pendulums, J; and J2 are the respective moments of inertia 
about axes through the centers of mass parallel to the 
actual axis of rotation, M; and M: are the respective total 
masses, and 7 is the common distance from each center of 
mass to the axis of rotation. By the conservation of energy 
we have 


$(11+ Mir’) w:? = Mighi, 
$(I2+ Mor*)wo? = Moghe, (2) 
where /, and he are the respective heights to which the 
centers of mass rise. Combining Eqs. (1) and (2), we find 
fy _ M2(I2-+Mor?) 
he Mi(i+ My,r*) 
The moments of inertia J; and J2 can be calculated from 


the geometry of the pendulums; or if they are constructed 
such that I1«M,r? and I2«Mor*, then Eq. (3) reduces 


(3) 


Fic. 2. A frame from movie film of a ballistic pendulum recoiling 
after being struck by a bullet. The numbers are retouched for clarity 
in reproduction; they show clearly in the original. 
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to 
hi/h2= M?/MY. 


The student may use his data to verify this prediction. 

One may also obtain the initial kinetic energies from 
Eq. (2), and thus study the sharing of energy in recoil. 
Agreement with expected results, to within about one 
percent, is obtained. 

Another laboratory study to which motion picture pho- 
tography lends itself is the measurement of rifle muzzle 
velocities with the ballistic pendulum. Although keenly 
interesting to most students, this experiment is often con- 
sidered too dangerous for class work; the probability of a 
serious accident is never negligible. However, if the experi- 
ment is photographed beforehand by the instructor, the 
student may study the film in safety. Figure 2 is a sample 
frame taken from a 64-frame-per-second photograph of a 
ballistic pendulum. The deflection of the block is measured 


(4) 
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The Demonstration of Nuclear 
Magnetic Resonance 


T does not seem to have been generally recognized that 
the demonstration of nuclear magnetic resonance re- 
quires only quite simple apparatus and should be within the 
scope of the average college physics department. In spite of 
the recent development of the subject there is a widespread 
recognition of its importance as a research technique but as 
yet there seems to be little realization that explanation of 
the basic principles coupled with an actual demonstration 
could form a useful addition to many undergraduate 
physics courses. 

The possible importance of the subject from a didactic 
point of view has been emphasized by Dr. K. K. Darrow in 
an address in which he reviewed the experimental work 
done to date and there is no doubt of the interest aroused 
when the phenomena are displayed and explained. In the 
case of proton resonance no elaborate equipment is needed 
to make possible such a display. 

Of the apparatus required the magnet itself is the major 
component. It should have a field strength of not less than 
4000 oersteds and a gap width of 2 cm or more with pole 
pieces large enough to ensure homogeneity over a volume of 
the order of 1 cm*. If such a magnet is available the re- 
mainder of the equipment may be assembled without 
difficulty. A 300-volt power supply and an oscilloscope are 
required and a small regenerative detector and amplifier 
circuit will need to be built. Constructional details for this 
are identical with those given elsewhere! for a magnetic field 
strength meter based on the same principle. They need not 
be repeated here. Some simplification is possible, however, 
for the present purpose since neither portability nor a broad 
range of application is required. The probe and amplifier 
described in the paper referred to may thus be incorporated 
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in degrees, the numbers being read from the pictures, and 
the height of rise calculated from the deflection and the 
length of the pendulum. 

Other possibilities in dynamics suited to photography 
suggest themselves, such as falling bodies and billiard-ball 
collisions. 

Some readers may believe that this method of laboratory 
study denies to the student the freedom of manipulating 
the apparatus. However, this disadvantage is minimized by 
having the student perform all operations except the actual 
photography (and, in the case of the ballistic pendulum, 
the actual firing). The disadvantage is further balanced 
by the improvement in accuracy which is gained when 
mechanical restrictions on the movement of the apparatus 
are removed. 


1Lemon and Ference, Analytical experimental physics (University 
of Chicago Press, 1946). 
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into a single chassis and the so-called ‘“‘wobbling”’ coils may 
be wound on the magnet pole pieces in any arrangement 
which will give a suitable wobbling amplitude. Also, the 
only sample coil required will be the one appropriate for the 
available magnetic field strength. 

In summary: it has been pointed out that with the con- 
struction of a small electronic circuit and the use of stand- 
ard laboratory equipment it is possible to display proton 
magnetic resonance. 


N. J. Hopkins 
Radiation Laboratory 


McGill University 
Montreal 


1 Hopkins, Rev. Sci. Inst. 20, 401 (1949). 


Velocity of a Longitudinal Wave by an 
Elementary Method 
SSUME that a prismatic bar, of cross-sectional area A, 
density p and Young’s modulus E£, is subjected to a 
sudden compressive force uniformly distributed over the 
face BC (Fig. 1), causing a displacement BB’ =s, and sup- 


BB' KK G 


ce «te H 
Fic. 1. Longitudinal wave traveling along a solid bar. 


pose that the disturbance is transmitted in an infinitesimal 
time ¢ to a section GH, at a distance BG =x with a velocity 
v, so that x =vt. If we call a the average acceleration of BC 
during the time ¢, s=at?/2, whence a=2s/f*. 
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The portion BGHC of the bar has the mass pAx, is com- 
pressed by s, and the net force acting on it is F= EAs/x. If 
we do not suppose v to be infinite, x is infinitesimal, and we 
may consider the displacements BB’, ---, KK’, --+ of the 
sections BC, ---, KL, --- as directly proportional to the 
distances BG, ---, KG, ---, respectively, while the dis- 
placement of section GH is zero. 

Consider the portion BGHC divided into n transverse 
layers of equal thickness and common mass uy, so that 
nu =pAx; then, we may think of F as resolved into n forces 
acting on the m layers and displacing them with the ac- 
celerations a, ---, ax, --+, 0 during the time ¢. These 
accelerations are directly proportional to the respective 
displacements; they form an arithmetical progression and 
their sum is na/2, while the sum of the components of F is 
pna/2=pAxa/2=pAxs/t?, and we can write F=pAxs/t?. 

Equalizing both expressions for F, pAxs/#?=EAs/x, or 
x?/t2= E/p. The quantity in the left side of this equation is 
v, and the formula v=(E/p)! is readily obtained. 


GuSTAVO VILLAR BUCETA 
a de Pinar del Rio 


Cuba 


An Instructional Technique for the General 
Physics Laboratory 


N instructional technique which has been remarkably 

successful in raising the level of student performance 

has been in use in the general physics laboratory at Temple 

University for a period of over twenty years and has not, to 

my knowledge, been described elsewhere. The technique 

involves the use, in a rather unconventional way, of 
undergraduate majors as laboratory assistants. 

General physics laboratory sections at Temple run to 
about eighteen students in size, and meet once a week for a 
period of three hours. To each section of this size there are 
assigned, in addition to the regular staff member in charge, 
two undergraduate laboratory assistants so that the staff- 
to-student ratio in the laboratory is about 1 to 6. 

Experiments are assigned on an individual basis, each 
student being responsible for obtaining his assignment from 
the instructor in charge a week before the date on which the 
experiment is to be performed. No attempt is made to have 
all students in a section perform the same experiment 
during a given laboratory period nor is there any effort 
toward maintaining the same pace for all students. Thus 
each student works at his own pace, being motivated so far 
as speed is concerned only by the requirement that a certain 
minimum number of experiments be completed during the 
semester. 

The uniqueness of the Temple scheme lies in the use 
made of the undergraduate laboratory assistants. These 
assistants, who are junior and senior majors in the depart- 
ment, spend their time quizzing the students, usually in 
groups of two or three at the blackboard, about the 
underlying principles and experimental techniques involved 
in the experiments on which the students are working. The 
questioning is generally of the Socratic type, the assistant 
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attempting to draw out of the students a connected account 
of the particular experiment. 

The members of the general physics class are told, of 
course, to expect these oral laboratory quizzes but what 
they do not know is when the quiz will occur in any par- 
ticular case. The student may be called to the blackboard at 
the beginning of the laboratory period before he has even 
had an opportunity to inspect his equipment or he may be 
called up when he is part way through the experiment or 
after he has completed it. Depending upon his preparation 
and his ability to answer the assistant’s questions, the quiz 
may be short, lasting only fifteen or twenty minutes, or it 
may take up the full three-hour laboratory period. After 
one experience with a three-hour quiz on an experiment or, 
if he is wise, before such an experience, the student will 
acquire the habit of coming to the laboratory thoroughly 
prepared on the experiment to which he has been assigned. 

The advantage of this scheme in increasing the effective- 
ness of the general physics laboratory is fairly evident. The 
student gets much more out of his experimental work be- 
cause he comes to the laboratory properly prepared and 
does not, therefore, perform his experiments in the more 
common cookbook style of alternating a reading of the 
instructions with an adjustment of some control or the 
reading of a scale. The advantages accruing to the depart- 
mental majors who serve as laboratory assistants are, how- 
ever, equally important. They feel that they are part of the 
instructional team of the department and thus acquire, 
early in their work, a professional awareness as physicists. 
An unusually large number of those who serve as under- 
graduate laboratory assistants continue their work in 
graduate school and eventually receive doctor’s degrees. 

The number of undergraduate assistants required at 
Temple to implement the scheme described is about equal 
to the number of departmental majors so that all majors 
who desire to work in this way may obtain employment. In 
general, each assistant devotes nine hours a week to this 


work and receives pay of approximately a dollar and forty 
cents per hour. 


BERNARD B. WATSON 
Division of Higher Education 


U. S. Office of Education 
Washington, D. C. 


Airy’s Theorem and the Improvement of Clocks 


HE article by Mr. Gaehr on Airy’s theorem in the 
American Journal of Physics (16, 336 (1948) ] con- 


tains some slips which, in the interests of accuracy, call for 
correction. 


The author of the Theorem is referred to on page 336 
as Lord Airy. Airy was never a peer, but died in his 90th 
year as Sir George Airy. 

Lord Grimthorpe, who was a peer, published the last 
edition of his treatise in 1903, not in 1872 as stated on 
page 336. 

The statements about Harrison’s chronometer, be- 
ginning at the foot of page 337, are also wide of the 
mark. The escapement is described and illustrated in 
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Britten’s Old Clocks and Watches (Spon, 6 ed.), page 
375. It was a crown wheel and verge escapement ex- 
actly as shown in Gaehr’s Fig. 1 and without question 
far inferior to the detached escapement invented by 
Mudge. 

On page 341, Mr. Gaehr uses a circle diagram to 
derive the change in time due to the impulse—in other 
words, the escapement error. I cannot quarrel with the 
general method, since I have published a derivation of 
Airy’s integral by a similar process (Science of Clocks 
and Watches, Pitman). Mr. Gaehr’s result does not 
agree with Airy; and in particular none of his expres- 
sions for AT on page 341 has the dimensions of a time 
as it should. 

The failure of Huygens’ cycloidal pendulum is 
ascribed by Mr. Gaehr to its introducing escapement 
errors. The ground for this assertion is not clear, since 
Mr. Gaehr’s expression for the escapement error has no 
term which a cycloidal pendulum could affect; and 
furthermore, Airy explicitly stated that his results 
applied to a cycloidal pendulum only. 


A. L. RAWLINGS 
33 Maynard Avenue 


Waterbury, Connecticut 


Airy’s Theorem and the Improvement 
of Clocks: a Reply 


INCE the appearance in this magazine last September 

of my article “‘Airy’s Theorem and the Improvement in 

Clocks,” a number of letters pointing out errors have come 

to me, among them the letter of Mr. A. L. Rawlings above. 

I should be grateful if, in the interests of accuracy, the 
corrections listed below be published. 


As may be inferred from my article, I did not have 
access to many books on the theory of clocks, although 
I have hunted through four large libraries. None of 
them included the excellent and recent book on Science 
of Clocks and Watches by A. L. Rawlings (Pitman, 
rev. ed., 1948), which should therefore have been in- 
cluded in the list of recommended references. Mr. 
Rawlings writes that he has there given a demonstra- 
tion of Airy’s equation by a method similar to mine. 

Lord Grimthorpe’s latest edition appeared in 1903, 
instead of in 1872. 

Lord Airy should read Sir George Airy. 

Page 338, first sentence of first column: the escape- 
ment of Harrison’s chronometer was not better than 
Mudge’s. 

Page 339, second column, line 10, et seg.: The whole 
duration of an impulse was about one-fourth of the 
swing from left to right. . . . the work of each im- 
pulse was calculated to be 25.6 ergs. In the 13th line of 
the first column of page 340, change 2.5 ergs to 20.8 
ergs. 

Near the bottom of the first column of page 341, and 
first line of the second column, the expressions 6;/z, 
62/x, (0:—62)/x should all be multiplied by 7/2. In the 
equations for AT, read 2-AT/T. 


. 
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In the second column of page 342, I wished to say 
that in eliminating one error, clock makers often 
introduced another, giving Huygen’s cycloidal pendu- 
lum as an example. One writer stated that Huygen’s 
pendulum failed because of escapement errors; another 
writer gave the friction of the cycloidal cheeks as the 
cause; and still another writer said that the bob was 
much too large to be considered as a point. 


PauL F, GAEHR 
Wells College 


Aurora, New York 


An Interesting Exhibit 


HE writer was recently asked to make up a group of 
exhibits for a college ‘‘open house” day. One of the 
exhibits which attracted considerable attention of both 
technical and nontechnical guests consisted of a fluid 
rotating at a speed of several revolutions per minute with- 
out a mechanical propellor. The fluid was propelled by the 
action of a vertical magnetic field on moving ions. 
The exhibit is shown in Fig. 1. The field was provided by 


Fic. 1. A current-carrying electrolyte in a magnetic field. 


a large magnetron permanent magnet M. A glass container 
G, approximately 15 cm in diameter and 5 cm deep, was 
placed between its, poles and contained a quantity of 
concentrated copper sulfate solution S. A copper wire W 
was suspended from the upper pole of the magnet and 
allowed to dip down into the container. A thin copper band 
B was inserted along the inner periphery of the container 
and a potential of 10 volts was applied between these 
electrodes. The direction of the current was reversed from 
time to time, which, of course, caused the direction of 
rotation of the fluid to change. Two copper electrodes are 
used for simplicity and to avoid the occurrence of chemical 
changes in the system. 

It is hoped by the author that this presentation will be of 
use to some poor soul faced with the dilemma of providing 
an entertaining public exhibit, while simultaneously 
satisfying technical people that he is not merely deceiving 


the public with spectacular nonsense. 


J. L. RYERSON 
Evansville College 
Evansville, Indiana 


Can We Fly to the Moon: Errata 


EVERAL errors appeared in the paper by Dr. Himpan 
and the writer entitled ‘‘Can We Fly to the Moon.’”! 
The list of corrections is as follows: 
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Page 257, Eq. (2) should read 
he = Wele—We(Ga— T/G) In(Ga/Ga—T) — fgzte. 


Page 258, left-hand column, lines 32 and 33 should 
read 


Ga=5X10°X5X10°N.=:-- 


Page 259, left-hand column, lines 14 and 15 should 
read 


Vet = 3540 m/sec. 


Page 259, left-hand column, line 32 should be 
a trifle more than the 5.07 X107 kg mthatis.. . 

Page 260, left-hand column, the line preceding Eq. 
(11) should read 
By integration of Eq. (10) it isfound .. . 

Page 260, Fig. 4: 

Values of the ordinate G must be multiplied by 10-°. 
Thus the ordinate-values shown in ascending order 
are: 10-4, 10-8, 10-2, 1071, 10°, 10', and 10%. 

Page 261, right-hand column, line 9 of text should 
read 
with an equivalent top speed of .. . 

Page 262, left-hand column, line 5 of item (b) should 
read 
The increase is proportional to (H7r)}. 


R. REICHEL 
Villa La Joliette 


St. Cloud, France 


1 Himpan and Reichel, Am. J. Physics 17, 251 (1949). 


The Dimensions of Angular Displacement 
and Torque 


N a recent article in this journal by P. Moon and D. E. 
Spencer! it was recommended that the linear dimension 
be split into two in order to provide a 1:1 relationship 
between concept and dimension in the case of angular dis- 
placement, moment of a force, and torque. I would like to 
suggest that the assignment of primary dimension to 
angular displacement (1) is more logical, and (2) is neces- 
sary in order that the dimension of a physical quantity 
yield correct information about it. 

Concerning item (1) it is generally contended that angu- 
lar displacement is a numeric since angle is the ratio of two 
lengths. But the relation @=s/r is a special case of the more 
general one 0,/02=5;/s2, or 0=ks/r where k is a propor- 
tionality constant and 6 may be measured in any unit. 
It is @/k which is dimensionless, not 6. To assume that 
choosing 6 in radians makes k dimensionless, then to use 
this assumption to prove that the resulting relation implies 
that @ is dimensionless, is circular reasoning. If a linear 
unit of such length were chosen that the gravitation con- 
stant G would be numerically 1, and if it were then as- 
sumed that this made it a numeric, then the resulting rela- 
tion between /, m, and t implied by the equation F= m,m2/d? 
would destroy most of the utility of dimension theory in 
mechanics. 

Another argument for assuming the radian to be dimen- 
sionless is that the Maclaurin’s expansion of sin 6 necessi- 
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tates that @ be a numeric. In Eshbach’s Handbook of 
Engineering Fundamentals, page 3-37, is the following 
statement: “Functional forms (like trigonometric or 
exponential functions) can have only numerical arguments; 
sin E, where E means energy, or voltage, has no sense 
because it violates dimensional homogeneity. This can be 
seen by expanding the function into the equivalent power 
series, 


sin E= E/1!— E3/3!4+E5/S!—---, 


where each term has different dimensions.” This argument 
is not valid. The n term of the Maclaurin expansion of a 
function is f"-!(0)@"—-!/(m— 1)! For any function the deriva- 
tive f("-)) (6) =d"— f(@) /d6"~ has the dimension of (6) /6"-1 
for any value of 6, and hence the general term has the 
dimension of [d"~1f(6)/d@"-*]e"—, that is, of f(6); this is 
independent of 7, so the series is dimensionally homogene- 
ous. A similar statement holds for the Taylor expansion. 
Concerning item (2), I would delay the discussion of the 
moment of a force until work and potential and kinetic 
energy are defined. Then a qualitative definition of the 
term can be stated, as usual, as the ability of a force to 
produce rotation of a body about an axis. A quantitative 
definition can then be arrived at by considering two parallel 
but oppositely directed forces F, and F2 which act per- 
pendicularly on a bar of constant mass (concentrated at 
points or distributed) which is free to turn about one end, 
the points of application of the forces being distant 7; and 
re units, respectively, from the axis. Then, instead of 
giving experiment as the basis for the relation Fir;= Fore 
(as is the practice in texts), the result can be obtained from 
the law of the conservation of energy, thus reducing by 
one the number of assumptions in elementary mechanics. 
Let XY be a line in the plane of the forces and passing 
through the axis, and @ the smallest positive angle through 
which the bar must be rotated in the direction of F, to bring 
it into coincidence with line X Y. Suppose the bar is in 
equilibrium and the magnitudes of the forces are constant. 
Then F; alone is potentially capable of rotating the bar 
through the angle @ (the direction of the force remaining 
perpendicular to the bar) thereby adding an amount of 
kinetic energy which must equal F,(r10/K) regardless of 
the masses involved or the final velocity. In like manner 
the potential energy of F2 with respect to the line XY is 
— F,(r26/K). Since the bar is assumed to be motionless with 
respect to the line X Y, its kinetic energy is zero, and from 
the conservation of energy the algebraic sum of the poten- 
tial energies is zero; i.e., Fir,10/K—For20/K=0, or Fini 
= Fyre. It should be emphasized that the direction of the 
distance through which the force acts at any instant is the 
direction of the force, and not that of the normal to the 
force; that is, the quantity Fir,0/K is work, and is meas- 
ured exactly as work is always measured. We now take 
as our measure of the moment of either force about the 
axis the energy which it (alone) is potentially capable of 
adding to the system per unit angular displacement. In 
any system of units in which the radian is the unit of 
angle, the moment M of a force is equal to Fr energy units 
per radian. In the gravitational system the unit is ft-lb 
per radian. The writing of the energy unit as lb-ft is un- 
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necessary and misleading, and the omission of the angular 
displacement dimension makes the meaning vague. 

The extension of the definition to a body in rotation 
readily leads to the conclusion that torque is the rate of 
change of kinetic energy per unit change in angular dis- 
placement (=d(}IJw*)/d@). For constant masses torque 
and force are thus the rates of change of rotational and 
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translational kinetic energy, respectively, the former with 
respect to unit angular displacement, and the latter with 
respect to unit linear displacement (=d(4mv*) /ds). 


Rosert S. SEAMONS 


Yakima Valley Jr. College 
Yakima, Washington 


1 Moon and Spencer, Am. J. Physics 17, 171-177 (1949). 
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Book Reviews 


Principles of Electricity and Magnetism. Second Ed., 
GayLorp P. HARNWELL. Pp. 670, Figs. 398, 1523 
cm. McGraw-Hill Book Company, Inc., New York, 
1948. Price $6.00. 


The first edition of this book appeared in 1938, and was 
reviewed in this Journal (7, 203 (1939)). The outline of 
Professor Harnwell’s text remains essentially unchanged in 
the new edition. However, the significant developments 
during the past ten years in the fields of atomic and nuclear 
physics and in the use of electromagnetic waves have 
influenced the emphasis on and the selection of topics in 
this new edition, e.g., the accounts of electronics and of 
radiation have been considerably expanded in the direction 
of higher frequencies, and a brief account of generators at 
higher frequencies and of resonant cavities and wave guides 
is now included, as well as some material on semiconductors. 

Throughout the book, the treatment is developed from 
the point of view of the physicist ; the fundamental aspects 
are emphasized, and the experimental basis of the con- 
siderations leading to Maxwell’s equations are clearly dis- 
played. The student using this text not only will learn about 
electricity, but will get an integrated picture of a great part 
of modern physics. Many applications from neighboring 
fields are given; those bearing on topics of current interest 
are stressed somewhat, as for example accelerators or the 
nuclear induction experiments. 

The text starts with electrostatics; vector analysis is 
freely used in the formulation of the general concepts, and 
complex variable methods are used in the discussion of 
special problems. The discussion of the electromagnetic 
effects is based on Ampére’s experiments. The electrical 
charge is considered as the fundamental electrical entity 
here as well as in the subsequent treatment of magnetism. 
Magnetic poles are used only later, as an aid in calculation. 

The chapter on nonohmic circuit elements is greatly 
enlarged and contains material not found in other texts. 

' The chapters on a.c. and vacuum tube circuits, and on 
coupled circuits and lines, contain a wide selection of 
topics; emphasis is placed on an understanding of phe- 


nomena, rather than the description of devices or their 
functioning. A concise chapter on the basic types of 
electrical machinery is also included. 

A summary of vector analysis and of frequently used 
material on differential equations is included in an ap- 
pendix. The mks system of units is used principally, but the 
older systems of units are mentioned and used occasionally. 
The text contains a wealth of problems, many of them new, 
and most of them having real physical meaning, like the 
problems in magnetism where pole pieces used in molecular 
beam experiments are discussed. 

The preface states that the text is intended for students 
having a sound preparation in general physics and in 
mathematics up to and including differential equations. 
Some previous acquaintance with topics usually covered in 
advanced calculus, seems also indicated, since div grad@ is 
used on p. 30 and log(x+iy) on p. 40. This may make the 
book a bit hard for many undergraduate students. 

Notwithstanding the use of mathematical tools, where 
they are a help in clarifying an argument, the physical 
problems and concepts and the experimental physicists 
approach always remain preeminent. 

In view of the many valuable features of Professor 
Harnwell’s book, and of its whole approach to intermediate 
teaching, it should be underscored that it could also serve 
in courses for students lacking a certain degree of mathe- 
matical maturity. This is possible because the more 
advanced mathematical tools are actually used only in the 
more theoretical chapters on electrostatics, electromag- 
netism, and radiation. For example, the more qualitative 
chapters on dielectrics and magnetic materials, or the 
chapter on photoelectric effect and related phenomena, 
with its discussion of thermionics and of semiconductors, 
could provide an excellent source of collateral reading 
material for any type of course in electricity. This book 
might prove to be an excellent text for an advanced physics 
course serving primarily the needs of engineering students. 
It certainly should provide a source of stimulating reading 
for any beginning physicists or their teachers. 


F. T. ADLER 
University of Wisconsin 
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Elementary Photography. Second Ed., GiLrorp G. 
Quar_es. Pp. 345, Figs. 42, Plates 78 (2 in color), 
649} in. McGraw-Hill Book Company, Inc., New 
York, 1948. Price $4.50. 


The general nature of this interesting elementary book on 
photography is very well described by the author’s remarks 
in the preface to the first edition. It is designed to be 
intermediate in character between very elementary books 
suitable for high school use and other quite advanced ones 
too technical for general reading. It was written primarily 
for use as a textbook for students without the background 
of physics and chemistry, but also with the desire to make 
it a useful tool for the average amateur photographer, or 
beginner. The author has succeeded very well in his 
project. 

The range of subject matter covered is unusual. To 
quote, the author ‘‘knows of no other book on the general 
subject (photography) which treats of pictorial composi- 
tion, projection control, picture making with paper nega- 
tives, toning, hand coloring and natural color in such a 
manner as to be usable by any serious worker, regardless of 
previous training.” 

In the preface to the new (second) edition the author 
states that the revision has not been made in order to 
incorporate new material, but rather to improve its 
“teachability” by making certain additions and rearrange- 
ments. The principal changes are the addition of study 
questions at the end of each chapter, and a rearrangement 
of the chapters. The study questions are quite numerous, 
and well chosen, and should direct the efforts of the reader 
in a very profitable manner. The rearrangement is typified 
by the fact that Chapter II is a short discussion of the most 
elementary optics of the camera, further discussion follow- 
ing in later chapters, and is followed immediately by ‘‘Film 
Sensitivity and Exposure,’”’ Chapter III, and by ‘“‘Develop- 
ment of Negatives,’’ Chapter IV, the main idea being to 
lead the student as soon as possible into the actual photo- 
graphic process. 

The subject matter which of necessity appears in any 
book on photography, however elementary, is very well 
presented. The discussion is always simple, complete 
enough to satisfy any beginner, never stilted. Chapters VII 
and VIII complete the discussion of optics, and include the 
aberrations of lenses and such items as the view finder and 
the range finder. Chapter X deals with the color response of 
the various common emulsion types and the use of correc- 
tion filters and contrast filters. Chapter XV gives the 
elements of projection printing, and Chapter XVIII 
describes the art of ‘‘finishing’’ the print. 

Certain chapters are unusual in an elementary treatise on 
photography. One usually hunts for such material in 
separate books. Here one finds a very nice discussion of 
“Composition,’’ Chapter XVI, the fine points of enlarging 
in “Projection Control,’ Chapter XVII, and ‘Picture 
Making with Paper Negatives,’’ Chapter XIX. ‘‘Natural 
Color Photography,’’ Chapter XXII, however, receives a 
rather sketchy treatment. 

The book ends with a suitable bibliography, Appendix A, 
and a formulary, Appendix B, which contains the essential 
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and standard formulae indispensable to any photographer, 
amateur or otherwise, who really practices the art. After 
reading the book, one has the feeling that here is a job well 
done. 


I, C. CorNoG 
University of Pennsylvania 


The Structure of Matter. Francis OWEN RICE AND 
Epwarp TELLER. Pp. 361+-xiii, Figs. 69, Tables 13, 
6X9} in. John Wiley & Sons, Inc., New York, 1949. 
Price $5.00. 

This book, the first of Wiley’s Structure of Matter Series, 
is a valuable addition to the already available books on the 
application of quantum theory to the interpretation of the 
behavior of matter. The reviewer’s very favorable reaction 
is based on the style of the book, rather than on the in- 
clusion of any new material. The book makes no pretense 
of adding to the known facts, but it does present these facts 
in a way that should appeal to a large audience, surely 
eager for such a style of exposition. 

Many textbooks and reference books require quite a lot 
of between-the-lines reading and interpretation. One of the 
main functions of a teacher is to simplify this task for the 
student ; and one of the chief advantages of taking a formal 
course in a subject, rather than depending on self-study, 
is this service rendered by the teacher. Much of this is 
accomplished in the classroom by what the reviewer likes 
to call the ‘‘hand-waving” technique. One talks in a general 
way and tries to create mental images that illuminate and 
clarify the subject, and make it more tangible and friendly 
and less formal and obscure. The Structure of Matter is one 
of the rare books that offer this service to the reader 
of other texts and reference books dealing with similar 
material. In addition, it presents a rather extensive set of 
phenomena associated with the properties of matter, all 
presented with clear qualitative discussions that give the 
book its particular quality. It thus serves as a useful 
source of material, as well as an interpreter of theories. 
The authors state as the purpose of the book the wish ‘‘to 
acquaint the reader with the scope of the phenomena that 
can be explained with the help of quantum mechanics.” 
In attempting to do this, most of the formal mathematical 
analysis is omitted and the qualitative results of the theory 
alone are presented. One may, at first thought, regret the 
omission of the mathematical details; but, as ope becomes 
better acquainted with the book, the advantages become 
evident. 

An important technique of thinking about physical 
phenomena is the use of reasonable estimates to suggest 
the behavior to be expected. The Structure of Matter makes 
frequent use of this technique, reminding the reader of its 
power. Unfortunately, the first example gives the precise 
rather than an approximate answer! This result occurs in 
the calculation (p. 7) of the energy of the ground state of 
the hydrogen atom. The uncertainty relation in the form 
ApAq~h/2r is used in a clever way to give the result. The 
reviewer objects to the combination of an accurate answer 
coupled with a lack of emphasis on the fact that a mere 
estimate is being made. If the authors had used ApAg 
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~h/4x, a common form, the result would have differed 
from the true value, but would have given its order of 
magnitude. The approximation would have emphasized 
the limitations as well as the advantages of the method. 
It is interesting to note that a similar calculation for the 
ground state of the linear harmonic oscillator, based on the 
form ApAg~h/2z, gives hvo, while the form ApAqg~h/4r 
leads to $hvo, the correct value. In later applications the 
authors make the approximate nature of their results more 
apparent. 

The general topics included are suggested by the chapter 
headings: Introduction (an introduction to quantum 
principles), The Hydrogen Atom, The Periodic System, 
Motion and Position of Nuclei in Molecules, Atoms and 
Molecules in Electric Fields, Van der Waals Attraction 
Forces, The Chemical Bond, Forces in the Solid State, 
Magnetic Properties of Matter, Molecular Vibrations, 
Electronic Spectra, Nuclear Chemistry, State of Matter in 
Stars. This list shows that the contents include the bulk of 
the results of the ‘‘classical’’ (that is, the pre-1940) quan- 
tum theory. An amply diverse lot of examples is given in 
each chapter. A specialist in a field may regret the omission 
of certain examples and may find minor errors or indirec- 
tions. For example, this reviewer was a bit unhappy about 
several aspects of the treatment of atomic spectra, and 
would have preferred a closer adherence to the conven- 
tional notation. The authors do not use a subscript to 
designate the value of J, and they use u and g to label odd 
and even states, contrary to the more common use of the 
symbol ° for odd states and omission of any symbol for 
even states. 

The use of diagrams has not been exploited in this book 
as fully as it might have been. While it is difficult to con- 
ceive a really adequate graphical scheme conveying to the 
reader the spatial form of wave-functions, and while there 
is need for the development of such a system, this reviewer 
feels that Rice and Teller have not improved over other 
schemes. One final objection is to the complete lack of 
references to original and comprehensive treatments of the 
various topics. Although the inclusion of such references 
might have required too many footnotes or too long a 
bibliography, this book should arouse such interest in the 
various examples, that the reader will almost certainly 
want to carry on his study and to refer to the more com- 
plete treatments. A bibliography would then be a great 
help. 

These criticisms are minor and should not mask the very 
many worthy features of the book, which far outweight 
such picayune objections. The Structure of Matier should 
be a useful book for nearly any physicist, regardless of his 
special field. It will remind him of the facts he has learned 
and may have forgotten; and it will probably introduce 
him to other fields of which he is unaware because removed 
from his main interests. 

Although students with little training in quantum me- 
chanics can get a great deal of good out of this book, it is 
really addressed to those who are already familiar with 
quantum mechanical techniques. Others may not be able 
to appreciate fully many of the points that are discussed. 

The Structure of Matter is basically sound, with only a 
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few minor deviations from customary views. It offers a 

new approach. It is a source of much useful information. 

It should prove to be a book that is referred to constantly. 
A. T. GoBLE 
Union College 


Elektrizitatslehre. R. W. Pont. Pp. 301, Figs. 497, 1725 
cm. Springer Verlag, Berlin-Géttingen-Heidelberg, 1949. 


The review of a German textbook in an American 
journal is justified only when it contains certain outstand- 
ing features, or gives opportunity to general reflections on 
educational problems which are under discussion on this 
side of the Atlantic. Both can be said about Professor 
Pohl’s Electricity, which forms the second volume of an 
Introduction to Physics, the first volume dealing with 
mechanics and the third with light. 

Pohl’s Electricity was published in an English version in 
1930, and issued several times thereafter (Blackie and Son 
Ltd., London and Glasgow). A comparison of the new book 
with the English translation shows considerable changes 
and in some instances definite improvements. 

The most attractive feature of the book has always been 
the great number of excellent figures, representing shadow 
projection pictures of instruments and apparatus as they 
are shown in Professor Pohl’s classes for beginners at 
Géttingen university. In central European universities it 
is the common practice in the first physics course to pro- 
duce ample experimental evidence through demonstration 
experiments. These courses are called ‘experimental 
physics,” and the name was often understood as an obliga- 
tion to proceed from experiment to experiment. In this way 
the lecture course became more and more a continuous 
show, and the intellectual level of instruction necessarily 
suffered. Lack of time hindered the lecturer to give theo- 
retical analytical demonstrations, or to point out the 
historical development of physics. The wealth of experi- 
ments in Professor Pohl’s lectures, I was told, was truly 
appreciated only by advanced students, who frequently 
returned to his classes on account of the excellent demon- 
strations of physical phenomena. The success of Pohl’s 
books was also greater among instructors and teachers in 
secondary schools (gymnasia), than among university 
students. Books and demonstration technique were also 
adapted for more elementary school physics courses. 

The new edition of. the Elektrizitdtslehre corrects to a 
certain degree the fault of neglecting mathematical demon- 
strations. Historical discussions are now more numerous, 
but still appear somewhat superimposed on the subject- 
matter. The book should now have greater usefulness for 
the younger student, but will still appeal more to the 
teacher, I believe. It is paradoxical that Pohl’s book can 
be considered as highly didactical, although it is not really 
pedagogical. 

At the time of its first publication Pohl’s Electricity was 
outstanding through two new features: the use of the 
practical electrical units throughout, and the change in 
the sequence of experimental demonstrations, on which 
Pohl prides himself particularly: ‘‘The starting point of our 
experiments are not rubbed amber and bits of paper, but 
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chemical sources of electricity as dry batteries and storage 
cells. . . . We measure electrical quantities in electrical 
units, as multiples of a unit-current (amp) and of a unit- 
tension (volt).”” Pohl considers electric current and voltage 
as terms with which the student is familiar, and uses 
ammeters and voltmeters from the beginning in the same 
way as one uses clocks and watches in mechanical experi- 
ments, i.e., without a previous investigation of their 
construction. 

For the consistent use of practical (now international) 
electrical units, and the rational system, Pohl’s book did 
pioneer work. He makes reference to the electrostatic and 
the electromagnetic systems only in passing. There is 
little doubt that for the engineer and applied physicist the 
use of practical electrical units of the meter-kilogram- 
second system (joule, watt, volt, ampere, ohm, coulomb, 
henry, farad, weber (=volt sec)) has great advantages. 
The two additional vacuum-constants, permittivity 
€o= 8.86 X10-" and permeability wo=1.256X10-§, cannot 
be considered to present serious difficulties. The appear- 
ance of the factor 1/47 in Coulomb's laws is compensated 
by the disappearance of this factor from other equations, 
e.g., for the capacitance of the plate condenser, and the 
energy of the electric and magnetic fields. 

However, from the didactic standpoint arises the ques- 
tion that we have to explain why the numerical values of €o 
and uo can be “theoretically”’ expressed as 1/(44% 9X 10"2) 
and 4rX10-?. Had the practical units been introduced 
without any reference to the absolute units, and had the 
constants arbitrary numerical values, there would be no 
difficulty. As it is, the thoughtful student wants an ex- 
plicit explanation, which he would not find in Pohl’s book. 
There is no American or English physics textbook, say 
on the sophomore level, which attempts an exposition of 
this part of an electricity course in a comprehensive and 
comprehensible way. To use different units and systems 
in different parts, and then to connect them in an appendix 
or a table, as is often done, is a sure way to bury the 
problem. 
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, One third of Pohl’s book is devoted to the principles of 
electric conduction, in all three states of matter; valuable 
material is included in this section which in America would 
be mostly treated in lectures on modern physics. There are 
also special chapters on Radioactivity, on Electric Waves, 
and on “Relativity as a Principle of Experience.” The book 
can well be called a complete work, which shows Professor 
Pohl’s ingenuity as an experimentalist, and also his ex- 
perience as a teacher by the lucid explanations of ques- 
tions of detail. 

There are some noticeably eccentric features in the book. 
Pohl does not use the name newton for the mks unit of 
force, but speaks about Grossdyn (Great-dyne). Instead of 
joule he writes for the energy unit Grossdyn-meter. In the 
English translation from an earlier German edition the 
force unit is referred to as ‘‘the nameless force unit,”’ but 
the translator of a new English version would be well 
advised to conform with English-American and _inter- 
national usage. This would also apply to the translation of 
the German Spannung either by tension or voltage, but not 
by pressure. 

More odd even than the apparent avoidance of newton 
and joule (but not of watt), is the introduction of new 
names for improvised units which are used for a certain 
demonstration—here insignificant—at the end of the 
book. Pohl gives them everyday names, namely Fritz for 
Widerstand (resistance), Max for Spannung (voltage) and 
Moritz for Strom (current). With the choice of the names 
Max and Moritz Professor Pohl refers to the classic 


German comic children’s book; he is evidently joking, and 
I think everybody can well see the humor of it. However, 
if we are allowed to recognize the great-in-the-small, does 
not Professor Pohl’s attitude in respect to the naming of 
physical units in a scientific text show signs of a fallacious 
isolationism and parochialism in contemporary German 
scientific education? 


Otto BLixn 
University of British Columbia 
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New Members of the Association 


The following persons have been made members or junior members (J) of the American Association of Physics 
Teachers since the publication of the preceding list [Am. J. Physics 17, 456 (1949)]. 


Baber, Stanley (J), 5002 S. ‘‘J’’ St., Tacoma 8, Wash. 

Badgley, Ralph Emerson, Eastern Oregon College of Edu- 
cation, La Grande, Ore. 

Beers, Yardley, 4 Ploughman’s Bush, Riverdale 63, N. Y. 

Bolen, Max C., Physics Department, The James Millikin 
University, Decatur, II. 

Bragg, Robert H. (J), 5623 Wabash Ave., Chicago 37, Ill. 

Burns, Jay III, 1712 W. Estes Ave., Chicago 26, Ill. 

Carlson, Henry Arthur (J), Rte No. 1, Blairstown, N. J. 

Causey, Robert L. (J), 936 Covington Ave., Bowling 
Green, Ky. . 


Coe, John E., 2024 Sunnyside Ave., Chicago 25, IIl. 

Coyle, Thomas (J), 910 Franko St., Throop, Pa. 

Darlington, Lillian E. (J), 6630 S. Peoria St., Chicago 21, 
Ill. 

Desjardins, Lucien Henri, 52 Assumption Ave., Worcester 
6, Mass. 

Ellison, Robert William (J), 717 S. 20 St., Lafayette, Ind. 

Forrester, Alvin Theodore, University of Southern Cali- 
fornia, Los Angeles 7, Calif. 

Garr, Carl Robert (J), 127 Pittsburgh St., Ravenna, O. 

Goldstein, Herbert, 39A Lee St., Cambridge 39, Mass. 
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Harper, William Irving, St. John’s College, Annapolis, Md. 

Hinkle, Mark L., Jr. (J), 603 West Fourth St., Marion, 
Ind. 

Honnold, Vincent Richard (J), Box 403, Notre Dame, Ind. 

Howard, Francis Joseph, 109 Torrey Pines Homes, La 
Jolla, Calif. 

Kelley, James B., 363 Clarendon Rd., Hempstead, N. Y. 

Kyle, Garland D., Department of Mathematics and 
Physics, Arkansas Agricultural, Mechanical and Normal 
College, Pine Bluff, Ark. 

McCormick, Thomas Herbert, Box 5, Tuskegee Institute, 
Ala. 

McHenry, James R. (J), 201 West Wood St., Norristown, 
re 

Moe, David Eimon (J), Physics Department, Washington 
University, St. Louis, Mo. 

Naegele, Eugene Lee (J), 4511 W. 228th St., Cleveland 16, 
‘Ohio. 

Peoples, J. Fred (J), 2734 N. Penn St., Indianapolis 5, Ind. 

Petrauskas, Alexander Anselmus, University of Notre 
Dame, Science Hall, Notre Dame, Ind. 

Pflasterer, Edward (J), 706-18th Ave. South, Nashville 4, 
Tenn. 

Pursell, Lyle Eugene (J), 323 West Lutz, West Lafayttte, 
Ind. 

Pyron, Berry Owen (J), 922 N. Highland Ave., N.E., 
Atlanta, Ga. 

Quinsler, Phillips Brooks, Jr. (J), 3840 Calle Guaymas, 
Tucson, Ariz. 

Rogers, Edward De Lancey Francis, Department of 
Physics, Stanford University, Stanford, Calif. 
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Sampias, Olga Ilona (J), Box 10, Barrackville, W. Va. 
Schein, Marcel, 1312 East 56th St., Chicago 37, Ill. 


Smith, Philip Bartlett, 507 E. Springfield Ave., Champaign, 
Ill. 

St. Clair, William Feagin, Box 97, Station A, Hattiesburg, 
Miss. 

Strowman, Philip M., 504 W. 112th St., New York 25, 
te. ee 

Swann, Edwin Garner, Western Geophysical Company, 
Party 51, Vernal, Utah. 

Terhune, John Alexander, 1250 Santa Barbara St., San 
Diego, Calif. 

Thayer, C. Richmond (J), Box 120, New Florence, Pa. 

Thiele, Frederick William, Box 61, City College, 139th St. 
and Convent Ave., New York, N. Y. 


Tubbs, Eldred Frank (J), Powers Road, Orchard Park, 
N. Y. 

Urban, Paul J., 105 West 16th St., Hays, Kans. 

Walker, Elbert Sterling (J), Rickwood Field, Birmingham 
8, Ala. 

Weaver, Ted L. (J), 318 N. ist West, Salt Lake City, 
Utah. 

Weddell, James Blout (J), 85 Sussex Ave., Morristown, 
N. J. 

Wheeler, George Frederick, 551 
Atlanta, Ga. 

Woods, Calvin Robert (J), 623 Edel Ave., Maywood, N. J. 

Wright, David Umstead, Jr., 7015 Torresdale Ave., 
Philadelphia 35, Pa. 
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Science Courses in General Education 


College science courses designed for general education 
exist in more than half of America’s four-year colleges 
and universities according to a survey recently completed 
by Robert A. Bullington of MacMurray College, Jackson- 
ville, Illinois. The study, which formed the basis of a Ph.D. 
dissertation at Northwestern University, was carried out 
under the sponsorship of the Cooperative Committee on 
Science and Mathematics Teaching of the American Asso- 
ciation for the Advancement of Science. 

Information concerning science offerings obtained from 
720 institutions revealed that 293 or 41 percent offer no 
science course specifically designed for students who are 
not majoring in science while 427 or 59 percent offer general 
education science courses of various types. 

Detailed information concerning the nature and objec- 
tives of 150 of these general education courses in science 
was obtained by means of questionnaires and by inspection 
of syllabi, manuals, reading lists, examinations and pub- 
lished articles, especially those in Science in General Educa- 
tion, edited by Earl J. McGrath. In addition, 17 schools 


were visited and information obtained firsthand about 
28 courses. 

The study brought to light some interesting facts con- 
cerning the teaching approach used in general education 
science courses. The subjectmatter survey approach is 
rapidly losing ground and is employed in fewer than one- 
third of the courses offered. In more than two-thirds of 
the courses the teaching approach involves an intensive 
study of selected units or problems or the historical 
development of principles and concepts. 

A number of trends are indicated by the study. One is 
in the direction away from the free election of general 
education science courses and toward the definite prescrip- 
tion of such courses; another is toward the single-subject 
course as opposed to the type in which material from a 
number of sciences is treated in the same course. 

There is no unanimity of opinion on the objectives to 
be achieved in these courses though there appears to be 
increasing emphasis on developing an appreciation of the 
scientific method and understanding of principles and a 
decreasing emphasis on wide coverage of facts.—B.B.W. 





